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1.  Introducción 

Este  es  un  manual  introductorio  al  entorno  de  cálculo  simbólico  Maxima,  directo  su¬ 
cesor  del  legendario  Macsyma. 

El  objetivo  del  manual  es  facilitar  el  acceso  a  este  programa  a  todas  aquellas  personas 
que  por  vez  primera  se  acercan  a  él. 

Maxima  es  un  programa  cuyo  objeto  es  la  realización  de  cálculos  matemáticos,  tanto 
simbólicos  como  numéricos;  es  capaz  de  manipular  expresiones  algebraicas  y  matriciales, 
derivar  e  integrar  funciones,  realizar  diversos  tipos  de  gráficos,  etc. 

Su  nombre  original  fue  Macsyma  ( MAC’s  SYmbolic  MAnipulation  System ,  donde 
MAC,  Machine  Aided  Cognition,  era  el  nombre  del  Laboratory  for  Computer  Science  del 
MIT  durante  la  fase  inicial  del  proyecto  Macsyma).  Se  desarrolló  en  estos  laboratorios  a 
partir  del  año  1969  con  fondos  aportados  por  varias  agencias  gubernamentales  norteame¬ 
ricanas  ( National  Aeronautics  and  Space  Administration,  Office  of  Naval  Research ,  U.S. 
Department  of  Energy  y  U.S.  Air  Forcé). 

El  concepto  y  la  organización  interna  del  programa  están  basados  en  la  tesis  doctoral 
que  Joel  Moses  elaboró  en  el  MIT  sobre  integración  simbólica.  Según  Marvin  Minsky, 
director  de  esta  tesis,  Macsyma  pretendía  automatizar  las  manipulaciones  simbólicas  que 
realizaban  los  matemáticos,  a  fin  de  entender  la  capacidad  de  los  ordenadores  para  actuar 
de  forma  inteligente. 

El  año  1982  es  clave.  El  MIT  transfiere  una  copia  de  Macsyma  a  la  empresa  Symbolics 
Inc.  para  su  explotación  económica,  haciendo  el  código  propietario,  y  otra  al  Departamen¬ 
to  de  Energía,  copia  ésta  que  será  conocida  con  el  nombre  de  DOE-Macsyma.  En  1992  la 
versión  comercial  de  Macsyma  sería  adquirida  por  una  empresa  que  se  llamaría  precisa¬ 
mente  Macsyma  Inc,  y  el  programa  iría  perdiendo  fuelle  progresivamente  ante  la  presencia 
en  el  mercado  de  otros  programas  similares  como  Maple  o  Mathematica,  ambos  los  dos 
inspirados  en  sus  orígenes  en  el  propio  Macsyma. 

Pero  ocurrieron  dos  historias  paralelas.  Desde  el  año  1982,  y  hasta  su  fallecimiento  en  el 
2001,  William  Schelter  de  la  Universidad  de  Texas  mantuvo  una  versión  de  este  programa 
adaptada  al  estándar  Common  Lisp,  la  cual  ya  se  conocía  con  el  nombre  de  Maxima 
para  diferenciarla  de  la  versión  comercial.  En  el  año  1998  Schelter  había  conseguido  del 
DOE  permiso  para  distribuir  Maxima  bajo  la  licencia  GNU-GPL  (www.gnu.org);  con 
este  paso,  muchas  más  personas  empezaron  a  dirigir  su  mirada  hacia  Maxima,  justo  en 
el  momento  en  el  que  la  versión  comercial  estaba  declinando.  Actualmente,  el  proyecto 
está  siendo  mantenido  por  un  grupo  de  desarrolladores  originarios  de  varios  países,  asistidos 
y  ayudados  por  otras  muchas  personas  interesadas  en  Maxima  y  que  mantienen  un  cauce 
de  comunicación  a  través  de  la  lista  de  correo 

http:/ / maxima.sourceforge.net  /  maximalist  .html 

Además,  desde  el  año  2006,  se  ha  activado  una  lista  de  correos  para  usuarios  de  habla 
hispana  en 


http:  /  /  sourceforge.net/mailarchive/forum.phpYforum  úd=50322 
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Puesto  que  Maxima  se  distribuye  bajo  la  licencia  GNU-GPL,  tanto  el  código  fuente 
como  los  manuales  son  de  libre  acceso  a  través  de  la  página  web  del  proyecto, 

http://maxima.sourceforge.net 

Uno  de  los  aspectos  más  relevantes  de  este  programa  es  su  naturaleza  libre;  la  licencia 
GPL  en  la  que  se  distribuye  brinda  al  usuario  ciertas  libertades: 

■  libertad  para  utilizarlo, 

■  libertad  para  modificarlo  y  adaptarlo  a  sus  propias  necesidades, 

■  libertad  para  distribuirlo, 

■  libertad  para  estudiarlo  y  aprender  su  funcionamiento. 

La  gratuidad  del  programa,  junto  con  las  libertades  recién  mencionadas,  hacen  de 
Maxima  una  formidable  herramienta  pedagógica,  de  investigación  y  de  cálculo  técnico, 
accesible  a  todos  los  presupuestos,  tanto  institucionales  como  individuales. 

No  quiero  terminar  esta  breve  introducción  sin  antes  recordar  y  agradecer  a  todo 
el  equipo  humano  que  ha  estado  y  está  detrás  de  este  proyecto,  desde  el  comienzo  de 
su  andadura  allá  por  los  años  70  hasta  el  momento  actual;  incluyo  también  a  todas  las 
personas  que,  aún  no  siendo  desarrolladoras  del  programa,  colaboran  haciendo  uso  de  él, 
presentando  y  aportando  continuamente  propuestas  e  ideas  a  través  de  las  listas  de  correos. 
Un  agradecimiento  especial  a  Robert  Dodier  y  James  Amundson,  quienes  me  brindaron 
la  oportunidad  de  formar  parte  del  equipo  de  desarrollo  de  Maxima  y  me  abrieron  las 
puertas  a  una  experiencia  verdaderamente  enriquecedora. 

Finalmente,  un  recuerdo  a  William  Schelter.  Gracias  a  él  disponemos  hoy  de  un  Maxima 
que,  paradójicamente,  por  ser  libre,  no  pertenece  a  nadie  pero  que  es  de  todos. 


Ferrol-A  Coruña 
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2.  Despegando  con  Maxima 

2.1.  Instalación 

Maxima  funciona  en  Windows,  Linux  y  Mac-OS. 

En  Windows,  la  instalación  consiste  en  descargar  el  binario  exe  desde  el  enlace  corres¬ 
pondiente  en  la  página  del  proyecto  y  ejecutarlo. 

En  Linux,  la  mayoría  de  las  distribuciones  tienen  ficheros  precompilados  en  sus  respec¬ 
tivos  repositorios  con  las  extensiones  rpm  o  deb,  según  el  caso.  Sin  embargo,  en  algunas 
distribuciones  las  versiones  oficiales  no  suelen  estar  muy  actualizadas,  por  lo  que  podría 
ser  recomendable  proceder  con  la  compilación  de  las  fuentes,  tal  como  se  expone  a  conti¬ 
nuación. 

Las  siguientes  indicaciones  hacen  referencia  al  sistema  operativo  Linux. 

Si  se  quiere  instalar  Maxima  en  su  estado  actual  de  desarrollo,  será  necesario  descar¬ 
gar  los  ficheros  fuente  completos  del  CVS  de  Sourceforge  y  proceder  posteriormente  a  su 
compilación.  Se  deberá  tener  operativo  un  entorno  Common  Lisp  en  la  máquina  (clisp, 
cmucl,  sbcl  o  gcl  son  todas  ellas  alternativas  válidas),  así  como  todos  los  programas  que 
permitan  ejecutar  las  instrucciones  que  se  indican  a  continuación,  incluidas  las  dependen¬ 
cias  tcl-tk  y  gnuplot.  Grosso  modo ,  estos  son  los  pasos  a  seguir  para  tener  Maxima 
operativo  en  el  sistema1: 

1 .  Descargar  las  fuentes  a  un  directorio  local  siguiendo  las  instrucciones  que  se  indican 
en  el  enlace  al  CVS  de  la  página  del  proyecto. 

2.  Acceder  a  la  carpeta  local  de  nombre  maxima  y  ejecutar  las  instrucciones 
. /bootstrap 

./configure  — enable-clisp  — enable-lang-es-utf 8 
make 

make  check 
sudo  make  install 
make  clean 

3.  Para  ejecutar  Maxima  desde  la  línea  de  comandos  (Figura  1)  se  deberá  escribir 
maxima 

si  se  quiere  trabajar  con  la  interfaz  gráfica  (Figura  2),  teclear 
xmaxima 

Otro  entorno  gráfico  es  Wxmaxima2  (Figura  3),  que  se  distribuye  conjuntamente  con 
el  ejecutable  para  Windows.  En  Linux,  una  vez  instalado  Maxima,  se  podrá  instalar  este 
entorno  separadamente. 

1Se  supone  que  se  trabaja  con  clisp,  con  codificación  Unicode  y  que  queremos  el  sistema  de  ayuda  en 
inglés  y  español. 

2http:  / /wxmaxima.sourceforge.net 
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Figura  1:  Maxima  desde  la  línea  de  comandos. 


Una  cuarta  alternativa,  que  resulta  muy  vistosa  por  hacer  uso  de  TgX,  es  la  de  ejecutar 
Maxima  desde  el  editor  Texmacs3  (Figura  4) . 

Emacs  es  un  editor  clásico  en  sistemas  Unix.  Aunque  hay  varios  formas  de  ejecutar  Ma¬ 
xima  desde  Emacs,  con  el  modo  Imaxima4  se  consigue  que  los  resultados  estén  formateados 
en  TeX. 

También  se  han  desarrollado  en  los  últimos  años  entornos  web  para  Maxima5. 

3http://www.  texmacs.org 

4http:  / /members3  .j  com.home.ne  .j  p /imaxima 

5http://maxima.sourceforge.net/relatedprojects.shtml 
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2.2.  Una  primera  sesión 

En  esta  sesión  se  intenta  realizar  un  primer  acercamiento  al  programa  a  fin  de  familia¬ 
rizarse  con  el  estilo  operativo  de  Maxima,  dejando  sus  habilidades  matemáticas  para  más 
adelante. 

Una  vez  accedemos  a  Maxima,  lo  primero  que  vamos  a  ver  será  la  siguiente  información: 

Maxima  5.14.0  http://maxima.sourceforge.net 
Using  Lisp  CLISP  2.38  (2006-01-24) 

Distribuí ed  under  the  GNU  Public  License.  See  the  file  COPYING. 

Dedicated  to  the  memory  of  William  Schelter. 

This  is  a  development  versión  of  Maxima.  The  function  bug_report() 
provides  bug  reporting  information. 

(%il) 

En  primer  lugar  nos  informa  sobre  la  versión  de  Maxima  que  se  está  ejecutando  y  la 
URL  del  proyecto,  a  continuación  indica  qué  entorno  Lisp  está  dándole  soporte  al  programa 
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Figura  3:  Wxmaxima,  el  entorno  gráfico  basado  en  wxwidgets. 


y  la  nota  legal  sobre  la  licencia,  por  supuesto  GNU-GPL.  Finaliza  con  un  recordatorio  a 
Schelter  y  cómo  debemos  proceder  si  encontramos  un  fallo  en  el  programa,  ayudando  así  a 
su  mejora.  Tras  esta  cabecera,  muy  importante,  aparece  el  indicador  ( °/0 i  1 )  esperando 
nuestra  primera  pregunta  (i  de  input).  Si  queremos  calcular  una  simple  suma  tecleamos 
la  operación  deseada  seguida  de  un  punto  y  coma  ( ; )  y  una  pulsación  de  la  tecla  retorno 

(°/„il)  45  +  23; 

(%ol)  68 

(*/.  12) 

Vemos  que  el  resultado  de  la  suma  está  etiquetado  con  el  símbolo  (°/0ol)  indicando 
que  es  la  salida  correspondiente  a  la  primera  entrada  (o  de  output).  A  continuación  (%i2) 
indica  que  Maxima  espera  nuestra  segunda  instrucción. 
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El  punto  y  coma  actúa  también  como  un  separador  cuando  escribimos  varias  instruc¬ 
ciones  en  un  mismo  renglón.  Nuestra  siguiente  operación  consistirá  en  asignar  el  valor 
34578  a  la  variable  x  y  el  984003  a  la  y,  solicitando  luego  su  producto, 

(0/„i2)  x: 34578;  y: 984003;  x*y; 

( %o2)  34578 

( %o3)  984003 

( %o4)  34024855734 

conviene  prestar  atención  al  hecho  de  que  la  asignación  de  un  valor  a  una  variable  se  hace 
con  los  dos  puntos  (:),  no  con  el  signo  de  igualdad,  que  se  reserva  para  las  ecuaciones. 

Es  posible  que  no  deseemos  los  resultados  intermedios  que  Maxima  va  calculando  o, 
como  en  este  caso,  las  asignaciones  a  las  variables  que  va  haciendo;  en  tales  situaciones 
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conviene  hacer  uso  del  delimitador  $,  que  no  devuelve  los  resultados  que  va  calculando. 
Repitiendo  el  cálculo  de  la  forma 

(0/„i5)  x: 34578$  y: 984003$  x*y; 

( %o7)  34024855734 

podemos  conseguir  una  salida  más  limpia.  Las  asignaciones  a  variables  se  mantienen  activas 
mientras  dure  la  sesión  con  Maxima,  por  lo  que  podemos  restar  las  variables  x  e  y  anteriores 

(0/.i8)  x-y; 


( %ol3)  -949425 

Esto  tiene  un  riesgo;  si  queremos  resolver  la  ecuación  x2  —  3x  +  1  =  0, 

(°/0i9)  solve(x~2-3*x+l=0,x)  ; 

A  number  was  found  where  a  variable  was  expected  -<solveJ 
—  an  error.  Quitting.  To  debug  this  try  debugmode(true) ; 

nos  devuelve  un  mensaje  de  error,  ya  que  donde  se  supone  que  hay  una  incógnita,  lo 
que  realmente  encuentra  es  un  número,  en  este  caso  34578.  Hay  dos  formas  de  resolver 
esta  situación;  la  primera  consiste  en  utilizar  el  operador  comilla  simple  (’),  que  evita  la 
evaluación  de  las  variables: 

(°/0i9)  solve(,x~2-3*’x+l=0,  ’x) ; 


(%o9) 
(•¿ilO)  x; 


-  _  V5  +  3 

X~  2  ,X~  2 


( %ol0)  34578 

como  se  ve,  x  mantiene  su  valor  original. 

La  segunda  forma  de  evitar  este  problema  consiste  en  vaciar  el  contenido  de  la  variable 
x  mediante  la  función  kill, 

(lili)  kill(x)$ 

(°/0il2)  solve(x~2-3*x+l=0,x)  ; 


(%ol2) 


-  V^"3  _  \/5  +  3 

X~  2  ,X~  2 


(°/.i!3)  x; 
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( %ol3)  x 

en  este  caso,  x  ha  perdido  su  valor. 

Las  etiquetas  (%in)  y  (7,on) ,  siendo  n  un  número  entero,  pueden  ser  utilizadas  a  lo 
largo  de  una  sesión  a  fin  de  evitar  tener  que  volver  a  escribir  expresiones;  así,  si  queremos 
calcular  el  cociente  con  x  =  34578  y  y  =  984003,  podemos  aprovechar  los  resultados 
(°/oo7)  y  C/0O8)  y  hacer 

(7.114)  %o7/°/„o8; 

*  11341618578 

<%°14>  -  316475 

Las  etiquetas  que  indican  las  entradas  y  salidas  pueden  ser  modificadas  a  voluntad 
haciendo  uso  de  las  variables  inchar  y  outchar, 

(°/0il5)  inchar; 

(%ol5)  %i 

(°/0il6)  outchar; 

(%ol6)  %o 

(°/0il7)  inchar:  ’menda; 

( %ol7)  menda 

(mendal8)  outchar:  ’ lerenda; 

(lerendal8)  lerenda 

(mendal9)  2+6; 


(lerendal9)  8 

(menda20)  inchar : ’°/0i$  outchar : ’°/0o$ 

(°/.i22) 

Al  final  se  han  restaurado  los  valores  por  defecto. 

En  caso  de  que  se  pretenda  realizar  nuevamente  un  cálculo  ya  indicado  deberán  escri¬ 
birse  dos  comillas  simples  ’),  no  comilla  doble,  junto  con  la  etiqueta  de  la  instrucción 
deseada 


(°/„i23)  ”7.i8; 
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( %o23)  x  -  984003 

Obsérvese  que  al  haber  dejado  a  x  sin  valor  numérico,  ahora  el  resultado  es  otro. 
Cuando  se  quiera  hacer  referencia  al  último  resultado  calculado  por  Maxima  puede  ser 
más  sencillo  hacer  uso  del  símbolo  de  forma  que  si  queremos  restarle  la  x  a  la  última 
expresión  podemos  hacer 

(%i24)  %-x; 

( %o24)  -984003 

Ciertas  constantes  matemáticas  de  uso  frecuente  tienen  un  símbolo  especial  en  Maxima: 
la  base  de  los  logaritmos  naturales  (e),  el  cociente  entre  la  longitud  de  una  circunferencia 
y  su  diámetro  (7r),  la  unidad  imaginaria  (i  =  y/—l),  la  constante  de  Euler-Mascheroni 
(7  =  —  e~x  In  xdx)  y  la  razón  aurea  (<p  =  1+2V^),  se  representarán  por  °/„e,  °/0pi,  °/„i, 
“/«gamma  y  °/„phi,  respectivamente. 

Es  muy  sencillo  solicitar  a  Maxima  que  nos  informe  sobre  alguna  función  de  su  lenguaje; 
la  técnica  nos  servirá  para  ampliar  información  sobre  cualquier  instrucción  a  la  que  se  haga 
referencia  en  este  manual,  por  ejemplo, 

(°/0i25)  describe  (sqrt) ; 

—  Función:  sqrt  (<x>) 

Raíz  cuadrada  de  <x>.  Se  representa  internamente  por  ‘<x>~(l/2)’. 
Véase  también  ‘rootscontract ’ . 

Si  la  variable  ‘radexpand’  vale  ‘true’  hará  que  las  raíces 
‘n'-ésimas  de  los  factores  de  un  producto  que  sean  potencias  de 
‘n'  sean  extraídas  del  radical;  por  ejemplo,  ‘ sqrt (16*x~2) '  se 
convertirá  en  ‘4*x’  sólo  si  ‘radexpand'  vale  ‘true’. 


There  are  also  some  inexact  matches  for  ‘sqrt’. 

Try  ‘??  sqrt'  to  see  them. 

(°/„o25)  true 

nos  instruye  sobre  la  utilización  de  la  función  sqrt.  Alternativamente,  se  puede  utilizar 
para  el  mismo  fin  el  operador  interrogación  (?), 

C/.Í26)  ?  sqrt 

—  Función:  sqrt  (<x>) 

Raíz  cuadrada  de  <x>.  Se  representa  internamente  por  ‘<x>~(l/2)’. 
Véase  también  ‘rootscontract ’ . 
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Si  la  variable  ‘radexpand’  vale  ‘true’  hará  que  las  raíces 
‘n’-ésimas  de  los  factores  de  un  producto  que  sean  potencias  de 
‘n'  sean  extraídas  del  radical;  por  ejemplo,  ‘ sqrt (16*x~2) '  se 
convertirá  en  ‘4*x’  sólo  si  ‘radexpand'  vale  ‘true’. 


There  are  also  some  inexact  matches  for  ‘sqrt’. 

Try  ‘??  sqrt'  to  see  them. 

(°/0o26)  true 

Cuando  no  conocemos  el  nombre  completo  de  la  función  sobre  la  que  necesitamos 
información  podremos  utilizar  el  operador  de  doble  interrogación  (??)  que  nos  hace  la 
búsqueda  de  funciones  que  contienen  en  su  nombre  cierta  subcadena, 

C/.Í27)  ??  sqrt 

0:  isqrt  (Operadores  generales) 

1:  sqrt  (Operadores  generales) 

2:  sqrtdenest  (Funciones  y  variables  para  simplif ication) 

3:  sqrtdispflag  (Operadores  generales) 

Enter  space-separated  numbers,  ‘all’  or  ‘none’ :  3 

—  Variable  opcional:  sqrtdispflag 
Valor  por  defecto:  ‘true' 

Si  ‘sqrtdispflag'  vale  ‘false’,  hará  que  ‘sqrt'  se  muestre  con  el 
exponente  1/2. 

(°/0o27)  true 

Inicialmente  muestra  un  sencillo  menú  con  las  funciones  que  contienen  la  subcadena 
sqrt  y  en  el  que  debemos  seleccionar  qué  información  concreta  deseamos;  en  esta  ocasión 
se  optó  por  3  que  es  el  número  asociado  a  la  variable  opcional  sqrtdispflag. 

A  fin  de  no  perder  los  resultados  de  una  sesión  con  Maxima,  quizás  con  el  objeto 
de  continuar  el  trabajo  en  próximas  jornadas,  podemos  guardar  en  un  archivo  aquellos 
valores  que  nos  puedan  interesar;  supongamos  que  necesitamos  almacenar  el  contenido  de 
la  variable  y,  tal  como  lo  definimos  en  la  entrada  (°/0i5) ,  así  como  el  resultado  de  la  salida 
(°/0o9),  junto  con  la  instrucción  que  la  generó,  la  entrada  (°/0i9), 

(°/0i28)  save  ( "mi .  sesión" , y ,  resultado=°/0o9 ,  ecuacion=°/0i9)  $ 


véase  que  el  primer  argumento  de  save  es  el  nombre  del  archivo,  junto  con  su  ruta 
si  se  considera  necesario,  que  la  variable  y  se  escribe  tal  cual,  pero  que  las  referencias  de 
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las  entradas  y  salidas  deben  ir  acompañadas  de  un  nombre  que  las  haga  posteriormente 
reconocibles. 

Por  último,  la  forma  correcta  de  abandonar  la  sesión  de  Maxima  es  mediante 
(°/0i29)  quitO; 

Abramos  ahora  una  nueva  sesión  con  Maxima  y  carguemos  en  memoria  los  resultados 
de  la  anterior  sesión, 

C7.il)  load("mi . sesion")$ 

(°/.i2)  y; 


( %o2)  984003 

(°/0i3)  ecuación; 

( %o3)  solve  (x2  —  3  x  +  1  =  0,  a?) 

(°/0i4)  resultado; 


( %o4) 


V5-3 

2 


V5  +  3’ 
2 


Si  ahora  quisiéramos  recalcular  las  soluciones  de  la  ecuación,  modificándola  previa¬ 
mente  de  manera  que  el  coeficiente  de  primer  grado  se  sustituyese  por  otro  número,  sim¬ 
plemente  haríamos 

(°/0i5)  subst  (5,3,  ecuación) ; 


( %o5)  solve  (x2  —  5x  +  1  =  0,x) 

(°/0i6)  ev('/„  nouns)  ; 


(%o6) 


\¡2A  —  5 
*= - 2— * 


V2T+5’ 


La  función  save  guarda  las  expresiones  matemáticas  en  formato  Lisp,  pero  a  nivel 
de  usuario  quizás  sea  más  claro  almacenarlas  en  el  propio  lenguaje  de  Maxima;  si  tal  es 
el  caso,  mejor  será  el  uso  de  stringout,  cuyo  resultado  será  un  fichero  que  también  se 
podrá  leer  mediante  load. 

Continuemos  experimentando  con  esta  segunda  sesión  de  Maxima  que  acabamos  de 
iniciar. 

Maxima  está  escrito  en  Lisp,  por  lo  que  ya  se  puede  intuir  su  potencia  a  la  hora  de 
trabajar  con  listas;  en  el  siguiente  diálogo,  el  texto  escrito  entre  las  marcas  /*  e  */  son 
comentarios  que  no  afectan  a  los  cálculos, 


2  DESPEGANDO  CON  MAXIMA 


14 


(°/0i7)  /*  Se  le  asigna  a  la  variable  x  una  lista  */ 

x:  [eos (%pi) ,  4/16,  [a,  b] ,  (-1)~3,  integrate(u~2,u)] ; 

(%°7) 

(°/0i8)  /*  Se  calcula  el  número  de  elementos  del  último  resultado  */ 
lengthO/o) ; 


( %o8)  5 

(°/0i9)  /*  Se  transforma  una  lista  en  conjunto,  */ 
/*  eliminando  redundancias.  Los  */ 

/*  corchetes  se  transforman  en  llaves  */ 
setify(x) ; 


(%»9)  {-tj.M.y} 

Vemos  a  continuación  cómo  se  procede  a  la  hora  de  hacer  sustituciones  en  una  expre¬ 
sión, 

(°/0ilO)  /*  Sustituciones  a  hacer  en  x  */ 
x,  u=2 ,  a=c; 


(lili)  /*  Forma  alternativa  para  hacer  lo  mismo  */ 
subst([u=2,  a=c] ,x) ; 

(%°n)  [-l.|,[e,b].-l,|] 

(°/0il2)  0/„o7[5],  u=[l,2,3]  ; 

( %0i2)  [ii» 

En  esta  última  entrada,  °/0o7  hace  referencia  al  séptimo  resultado  devuelto  por  Maxima, 
que  es  la  lista  guardada  en  la  variable  x,  de  modo  que  %o7  [5]  es  el  quinto  elemento  de 
esta  lista,  por  lo  que  esta  expresión  es  equivalente  a  a  x  [5] ;  después,  igualando  u  a  la 
lista  [1,2,3]  este  quinto  elemento  de  la  lista  es  sustituido  sucesivamente  por  las  tres 
cantidades.  Tal  como  ya  se  ha  comentado,  Maxima  utiliza  los  dos  puntos  ( : )  para  hacer 
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asignaciones  a  variables,  mientras  que  el  símbolo  de  igualdad  se  reserva  para  la  construcción 
de  ecuaciones. 

Un  aspecto  a  tener  en  cuenta  es  que  el  comportamiento  de  Maxima  está  controlado 
por  los  valores  que  se  le  asignen  a  ciertas  variables  globales  del  sistema.  Una  de  ellas  es 
la  variable  numer,  que  por  defecto  toma  el  valor  lógico  false,  lo  que  indica  que  Maxima 
evitará  dar  resultados  en  formato  decimal,  prefiriendo  expresiones  simbólicas  y  fracciona¬ 
rias, 

(°/0il3)  numer; 

(%ol3) 

C/.Í14)  sqrt  (8)  / 12 ; 

(%ol4) 

(°/0il5)  numer :true$ 

(0/.il6)  sqrt  (8)  / 12 ; 

(%ol6)  0.2357022603955159 

(°/0il7)  numer :false$  /*se  reinstaura  valor  por  defecto  */ 

Sin  embargo,  en  ciertos  contextos  será  preferible  el  uso  de  la  función  f  loat, 

(°/0il8)  float(sqrt(8)/12) ; 

(%ol8)  0.2357022603955159 

Las  matrices  también  tienen  su  lugar  en  Maxima;  la  manera  más  inmediata  de  cons¬ 
truirlas  es  con  la  instrucción  matrix, 

(°/0il9)  A:  matrix(  [sin(x)  ,2,°/0pi]  ,  [0 ,y~2,5/8] ) ; 


(%ol9) 

(°/0i20)  B:  matrix ( [1,2] ,  [0,2] , 


(sin  x  2  7r\ 

o  y 2  ¡) 

[-5 , integrate(x~2 ,x,0, 1)]  ) ; 


false 


2“$ 

~3~ 


(%o20) 


'  1  2' 
0  2 
i 
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(°/0i21)  A.B;  /*  producto  matricial  */ 


(%o21) 


r  +  f  + 


En  cuanto  a  las  capacidades  gráficas  de  Maxima,  se  remite  al  lector  a  la  sección  corres¬ 
pondiente. 

A  fin  de  preparar  publicaciones  científicas,  si  se  quiere  redactar  documentos  en  formato 
UTjvX,  la  función  tex  de  Maxima  será  una  útil  compañera;  en  el  siguiente  ejemplo  se  calcula 
una  integral  y  a  continuación  se  pide  la  expresión  resultante  en  formato  T^X: 


(°/0i22)  ’integrate(sqrt(a+x)/x~5,x,  1,2)  =  intégrate  (sqrt(2+x)/x~5,x,  1,2)$ 


(°/.i23)  tex('/,); 

$$\int_{l}~{2}{{{\sqrt{x+4»\over{x~5»\ ;  dx}={{5\ ,  \sqrt{2}\ ,  \log 
\left  (5-2\ ,  \sqrt{2}\ ,  \sqrt{3}\right)  +548\ ,  \sqrt{3»\over{2048»-{{15 
\ , \sqrt{2>\ , \log  \lef t (3-2\ , \sqrt{2}\right) +244}\over{6144}}$$ 

Al  pegar  y  copiar  el  código  encerrado  entre  los  dobles  símbolos  de  dólar  a  un  documento 
ET^X,  el  resultado  es  el  siguiente: 


I2  x/aH-a 

Ji 


'x  + a  5  y/2  log  (5  —  2  y/2  y/3)  +  548  y/3 


15  y/2  log  (3  —  2  y/2)  +  244 
6144 


Un  comentario  sobre  la  entrada  0/„i22.  Se  observa  en  ella  que  se  ha  escrito  dos  veces 
el  mismo  código  a  ambos  lados  del  signo  de  igualdad.  La  única  diferencia  entre  ambos  es 
el  apóstrofo  que  antecede  al  primer  intégrate;  éste  es  el  operador  de  comilla  simple,  ya 
anteriormente  aparecido,  que  Maxima  utiliza  para  evitar  la  ejecución  de  una  instrucción, 
de  forma  que  en  el  resultado  la  primera  integral  no  se  calcula,  pero  sí  la  segunda,  dando 
lugar  a  la  igualdad  que  se  obtiene  como  resultado  final.  Las  expresiones  que  en  Maxima 
se  marcan  con  la  comilla  simple  para  no  ser  evaluadas  reciben  el  nombre  de  expresiones 
nominales,  en  contraposición  a  las  otras,  conocidas  como  expresiones  verbales. 

Ya  se  ha  comentado  más  arriba  cómo  solicitar  ayuda  de  Maxima.  Algunos  usuarios 
encuentran  engorroso  que  la  ayuda  se  interfiera  con  los  cálculos;  esto  se  puede  arreglar  con 
un  mínimo  de  bricolage  informático.  Al  tiempo  que  se  instala  Maxima,  se  almacena  tam¬ 
bién  el  sistema  de  ayuda  en  formato  html  en  una  determinada  carpeta;  el  usuario  tan  sólo 
tendrá  que  buscar  esta  carpeta,  en  la  que  se  encontrará  la  página  principal  maxima.html, 
y  hacer  uso  de  la  función  system  de  Maxima  que  ejecuta  comandos  del  sistema,  la  cual 
admite  como  argumento  una  cadena  alfanumérica  con  el  nombre  del  navegador  (en  el 
ejemplo,  mozilla)  seguido  de  la  ruta  hacia  el  documento  maxima.html, 


(0/.i24) 

system ( " (mozilla  /usr/local/ share/maxima/5 . 14 . 0/doc/html/ es . utf 8/maxima . html) &" ) $ 
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Este  proceso  se  puede  automatizar  si  la  instrucción  anterior  se  guarda  en  un  fichero 
de  nombre  maxima- init  .mac  y  se  guarda  en  una  carpeta  oculta  de  nombre  .maxima 
en  el  directorio  principal  del  usuario6.  Así,  cada  vez  que  arranque  Maxima,  se  abrirá  el 
navegador  con  el  manual  de  referencia  del  programa. 


'En  Linux 
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3.  Números,  estructuras  de  datos  y  sus  operaciones 

3.1.  Operaciones  aritméticas 

Maxima  puede  trabajar  con  números  enteros  tan  grandes  como  sea  necesario, 
C/.il)  200!; 

(%ol) 

7886578673647905035523632139321850622951359776871732632947425332443594 

4996340334292030428401198462390417721213891963883025764279024263710506 

1926624952829931113462857270763317237396988943922445621451664240254033 

2918641312274282948532775242424075739032403212574055795686602260319041 

7032406235170085879617892222278962370389737472000000000000000000000000 

0000000000000000000000000 

Los  operadores  aritméticos  más  comunes  son: 


+ 

suma 

- 

resta 

* 

producto 

/ 

división 

~  o  ** 

potencia 

sin  aproximaciones  decimales;  así  te- 

3 

haremos 

(0/.i2)  ( (3~7/ (4+3/5)  )  ~  (-i) +7/9)  ~3 ; 

„  620214013952 

^  °°  ’  1307544150375 

obteniendo  el  resultado  en  forma  de  fracción  simplificada. 

De  todos  modos  podemos  solicitar  el  resultado  en  forma  decimal;  por  ejemplo,  si  que¬ 
remos  la  expresión  decimal  del  último  resultado, 

(0/.i3)  floatC/D; 

( %o3)  0.4743350454163436 

Maxima  puede  trabajar  con  precisión  arbitraria.  Para  calcular  el  valor  del  cociente  ^ 
con  100  cifras  decimales,  deberemos  especificar  primero  la  precisión  requerida  asignándole 
a  la  variable  f  pprec  el  valor  100  y  a  continuación  realizar  el  cálculo,  solicitando  la  expresión 
decimal  con  una  llamada  a  la  función  bfloat: 


Maxima  devuelve  los  resultados  de  forma  exacta, 
nemos  que  para  realizar  el  cálculo  de  la  expresión 
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(°/0i4)  fpprec:100$  bfloat(%e  /  °/0pi) ; 

(%o5) 

8.65255979432265087217774789646089617428744623908515539454330288948045 
0445706770586319246625161845173s  x  10"1 

Nótese  que  cuando  un  resultado  es  devuelto  en  el  formato  bf  loat  se  escribe  la  letra  B 
en  lugar  de  la  clásica  E  para  indicar  el  exponente. 

En  la  instrucción  °/„i4  se  establece  la  precisión  deseada  y  luego  se  transforma  el  valor 
simbólico  °/0e  /  °/„pi  mediante  la  función  bfloat.  Recuérdese  que  el  símbolo  $  sirve  como 
delimitador  entre  instrucciones. 

Veamos  cómo  factorizar  un  números  tal  como  200!  en  sus  factores  primos, 

(°/0i6)  factor (°/„ol) ; 

2197  397  549  732  ll19  1316  1711  1910  238  296  316  375  414  434  474  533  593  613  672  712 
( %o6)  732  792  832  892  972  101 103 107 109 113 127 131  137 139 149 151  157 163 167 173 

179181191193197199 

En  la  siguiente  factorización  hacemos  uso  de  la  variable  global  showtime  para  saber  el 
tiempo  que  le  lleva  ejecutar  el  cálculo, 

(°/0i7)  showtime  :true$ 

Evaluation  took  0.00  seconds  (0.00  elapsed)  using  80  bytes. 

(0/„i8)  factor (2~300-l) ; 

Evaluation  took  0.86  seconds  (0.85  elapsed)  using  48.571  MB. 

32  53  7  H  13  31 41 61 101  151  251 331 601 1201 1321  1801 4051  8101 63901 100801 
^/o°  ’  2685011056720113334701  1182468601  1133836730401 

(°/0i9)  showtime  :false$  /*  desactivamos  la  información  de  los  tiempos  */ 

El  texto  que  se  escribe  entre  /*  y  */  son  comentarios  que  Maxima  ignora. 

Además  de  la  variable  opcional  showtime,  otro  mecanismo  de  control  de  los  tiempos 
de  ejecución  es  la  función  time, 

(°/0ilO)  time(°/„o8); 

(%ol0)  [1.7161081] 

En  relación  con  los  números  primos,  para  saber  si  un  número  lo  es  o  no, 

(*/,ill)  primep(3~56-l) ; 

(%oll)  false 

Y  para  solicitarle  a  Maxima  que  compruebe  si  un  número  es  par  o  impar  necesitamos 
echar  mano  de  las  funciones  evenp  or  oddp,  respectivamente, 
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(°/0il2)  evenp(42); 

( %ol2)  true 

(0/.il3)  oddp(31); 

( %ol3)  true 

Le  solicitamos  a  Maxima  que  nos  haga  una  lista  con  todos  los  divisores  de  un  número, 


(°/0il4)  divisors (100) ; 

( %ol4)  {1,2, 4,  5, 10, 20, 25, 50, 100} 

En  lo  que  concierne  a  la  división,  puede  ser  necesario  conocer  el  cociente  entero  y  el 
resto  correspondiente;  para  ello  se  dispone  de  las  funciones  quotient  y  remainder, 

(°/0il5)  quotient  (5689, 214) ; 

( %ol5)  26 

(°/0il6)  remainder (5689, 214)  ; 


3.2.  Números  complejos 

Como  ya  se  comentó  más  arriba,  la  unidad  imaginaria  \f— T  se  representa  en  Maxima 
mediante  el  símbolo  %i, 

(°/.il)  °/„i~2; 

(%ol)  -1 

Se  soportan  las  operaciones  básicas  con  números  complejos, 

(°/0i2)  zl :  3+5*%i$  z2 :  l/2-4*°/0i$  /* zl  y  z2*/ 

(°/0i4)  zl  +  z2;  /*suma*/ 

( %o4)  Ú7- 


(°/0i5)  zl  -  z2; 


/♦resta*/ 
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(%o5) 


(°/0i6)  zl  *  z2;  /^multiplicación*/ 


(%o6)  (| -4¿^)  (5¿  +  3) 

(°/0i7)  zl  /  z2;  /*  división  */ 


(%o25) 


Es  posible  que  estos  dos  últimos  resultados  nos  parezcan  frustrantes  y  los  deseemos  en 
otro  formato;  a  tal  efecto  podemos  pedir  a  Maxima  que  nos  devuelva  (°/„o6)  y  (°/0o7)  en 
forma  cartesiana, 

(°/0i8)  rectf  orm(°/0o6)  ; 


( %o8) 

(°/0i9)  rectf  orm(70o7)  ; 


43 

~2 


19  i 

~Y 


( %o9) 


58  ¿  _  74 
"65"  ~~  65 


Las  funciones  realpart  e  imagpart  extraen  del  número  complejo  sus  partes  real  e 
imaginaria,  respectivamente, 

(°/0ilO)  realpart  (°/0o7) ; 


(%ol0) 

(Xill)  imagpart  (°/0o7) ; 


_74 

65 


( %oll) 


58 

65 


Antes  hemos  optado  por  los  resultados  en  formato  cartesiano;  pero  también  es  posible 
solicitarlos  en  su  forma  polar, 


(°/0il2)  polarf  orm(%o7) ; 


2  \/34  ú  (7r-arctan(  §f ) ) 


( %ol2) 


V65 
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(°/0il3)  rectf orm(°/0)  ;  /*  forma  rectangular  */ 


( %ol3) 


58^34* 


74^34 


(°/0il4)  float(%);  /*  formato  decimal  */ 


( %ol4)  0.8923076923076924  i  -  1.138461538461538 

La  norma  de  un  número  complejo  se  calcula  con  la  función  abs  y  admite  el  argumento 
tanto  en  forma  cartesiana  como  polar, 

(°/0il5)  abs(°/„o9); 


(%ol5) 

(°/0il6)  abs(°/„ol2) ; 


2  VM 
Vfó 


(%ol6) 


2ÚM 

\/65 


Por  último,  el  conjugado  de  un  número  complejo  se  calcula  con  la  función  conjúgate. 
Como  se  ve  en  este  ejemplo,  el  resultado  dependerá  del  formato  del  argumento,  cartesiano 
o  polar, 


(°/0il7)  conjúgate (°/0o9) ;  /*forma  cartesiana*/ 


( %ol7) 


58 ¿  _  74 
"65"  _  65 


(°/0il8)  conjúgate (°/„ol2) ;  /*forma  polar*/ 


(%ol8) 


2  \/34  é  arctan(§?) 
V65 


3.3.  Listas 

Las  listas  son  objetos  muy  potentes  a  la  hora  de  representar  estructuras  de  datos;  de 
hecho,  toda  expresión  de  Maxima  se  representa  internamente  como  una  lista,  lo  que  no 
es  de  extrañar  habida  cuenta  de  que  Maxima  está  programado  en  Lisp  ( List  Processing). 
Veamos  cómo  podemos  ver  la  representación  interna,  esto  es  en  Lisp,  de  una  sencilla 
expresión  tal  como  1  +  3a, 


(°/„il)  :  lisp  #$l+3*a$ 

( (MPLUS  SIMP)  1  ( (MTIMES  SIMP)  3  $a) ) 
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Nótese  que  el  formato  general  es  de  la  forma  :  lisp  #$expr$,  siendo  expr  una  expresión 
cualquiera  en  Maxima. 

Pero  al  nivel  del  usuario  que  no  está  interesado  en  las  interioridades  de  Maxima, 
también  se  puede  trabajar  con  listas  como  las  definidas  a  continuación,  siempre  encerradas 
entre  corchetes, 

(°/0i2)  r :  [1,  [a,  3]  , sqrt(3)/2,  "Don  Quijote"]; 


( %o2)  |^1,  [o,  3] ,  Don  Quijote 

Vemos  que  los  elementos  de  una  lista  pueden  a  su  vez  ser  también  listas,  expresiones 
matemáticas  o  cadenas  de  caracteres  incluidas  entre  comillas  dobles,  lo  que  puede  ser 
aprovechado  para  la  construcción  y  manipulación  de  estructuras  más  o  menos  complejas. 
Extraigamos  a  continuación  alguna  información  de  las  listas  anteriores, 

(°/0i3)  listp(r);  /*  es  r  una  lista?  */ 

( %o3)  true 

(°/0i4)  first(r);  /*  primer  elemento  */ 

(%o4)  1 

(°/0i5)  second(r) ;  /*  segundo  elemento  */ 


( %o5)  [a,  3] 

(°/0i6)  third(r) ;  /*  ...hasta  tenth  */ 


( %o6) 


V3 

2 


(°/0i7)  last(r);  /*  el  último  de  la  lista  */ 


( %o7)  Don  Quijote 

(°/0i8)  rest(r);  /*  todos  menos  el  primero  */ 

( %o8)  |[a,  3] ,  ^ ,  Don  Quijote 

(°/0i9)  part(r,3);  /*  pido  el  que  quiero  */ 
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( %o9)  A 

(°/0ilO)  length(r) ;  /*  cuantos  hay?  */ 

( %ol0)  4 

(Xill)  reverse(r);  /*  le  damos  la  vuelta  */ 

( %oll)  |üon  Quijote,  [a, ; 

(°/0il2)  member(a,r);  /*  es  a  un  elemento?*/ 


( %ol2)  false 

(°/0il3)  member  ( [a, 3]  ,r) ;  /*  lo  es  [a, 3]?  */ 

( %ol3)  true 

(°/0il4)  sort(q:  [7, a, 1 ,d, 5,3,b] ) ;  /*  asigno  valor  a  q  y  ordeno  */ 

(%ol4)  [1,  3,  5,  7,  o,  b,  d\ 

(°/0il5)  delete(  [a, 3]  ,r) ;  /*  borro  elemento  */ 

(%ol5)  |^1,  Don  Quijote 

Nótese  que  en  todo  este  tiempo  la  lista  r  no  se  ha  visto  alterada, 

(0/.il6)  r; 

(%ol6)  |-1,  [o,  3] ,  Don  Quijote 

Algunas  funciones  de  Maxima  permiten  añadir  nuevos  elementos  a  una  lista,  tanto  al 
principio  como  al  final, 

(°/0il7)  cons(l+2,q); 


( %ol7) 

(°/0il8)  endcons(x,q)  ; 


[3, 7, a,  l,d, 5, 3, 6] 
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(%ol8)  [7,  a,l,d,  5,  3,  b,x\ 

(°/.il9)  q; 

(%ol9)  [7,  a,l,d,  5,  3,  b] 

En  este  ejemplo  hemos  observado  también  que  la  lista  q  no  ha  cambiado.  Si  lo  que 
queremos  es  actualizar  su  contenido, 

(°/0i20)  q:  endcons(x,cons(l+2,q))$ 

(°/.i21)  q; 

(%o21)  [3, 7,  a,  l,d,  5, 3,  b,  x] 

Es  posible  unir  dos  listas, 

(°/0i22)  append(r,q); 

( %o22)  |l,  [a,  3] ,  ^,Don  Quijote,  3, 7,  a,  1,  d,  5, 3,  b,  xj 

Cuando  los  elementos  de  una  lista  van  a  obedecer  un  cierto  criterio  de  construcción, 
podemos  utilizar  la  función  makelist  para  construirla, 

(°/0i23)  s : makelist (2+k*2,k, 0, 10)  ; 

( %o23)  [2, 4,  6,  8, 10, 12, 14, 16, 18, 20, 22] 

donde  le  hemos  indicado  a  Maxima  que  nos  construya  una  lista  con  elementos  de  la  forma 
2+2*k,  de  modo  que  k  tome  valores  enteros  de  0  a  10. 

La  función  apply  permite  suministrar  a  otra  función  todos  los  elementos  de  una  lista 
como  argumentos,  así  podemos  sumar  o  multiplicar  todos  los  elementos  de  la  lista  s  recién 
creada, 

(0/„i24)  apply ("+" ,s) ; 

( %o24)  132 

(°/0i25)  apply ("*"  ,s)  ; 

( %o25)  81749606400 

aunque  estas  dos  operaciones  hubiese  sido  quizás  más  apropiado  haberlas  realizado  con 
las  funciones  sum  y  product. 

A  veces  interesará  aplicar  una  misma  función  a  varios  elementos  de  una  lista  de  forma 
independiente,  para  lo  que  haremos  uso  de  map;  a  continuación  un  ejemplo  de  cálculo  de 
factoriales  y  otro  trigonométrico, 
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(°/0i26)  map(" ! " ,s) ; 

[2, 24, 720, 40320, 3628800, 479001600, 87178291200, 20922789888000, 

^  /o°  '  6402373705728000,2432902008176640000, 1124000727777607680000] 

(°/0i27)  map(sin, s) ; 

( %o27)  [sin  2,  sin  4,  sin  6,  sin  8,  sin  10,  sin  12,  sin  14,  sin  16,  sin  18,  sin  20,  sin  22] 

A  veces  se  quiere  utilizar  en  map  una  función  que  no  está  definida  y  que  no  va  a 
ser  utilizada  más  que  en  esta  llamada  a  map.  En  tales  casos  se  puede  hacer  uso  de  las 
funciones  lambda-,  supóngase  que  a  cada  elemento  x  de  una  lista  se  le  quiere  aplicar  la 
función  f(x)  =  sin(x')  +  ^,  sin  necesidad  de  haberla  definido  previamente, 

(°/0i28)  map(lambda(  [x]  ,  sin(x)  +  1/2) ,  [7, 3,4,9] ) ; 


( %o28)  |sin7  +  sin  3  +  ^,sin4  +  ^,sin9  +  ^ 

Como  se  ve,  la  función  lambda  admite  dos  argumentos;  el  primero  es  una  lista  (no  se 
puede  evitar  la  redundancia)  de  argumentos,  siendo  el  segundo  la  expresión  que  indica 
qué  se  hace  con  los  elementos  de  la  lista  anterior. 

Por  último,  las  listas  también  pueden  ser  utilizadas  en  operaciones  aritméticas, 

C/.Í29)  [l,2,3]  +  [a,b,c]  ; 

(%o29)  [a+l,6  +  2,c+3] 

C/.Í30)  [l,2,3]*[a,b,c]  ; 


( %o30)  [a,  2  b,  3  c] 

(°/0i31)  [1,2,3]/  [a,b, c]  ; 


(%o31) 

(°/0i32)  [1,2,3]-  [a,b, c]  ; 


1  2  3' 

a’b’c 


( %o32)  [1  -  o,  2  -  b,  3  -  c] 

(°/0i33)  [1,2,3]  .  [a,b,c]  ;  /*  producto  escalar  */ 


( %o33) 


3c  +  26  +  a 
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(°/0i34)  [a,b,c]~3; 

( %o34)  [a3,  b3,  c3] 

(0/.i35)  3“[a,b,c]; 

(%o35)  [3a,3b,3c] 

En  Maxima,  las  listas  se  pueden  interpretar  como  vectores;  para  más  información, 
acódase  a  la  Sección  5.6 

Para  que  estas  operaciones  puedan  realizarse  sin  problemas,  la  variable  global  listarith 
debe  tomar  el  valor  true,  en  caso  contrario  el  resultado  será  bien  distinto, 

(°/0i36)  listarith  :false$ 

(0/0i37)  [1,2, 3] +  [4,5, 6]  ; 

(%o37)  [4,5,6]  +  [1,2,3] 

(°/0i38)  listarith  :true$ 

Como  ya  se  vió  al  comienzo  de  esta  sección,  una  lista  puede  ser  elemento  de  otra  lista, 
si  queremos  deshacer  todas  las  listas  interiores  para  que  sus  elementos  pasen  a  formar 
parte  de  la  exterior, 

(°/0i39)  flatten(  [1,  [a,b]  ,2,3,  [c,  [d,e]]] ) ; 

(%o39)  [1,  o,  b,  2, 3,  c,  d,  e] 

El  cálculo  del  módulo  de  un  vector  se  puede  hacer  mediante  la  definición  previa  de  una 
función  al  efecto: 

(°/0i40)  modulo  (v)  :  = 

if  listp(v) 

then  sqrt (apply ("+" ,v~2) ) 

else  error ("Mucho  ojito:  ",  v,  "  no  es  un  vector  !!!!")$ 

(°/0i41)  xx :  [a,b,  c  ,d,  e]  $ 

(°/„i42)  yy :  [3 , 4 ,  -6 , 0 , 4/5]  $ 

(°/0i43)  modulo  (xx-yy) ; 


( %o43) 


e  —  - 


+  cP  +  (c  +  6)2  +  (b  —  4) 2  +  (a  —  3)' 
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3.4.  Arrays 

Los  arrays  son  estructuras  que  almacenan  datos  secuencial  o  matricialmente.  Se  di¬ 
ferencian  de  las  listas  en  que  éstas  pueden  alterar  su  estructura  interna  y  los  arrays  no; 
pero  la  ventaja  frente  a  las  listas  estriba  en  que  el  acceso  a  un  dato  cualquiera  es  directo, 
mientras  que  en  las  listas,  al  ser  estructuras  dinámicas,  para  acceder  al  n-ésimo  dato,  un 
puntero  debe  recorrer  todos  los  elementos  de  la  lista  desde  el  que  ocupa  el  primer  lugar. 
Así  pues,  se  recomienda  el  uso  de  arrays  cuando  se  vaya  a  utilizar  una  secuencia  larga 
de  cantidades  cuya  estructura  se  va  a  mantener  inalterable;  en  general,  se  mejorarán  los 
tiempos  de  procesamiento. 

Si  queremos  construir  un  array,  lo  debemos  declarar  previamente  mediante  la  función 
array,  indicando  su  nombre  y  dimensiones;  a  continuación  le  asignaríamos  valores  a  sus 
elementos, 

array (arreglillo ,2,2); 

arreglillo 

arreglillo [0,0]  :  34$ 
arreglillo [1 ,2] :x+y$ 
arreglillo  [2,2]  :  f  loat  (°/0pi)  $ 
listarray (arreglillo) ; 

[34,  #####,  #####,  #####,  #####, 

V  +  x,  #####,  #####,  3.141592653589793] 

En  el  ejemplo,  hemos  definido  un  array  de  dimensiones  3  por  3,  ya  que  la  indexación  de 
los  arrays  comienza  por  el  cero,  no  por  el  uno  como  en  el  caso  de  las  listas;  a  continuación 
hemos  asignado  valores  a  tres  elementos  del  array:  un  entero,  una  expresión  simbólica  y 
un  número  decimal  en  coma  flotante.  Finalmente,  la  función  listarray  nos  permite  ver 
el  contenido  del  array,  utilizando  el  símbolo  #####  para  aquellos  elementos  a  los  que  aún 
no  se  les  ha  asignado  valor  alguno. 

Si  ya  se  sabe  de  antemano  que  el  array  va  a  contener  únicamente  números  enteros  o 
números  decimales,  podemos  declararlo  al  definirlo  para  que  Maxima  optimice  el  procesa¬ 
miento  de  los  datos.  A  continuación  se  declara  un  array  de  números  enteros  haciendo  uso 
del  símbolo  f  ixnum  (en  el  caso  de  números  decimales  en  coma  flotante  se  utilizaría  f  lonum) 
y  a  continuación  la  función  fillarray  rellenará  el  array  con  los  elementos  previamente 
almacenados  en  una  lista, 

(°/.i6)  lis :  [6,0,-4,456,56, -99] ; 


(°/.il) 


(%ol) 

(°/.i2) 

(°/.i3) 

(°/.i4) 

(0/.i5) 


(%o5) 


[6,0,-4,456,56,-99] 


(°/0i7)  array(arr ,f ixnum, 5) ; 
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( %o7)  arr 

(°/0i8)  f illarray (arr , lis) ; 


(0/„i9)  arr  [4]  ; 


( %o9)  56 

(°/0ilO)  listarray(arr) ; 

( %olO)  [6,0,-4,456,56,-99] 

Hasta  ahora  hemos  hecho  uso  de  los  llamados  arr  ay  declarados.  Si  no  hacemos  la 
declaración  pertinente  mediante  la  función  array,  y  en  su  lugar  realizamos  asignaciones 
directas  a  los  elementos  de  la  estructura,  creamos  lo  que  se  conoce  como  array  de  claves 
( hash  array )  o,  en  la  terminología  de  Maxima,  array  no  declarado;  estos  arrays  son  mucho 
más  generales  que  los  declarados,  puesto  que  pueden  ir  creciendo  de  forma  dinámica  y  los 
índices  no  tienen  que  ser  necesariamente  números  enteros, 

(°/„ill)  tab[9]  :76; 

(%oll)  76 

(°/0il2)  tab[-58]  :2/9; 

( %0l2)  l 

(°/0il3)  tab  [juan]  :  sqrt  (8) ; 

( %ol3)  2  § 

(°/0il4)  tab  ["cadena"]  :  a-b; 

( %ol4)  a  —  b 

(°/0il5)  tab [x+y]  :5; 


( %ol5) 

(°/0il6)  listarray(tab) ; 


5 
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(%ol6)  ^,76,a-6,2Í,jj 

(°/0il7)  tab ["cadena"]  ; 

( %ol7)  a  -  b 

Esta  generalidad  en  la  estructura  de  los  arrays  no  declarados  tiene  la  contrapartida  de 
restarle  eficiencia  en  el  procesamiento  de  los  datos  almacenados,  siendo  preferible,  siempre 
que  sea  posible,  la  utilización  de  los  declarados. 

Ya  por  último,  junto  con  los  arrays  declarados  y  no  declarados  mencionados,  es  posible 
trabajar  a  nivel  de  Maxima  con  arrays  de  lisp.  En  este  caso  su  construcción  se  materializa 
con  la  función  make_array, 

(°/0il8)  a:  make_array(flonum,5) ; 

(%ol8)  {Array  :  #  (0.0  0.0  0.0  0.0  0.0)} 

(0/.il9)  a[0]  :3.4; 

(%ol9)  3.4 

(°/0i20)  a  [4]  :  6 . 89 ; 

( %o20)  6.89 

(°/.i21)  a; 

(%o21)  {Array  :  #  (3.4  0.0  0.0  0.0  6.89)} 

Nótense  aquí  un  par  de  diferencias  respecto  de  la  función  array:  el  número  5  indica  el 
número  de  elementos  (indexados  de  0  a  4)  y  no  el  índice  máximo  admisible;  por  otro  lado, 
la  función  no  devuelve  el  nombre  del  objeto  array,  sino  el  propio  objeto. 

3.5.  Cadenas 

Una  cadena  se  construye  escribiendo  un  texto  o  cualquier  otra  expresión  entre  comillas 
dobles.  Vamos  a  hacer  algunos  ejemplos  básicos, 

c°/.il)  xx :  "Los 


(%ol) 

(°/0i2)  yy:  "  cerditos"; 


Los 
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( %o2)  cerditos 

(°/0i3)  zz :  concat  (xx ,  sqrt  (9)  ,  yy) ; 

( %o3)  Los  3  cerditos 

la  función  concat  admite  varios  parámetros,  los  cuales  evalúa  y  luego  concatena  convir¬ 
tiéndolos  en  una  cadena. 

La  función  stringp  nos  indica  si  su  argumento  es  o  no  una  cadena, 

(0/„i4)  stringp  (zz) ; 

( %o4)  true 

Para  saber  el  número  de  caracteres  de  una  cadena  haremos  uso  de  slength 
(°/0i5)  slength  (zz) ; 

( %o5)  14 

Una  vez  transformada  en  cadena  una  expresión  de  Maxima,  deja  de  ser  evaluable, 
(%i6)  concat (sqrt (25) ) ; 

( %o6)  5 

(°/.i7)  */.  +  2; 

( %o7)  5  +  2 

(°/0i8)  stringp  (°/0) ; 

( %o8)  false 

el  resultado  no  se  simplifica  porque  uno  de  los  sumandos  es  una  cadena  y  no  se  reconoce 
como  cantidad  numérica. 

La  función  charat  devuelve  el  carácter  que  ocupa  una  posición  determinada, 

(°/0i9)  charat  (zz, 7)  ; 

( %o9)  c 

También  podemos  solicitar  que  Maxima  nos  devuelva  una  lista  con  todos  los  caracteres 
o  palabras  que  forman  una  cadena, 
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(°/0ilO)  charlist  (zz)  ; 

(%olO)  [L,o,s,  ,3,  .  c,e,r,d,  i.t.  o.s] 

(Xill)  tokens(zz); 

(%oll)  [Los,  3,cerditos] 

O  podemos  transformar  letras  minúsculas  en  mayúsculas  y  viceversa,  todas  las  letras 
o  sólo  un  rango  de  ellas, 

(°/0il2)  /*  todo  a  mayúsculas  */ 
supcase (zz) ; 


( %ol2)  LOS  3  CERDITOS 

(°/0il3)  /*  a  minúsculas  entre  la  9a  y  12a  */ 
sdowncase (% , 9 , 12) ; 


( %ol3)  LOS  3  CErdiTOS 

(°/0il4)  /*  a  minúsculas  de  la  13a  hasta  el  final  */ 
sdowncase (% , 13) ; 


( %ol4)  LOS  3  CErdiTos 

(°/0il5)  /*  todo  a  minúsculas  */ 
sdowncase (%) ; 


( %o72)  los  3  cerditos 

(°/0il6)  /*  a  mayúsculas  entre  la  la  y  2a  */ 
supcase (7o,  1 ,2) ; 

( %o74)  Los  3  cerditos 

Para  comparar  si  dos  cadenas  son  iguales  disponemos  de  la  función  sequal, 
(70il7)  tt :  zz; 


(%ol7) 

(70il8)  sequal  (tt,  zz) ; 


Los  3  cerditos 
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( %  0I8)  true 

(°/0il9)  tt:  supcase(zz); 

( %ol9)  LOS  3  CERDITOS 

(°/0i20)  sequal(tt ,zz) ; 

( %o20)  false 

Podemos  buscar  subcadenas  dentro  de  otras  cadenas  mediante  la  función  ssearch;  los 
siguientes  ejemplos  pretenden  aclarar  su  uso, 

(0/.i21)  zz; 

( %o21)  Los  3  cerditos 

(°/0i22)  /*  busca  la  subcadena  "cer"  */ 
ssearch("cer" ,zz) ; 


(%o22)  7 

(°/0i23)  /*  busca  la  subcadena  "Cer"  */ 
ssearch("Cer" ,zz) ; 


( %o23)  false 

(°/0i24)  /*  busca  "Cer"  ignorando  tamaño  */ 
ssearch("Cer" ,zz, ’ sequalignore) ; 

(%o24)  7 

(°/0i25)  /*  busca  sólo  entre  caracteres  1  y  5  */ 
ssearch ("Cer" ,zz, ’ sequalignore ,1,5) ; 

( %o25)  false 

Igual  que  buscamos  subcadenas  dentro  de  otra  cadena,  también  podemos  hacer  susti¬ 
tuciones;  en  este  caso  la  función  a  utilizar  es  ssubst, 

(°/0i26)  ssubst  ( "mosqueteros " , "  cerditos " ,  zz) ; 
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( %o26)  Los  3  mosqueteros 

Esta  función  tiene  como  mínimo  tres  argumentos,  el  primero  indica  la  nueva  cadena 
a  introducir,  el  segundo  es  la  subcadena  a  ser  sustituida  y  el  tercero  la  cadena  en  la  que 
hay  que  realizar  la  sustitución.  Es  posible  utilizar  esta  función  con  más  argumentos;  como 
siempre,  para  más  información,  ?  ssubst. 

Con  la  función  substring  podemos  obtener  una  subcadena  a  partir  de  otra,  para  lo 
cual  debemos  indicar  los  extremos  de  la  subcadena, 

(°/0i27)  substring(zz ,7) ; 

( %o27)  cerditos 

(°/0i28)  substring(zz,3,8) ; 

( %o28)  s  3  c 

Como  se  ve,  el  segundo  argumento  indica  la  posición  a  partir  de  la  cual  queremos  la 
subcadena  y  el  tercero  la  última  posición;  si  no  se  incluye  el  tercer  argumento,  se  entiende 
que  la  subcadena  se  extiende  hasta  el  final  de  la  de  referencia. 

En  algunos  contextos  el  usuario  tiene  almacenada  una  expresión  sintácticamente  correc¬ 
ta  de  Maxima  en  formato  de  cadena  y  pretende  evaluarla,  en  cuyo  caso  la  función  eval_string 
la  analizará  y  evaluará,  en  tanto  que  parse_string  la  analiza  pero  no  la  evalúa,  limitándo¬ 
se  a  transformarla  a  una  expresión  nominal  de  Maxima, 

(°/0i29)  /*  una  cadena  con  una  expresión  de  Maxima  */ 
expr:  " intégrate (x~5 ,x) " ; 

(%o29)  intégrate  (xA5,x) 

(°/0i30)  /*  se  evalúa  la  cadena  */ 
eval_string(expr) ; 

(%030)  y 

(°/0i31)  /*  se  transforma  en  una  expresión  nominal  */ 
parse_string(expr) ; 

(%o31)  intégrate  (x5 ,  x) 

C/.Í32) 

(%o32)  y 

Véase  que  el  resultado  de  parse_string  ya  no  es  una  cadena,  sino  una  expresión 
correcta  de  Maxima  sin  evaluar,  cuya  evaluación  se  solicita  a  continuación  con  la  doble 
comilla  simple. 
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3.6.  Conjuntos 

Se  define  a  continuación  un  conjunto  mediante  la  función  set, 

C7.il)  el :  set  (a,  [2  ,k]  ,b, sqrt  (2) , a,  set  (a,b)  , 

3, "Sancho" ,set() ,b,sqrt(2) ,a) ; 

( %ol)  |3,  \/2,  {}  ,  [2,  k]  ,  Sancho,  a,  {a,  b}  ,  6  j 

Como  se  ve,  se  admiten  objetos  de  muy  diversa  naturaleza  como  elementos  de  un 
conjunto:  números,  expresiones,  el  conjunto  vacío  ({}),  listas,  otros  conjuntos  o  cadenas  de 
caracteres.  Cuando  se  trabaja  con  listas,  puede  ser  de  utilidad  considerar  sus  componentes 
como  elementos  de  un  conjunto,  luego  se  necesita  una  función  que  nos  transforme  una  lista 
en  conjunto, 

(°/oi2)  [ [2 , k]  ,sqrt(2)  ,set(b,a) ,  [k,2]  ,  "Panza"]  ; 

( %o2)  |j2,  k] ,  \/2,  {a,  b}  ,  [ k ,  2] ,  PanzaJ 

(°/0i3)  c2 :  setify  (°/0)  ; 

( %o3)  {"\/2,  [2,  k]  ,  Panza,  {o,  b}  ,  [ k ,  2]  j 

el  cambio  en  la  naturaleza  de  estas  dos  colecciones  de  objetos  se  aprecia  en  la  presencia  de 
llaves  en  lugar  de  corchetes.  De  igual  manera,  podemos  transformar  un  conjunto  en  lista, 

(°/0i4)  listify (70ol)  ; 

( %o4)  |^3,  \/2,  {}  ,  [2,  k]  ,  Sancho,  a,  {a,  b}  ,  bj 

Comprobemos  de  paso  que  {}  representa  al  conjunto  vacío, 

(0/.i5)  emptyp  ('/.  [3]  )  ; 

( %o5)  true 

Recuérdese  que  °/0  sustituye  a  la  última  respuesta  dada  por  Maxima,  que  en  este  caso 
había  sido  una  lista,  por  lo  que  7.  [3]  hace  referencia  a  su  tercera  componente. 

Para  comprobar  si  un  cierto  objeto  forma  parte  de  un  conjunto  hacemos  uso  de  la 
instrucción  elementp, 

(70i6)  elementp  (sqrt  (2) ,  el)  ; 

( %o6)  true 

Es  posible  extraer  un  elemento  de  un  conjunto  y  luego  añadirle  otro  distinto 


3  NÚMEROS,  ESTRUCTURAS  DE  DATOS  Y  SUS  OPERACIONES 


36 


(°/0i7)  el:  disjoin(sqrt (2) , el) ;  /*  sqrt(2)  fuera  */ 

( %o7)  {3,  {}  ,  [2,  k ] ,  Sancho,  a,  {a,  b}  ,  6} 

(°/0i8)  el:  adjoin(sqrt(3) ,cl) ;  /*  sqrt(3)  dentro  */ 

( %o8)  |3,  V3,  {}  ,  [2,  k]  ,  Sancho,  a,  {a,  b}  ,  6  j 

La  sustitución  que  se  acaba  de  realizar  se  pudo  haber  hecho  con  la  función  subst, 
(°/0i9)  /*  nuevamente  a  poner  sqrt(2)  */ 
subst(sqrt(2) , sqrt(3) ,  el) ; 

( %o9)  |3,  V2 ,  {}  ,  [2,  k]  ,  Sancho,  a,  {a,  b}  ,  6  j 

La  comprobación  de  si  un  conjunto  es  subconjunto  de  otro  se  hace  con  la  fución 
subsetp, 

(°/0ilO)  subsetp  (set  ( [k,2]  ,  "Panza")  ,c2) ; 

( %ol0)  true 

A  continuación  algunos  ejemplos  de  operaciones  con  conjuntos, 

(lili)  union(cl , c2) ; 

( %oll)  |3,  \/2,  Vs,  {}  ,  [2,  k] ,  Panza,  Sancho,  o,  {a,  b}  ,  b,  [ k ,  2]  j 

(°/0il2)  intersection(cl,c2)  ; 

(%ol2)  {[2  ,k},{a,b}} 

(°/0il3)  setdif  f  erence  (el ,  c2)  ; 

(%ol3)  |3,  V3,  {}  ,  Sancho,  a,  &j 

(°/0il4)  cardinality  00  ; 

(%ol4)  6 

Vemos  aquí  también  cómo  pedir  el  cardinal  de  un  conjunto. 

Igual  que  se  ha  visto  cómo  aplicar  una  función  a  todos  los  elementos  de  una  lista, 
podemos  hacer  lo  mismo  con  los  de  un  conjunto, 

(°/0il5)  map (sin, set (1,2, 3, 4, 5)) ; 

( %ol5)  {sin  1 ,  sin  2,  sin  3,  sin  4,  sin  5} 

Por  último  ya,  el  producto  cartesiano  de  tres  conjuntos, 

(°/0il6)  cartesian_product(set(l,2)  ,set(a,b,c)  ,set(x,y)) ; 


(%ol6) 


{[1,  a,  x]  ,  [1,  o,  y] ,  [1,  b,  x ] ,  [1,6,  y] ,  [1,  c,  x\ ,  [1,  c,  y]  , 
[2,  a,  x] ,  [2,  a,  y] ,  [2,  b,  x]  ,  [2,  b,  y] ,  [2,  c,  z] ,  [2,  c,  y]} 
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3.7.  Grafos 

El  soporte  para  el  desarrollo  de  algoritmos  en  el  ámbito  de  la  teoría  de  grafos  se 
encuentra  en  el  módulo  graphs.  Se  trata  de  un  paquete  de  reciente  creación  que  aporta  la 
infraestructura  básica  para  el  trabajo  con  grafos,  orientados  o  no.  No  se  admiten  bucles 
ni  aristas  múltiples. 

C/.il)  load  (graphs)  $ 

(°/0i2)  g:  create_graph( 

[1,2, 3, 4,5, 6, 7] , 

[  [1,2]  ,[2, 3]  ,[2, 4],  [3, 4], 

[4, 5],  [5, 6],  [4, 6],  [6, 7]  ]); 

( %o2)  GRAPH 

Como  vemos,  la  respuesta  de  la  función  create_graph  no  es  muy  informativa,  pero 
nos  podemos  hacer  una  idea  de  la  estructura  del  grafo  haciendo 

(°/.i3)  print_graph(g)  $ 

Graph  on  7  vértices  with  8  edges . 

Adjacencies: 

7  :  6 

6:745 
5:64 
4  :  6  5  3  2 

3:42 
2:431 
1  :  2 

Vemos  que,  por  ejemplo,  el  vértice  número  3  comparte  aristas  (no  orientadas)  con  los 
vértices  4  y  2. 

En  caso  de  que  el  grafo  anterior  fuese  orientado  (digrafo)  su  definición  sería: 

(°/0i4)  dg:  create_graph( 

[1,2,3,4,5,6,71 , 

[  [1,2],  [2, 3],  [2, 4],  [3, 4], 

[4, 5],  [5, 6],  [4, 6],  [6, 7]  ], 
directed  =  true)$ 

(0/.i5)  print_graph(dg)  $ 

Digraph  on  7  vértices  with  8  ares. 

Adjacencies : 

7  : 

6  :  7 
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Figura  5:  Grafos:  a)  no  orientado;  b)  orientado;  c)  complemetario  del  no  orientado;  d) 
camino  más  corto. 


5  :  6 

4:65 
3  :  4 

2:43 
1  :  2 

En  este  caso,  print_graph  nos  muestra  los  nodos  hijos  a  la  derecha  de  los  padres.  En 
particular,  el  nodo  7  no  es  origen  de  ningún  arco,  pero  es  hijo  del  nodo  6.  Otra  manera 
alternativa  de  mostrar  la  estructura  de  un  grafo  es  mediante  la  función  draw_graph 

(°/0i6)  draw_graph(g,  show_id=true ,  terminal=eps)$ 

(°/0i7)  draw_graph(dg,  show_id=true ,  head_length=0 . 05 ,  terminal=eps)$ 

Cuyos  aspectos  se  pueden  observar  en  los  apartados  a)  y  b)  de  la  Figura  5. 
Continuemos  con  el  grafo  no  orientado  g  definido  más  arriba  y  veamos  cómo  extraer  de 
él  cierta  información:  el  número  cromático,  matriz  de  adyacencia,  grafo  complementario, 
si  es  conexo,  o  planar,  su  diámetro,  etc. 

(°/0i8)  chromatic_number  (g) ; 


(%o8) 


3 
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(°/0i9)  adjacency_matrix(g) ; 

/O  1  O  O  o  o  0\ 

10  110  0  0 
0  10  10  0  0 
( %o9)  0  110  110 

0  0  0  1  0  1  0 

0  0  0  1  1  0  1 

\0  O  O  O  O  1  0/ 

(°/0ilO)  ge :  complement_graph(g) ; 

( %olO)  GRAPH 

(lili)  print_graph(°/„)  $ 

Graph  on  7  vértices  with  13  edges . 

Adjacencies : 

1  :  3  4  5  6  7 

2:567 

3  :  15  6  7 

4:17 

5  :  12  3  7 

6:123 

7  :  1  2  3  4  5 

Pedimos  ahora  a  Maxima  que  nos  represente  gráficamente  el  complementario  de  g, 
cuyo  aspecto  es  del  apartado  c)  de  la  Figura  5. 

(°/0il2)  draw_graph(gc,  show_id=true,  terminal=eps)$ 

Continuamos  con  algunas  otras  propiedades: 

(°/0il3)  is_connected(g) ; 

( %ol3)  true 


(°/.il4)  is_planar  (g)  ; 


(%ol4) 


true 
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(°/0il5)  diaineter  (g)  ; 

( %ol5)  4 

Junto  con  la  definición  manual,  también  es  posible  generar  determinados  gratos  me¬ 
diante  funciones,  como  el  cíclico,  el  completo,  el  dodecaedro,  un  grafo  aleatorio,  etc. 

(°/0il6)  print_graph(gl :  cycle_graph(5)  )$ 

Graph  on  5  vértices  with  5  edges . 

Adjacencies: 

4:03 

3:42 

2:31 

1:20 

0:41 


(°/0il7)  print_graph(g2:  complete_graph(4))$ 

Graph  on  4  vértices  with  6  edges . 
Adjacencies : 

3:210 

2:310 

1:320 

0:321 


(°/0il8)  print_graph(g3 :  dodecahedron_graph()  )$ 

Graph  on  20  vértices  with  30  edges. 

Adjacencies: 

19  :  14  13  1 

18  :  13  12  2 

17  :  12  11  3 

16  :  11  10  4 

15  :  14  10  0 

14  :  5  19  15 

13  :  7  19  18 

12  :  9  18  17 

11  :  8  17  16 

10  :  6  16  15 

9  :  7  8  12 

8  :  9  6  11 

7  :  5  9  13 
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6  :  8  5  10 

5  :  7  6  14 

4  :  16  0  3 

3  :  17  4  2 

2  :  18  3  1 

1  :  19  2  0 

0  :  15  4  1 


(°/0il9)  /*  0.4  de  probabilidad  para  cada  arco  */ 
print_graph (g4 :  random_graph (10,0.4))$ 

Graph  on  10  vértices  with  15  edges. 

Adjacencies: 

9  :  8  5  2  1  0 

8:963 

7:521 

6:83 

5:97 

4  :  2 

3:860 
2  :  9  7  4  0 

1:970 
0  :  9  3  2  1 

Por  último  ya,  un  ejemplo  de  cálculo  de  la  ruta  más  corta.  Se  genera  un  grafo  orientado 
aleatorio,  a  continuación  se  calcula  la  ruta  más  corta  entre  dos  vértices  dados  y  finalmente 
se  hace  una  representación  gráfica  del  cálculo  realizado,  la  cual  se  puede  ver  en  el  apartado 
d)  de  la  Figura  5. 

(°/0i20)  r  :  random_digraph(20,  0.25)$ 


(°/0i21)  /*  Una  lista  de  los  vértices  del  digrafo  */ 
vértices (r) ; 

( %o21)  [19, 18, 17, 16, 15, 14, 13, 12, 11, 10, 9, 8, 7, 6, 5, 4, 3, 2, 1, 0] 


(°/0i22)  sh:  shortest_path(0 , 14, r) ; 


( %o22) 


[0,4,1,17,14] 
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(°/0i23)  draw_graph( 
r , 

show_edges 

show_id 

head_length 

terminal 


vertices_to_path(sh) , 
true , 

0.05, 

eps)$ 


El  paquete  graphs  también  permite  el  uso  de  etiquetas  en  los  nodos  y  ponderaciones  en 
las  aristas  y  arcos.  El  lector  interesado  podrá  consultar  la  documentación  correspondiente7. 


rA  partir  de  la  versión  5.14  de  Maxima 
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4.  Álgebra 

4.1.  Transformaciones  simbólicas 

Sin  duda,  una  de  las  capacidades  más  destacables  de  Maxima  es  su  habilidad  para 
manipular  expresiones  algebraicas.  Desarrollemos  un  ejemplo  que  empieza  por  asignar  a 
la  variable  q  una  expresión  literal: 

(°/.il)  q:  (x+3)~5-(x-a)~3+(x+b)~(-l)  +  (x-l/4)~(-5)  ; 

(%ol)  1  -(ifc_o)3  +  (l  +  3)S  +  _L_ 

Se  observa  que  en  principio  Maxima  no  realiza  ningún  cálculo.  La  función  expand  se 
encarga  de  desarrollar  las  potencias, 

(°/0i2)  expand (q) ; 


- .4  ,  3  \  a  , - r-  +  -TT  +  X5  +  15  X4  +  89  X3  +  3  a  x2  +  270  x2- 

( %o2)  -6-^  +  ^-^  +  Ü-uk  *+b 

3  a2  x  +  405  x  +  a3  +  243 

No  obstante  es  posible  que  no  nos  interese  desplegar  toda  la  expresión,  entre  otras 
cosas  para  evitar  una  respuesta  farragosa  y  difícil  de  interpretar;  en  tal  caso  podemos 
utilizar  expand  añadiéndole  dos  argumentos  y  operar  de  la  manera  siguiente 

(°/0i3)  q,expand(3,2) ; 


( %o3) 


-  +  (x  +  3)5  —  x3  +  3  a  x2  —  3  a2  x  + 


(*-if 


Con  el  primer  argumento  indicamos  que  queremos  la  expansión  de  todas  aquellas 
potencias  con  exponente  positivo  menor  o  igual  a  3  y  de  las  que  teniendo  el  exponente 
negativo  no  excedan  en  valor  absoluto  de  2. 

Dada  una  expresión  con  valores  literales,  podemos  desear  sustituir  alguna  letra  por 
otra  expresión;  por  ejemplo,  si  queremos  hacer  los  cambios  a  =  2,  b  =  2c  en  el  último 
resultado  obtenido, 


(°/0i4)  °/0,a=2,b=2*c; 


( %o4)  — 3_  +  (x  +  3)5-x3  +  6x2-12x+ - Vt  +  8 

x  +  ¿c  (.x-  -  4) 

En  estos  dos  últimos  ejemplos  (°/0i3  y  °/0i4)  se  presentaron  sentencias  en  las  que  había 
elementos  separados  por  comas  (,).  Esta  es  una  forma  simplificada  de  utilizar  la  función 
ev,  que  evalúa  la  primera  expresión  asignando  los  valores  que  se  le  van  indicando  a  con¬ 
tinuación;  por  ejemplo,  (%i3)  se  podía  haber  escrito  de  la  forma  ev(q,expand(3,2))  o 
expand(q,3,2)  y  (°/0i4)  como  ev (°/0 , a=2 , b=2*c) .  El  uso  de  la  variante  con  ev  está  más 
indicado  para  ser  utilizado  dentro  de  expresiones  más  amplias.  Obsérvese  el  resultado 
siguiente 
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(°/0i5)  3*x~2  +  ev(x~4,x=5); 


(%o5)  3x2  +  625 

donde  la  sustitución  x  =  5  se  ha  realizado  exclusivamente  dentro  del  entorno  delimitado 
por  la  función  ev. 

De  forma  más  general,  la  función  subst  sustituye  subexpresiones  enteras.  En  el  siguien¬ 
te  ejemplo,  introducimos  una  expresión  algebraica  y  a  continuación  sustituimos  todos  los 
binomios  x+y  por  la  letra  k, 

0/.i6)  1/ (x+y)-(y+x)/z+(x+y)~2; 


( %o6) 

(°/0i7)  subst  (k ,  x+y ,  °/„) ; 


-{y  +  x) 


1 

y  +  x 


( %o44)  _^  +  A;2+l 

No  obstante,  el  siguiente  resultado  nos  sugiere  que  debemos  ser  precavidos  con  el  uso 
de  esta  función,  ya  que  Maxima  no  siempre  interpretará  como  subexpresión  aquella  que 
para  nosotros  sí  lo  es: 

(°/0i8)  subst  (sqrt(k)  ,x+y,  (x+y)  ~2+(x+y) ) ; 

( %o8)  y  +  x  +  k 

Como  una  aplicación  práctica  de  subst,  veamos  cómo  podemos  utilizarla  para  obtener 
el  conjugado  de  un  número  complejo, 

(°/0i9)  subst  (-°/„i ,  '/.i ,  a+b*°/0i)  ; 


( %o9)  a  —  ib 

La  operación  inversa  de  la  expansión  es  la  factorización.  Expandamos  y  factoricemos 
sucesivamente  un  polinomio  para  comprobar  los  resultados, 

(°/0ilO)  expandí  (a-2)*(b+l)~2*(a+b)~5) ; 

(%ol0) 

a  b7  -  2  b7  +  5  a2  b6  -  8  a  b6  -  4  b6  +  10  a3  65  -  10  a2  b5  -  19  a  b5  -  2  b5  +  10  a4  b4- 
35  a2  b4  -  10  a  b4  +  5  o5  b3  +  10  o4  b 3  -  30  a3  b3  -  20  a2  b3  +  a6  b2  +  8  a5  b2  -  10  o4  b2  - 
20  a3  b2  +  2  o6  b+  a5  b-  10  a4  b  +  o6  —  2  a5 


(Xill)  factor (%) ; 
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(%oll)  (a -2)  (b  +  lf  (b  +  a)5 

El  máximo  común  divisor  de  un  conjunto  de  polinomios  se  calcula  con  la  función  gcd 
y  el  mínimo  común  múltiplo  con  lcm 

(0/„il2)  pl:  x~7-4*x~6-7*x~5+8*x~4+ll*x~3-4*x''2-5*x; 

( %ol2)  x7  -  4  x6  -  7  x5  +  8  x4  +  11  x3  -  4  x2  -  5  x 

(0/„il3)  p2 :  x~4-2*x~3-4*x~2+2*x+3; 

(%ol3)  x4  —  2x3  —  4x2  +  2x  +  3 

(°/.il4)  gcd(pl ,p2) ; 

( %ol4)  x3  +  x2  -  x  -  1 

(°/0il5)  load(f\mcts)$ 

(°/0il6)  lcm(pl  ,p2)  ; 

( %ol6)  (x  —  5)  (x  —  3)  (x  —  l)2  x  (x  +  l)3 

En  (%i  12)  y  (°/0il3)  definimos  los  polinomios  pl  y  p2,  a  continuación  calculamos  su 
máximo  común  divisor  (mcd)  en  (°/0il4)  y  antes  de  pedir  el  mínimo  común  múltiplo  en 
(°/0il6)  cargamos  el  paquete  functs  en  el  que  se  encuentra  definida  la  función  lcm.  Es 
posible  que  deseemos  disponer  del  mcd  factorizado,  por  lo  que  hacemos 

(°/0il7)  factor (0/0ol4)  ; 

(%ol7)  (x-l)(:r  +  l)2 

Si  en  algún  momento  queremos  inhibir  las  simplificaciones,  podemos  hacer  uso  de  la 
variable  global  simp,  cuyo  valor  por  defecto  es  true, 

(°/0il8)  simp:  false$ 

(0/„il9)  2  +  2  +  a  +  a; 

( %ol9)  2  +  2  +  o  +  o 

(°/0i20)  simp:  true$ 

(°/0i21)  7„ol9; 

(%o21)  2o +  4 

Veamos  ahora  algunas  transformaciones  que  se  pueden  hacer  con  fracciones  algebraicas. 
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(°/0i22)  expr :  (x“2-x) / (x~2+x-6) -5/ (x~2-4) ; 


( %o22) 

Podemos  factorizar  esta  expresión. 
(°/0i23)  factor  (expr)  ; 


x  —  x 
x2  +  x  —  6 


5 

x 2  —  4 


( %o23) 


(x  —  3)  (x2  +  4x  +  5) 
(x  -  2)  (x  +  2)  (x  +  3) 
La  expansión  devuelve  el  siguiente  resultado, 

(°/0i24)  expand(expr)  ; 


( %o24) 


5 

'  ^234 


x2  +  x  —  6  x2  +  x  —  6 
Ahora  descomponemos  la  fracción  algebraica  expr  en  fracciones  simples, 


(°/0i25)  partf  rae  (expr,  x) ; 


(%o25) 


12  5  17  1 

5  (x  +  3)  +  4  (x  +  2)  20  (x  -  2)  + 


Por  último,  transformamos  la  misma  expresión  a  su  forma  canónica  CRE  ( Cannonical 
Rational  Expression ),  que  es  un  formato  que  Maxima  utiliza  para  reducir  expresiones 
algebraicas  equivalentes  a  una  forma  única, 


(°/0i26)  radcan(expr)  ; 


(%o26) 


x3  +  x2  —  7  x  —  15 
x3  +  3x2  —  4x  —  12 


La  función  ratsimp  también  simplifica  cualquier  expresión  racional,  así  como  las  sub¬ 
expresiones  racionales  que  son  argumentos  de  funciones  cualesquiera.  El  resultado  se  de¬ 
vuelve  como  el  cociente  de  dos  polinomios.  En  ocasiones  no  es  suficiente  con  una  sola 
ejecución  de  ratsimp,  por  lo  que  será  necesario  aplicarla  más  veces,  esto  es  lo  que  hace 
precisamente  la  función  fullratsimp;  concretemos  esto  con  un  ejemplo: 


(°/„i27)  (x~(a/2)-l)~2*(x~(a/2)+l)~2  /  (x~a-l) ; 


(%o27) 


xa  —  1 
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(°/0i28)  ratsimp(yo) ;  /*  simplificamos  una  vez  */ 


(%o28) 


x2a  -2xa  +  1 
xa  —  1 


(°/0i29)  ratsimp(7o) ;  /*  simplificamos  otra  vez  */ 


( %o29)  Xa  -  1 

(°/0i30)  fullratsimp(°/„o27) ;  /*  simplificamos  todo  de  una  vez!  */ 


( %o30)  xa  -  1 

Dada  una  fracción  algebraica,  podemos  obtener  separadamente  el  numerador  y  el  de¬ 
nominador, 

(°/„i31)  fr:  (x~3-4*x~2+4*x-2)/(x~2+x+l) ; 

( %o31) 

(°/0i32)  num(fr) ; 

( %o32)  x3  -4x2  +  4x-2 

(°/0i33)  denom(fr) ; 

( %o33)  x2  +  x  +  1 

Maxima  nunca  puede  adivinar  a  priori  las  intenciones  del  usuario  para  con  las  expre¬ 
siones  con  las  que  está  trabajando,  por  ello  sigue  la  política  de  no  tomar  ciertas  decisiones, 
por  lo  que  a  veces  puede  parecer  que  no  hace  lo  suficiente  o  que  no  sabe  hacerlo.  Esta  forma 
vaga  de  proceder  la  suple  con  una  serie  de  funciones  y  variables  globales  que  permitirán 
al  usuario  darle  al  motor  simbólico  ciertas  directrices  sobre  qué  debe  hacer  con  las  ex¬ 
presiones,  cómo  reducirlas  o  transformarlas.  Sin  ánimo  de  ser  exhaustivos,  siguen  algunos 
ejemplos  en  los  que  se  controlan  transformaciones  de  tipo  logarítmico  y  trigonométrico. 

Las  expresiones  que  contienen  logaritmos  suelen  ser  ambiguas  a  la  hora  de  reducirlas. 
Transformamos  a  continuación  la  misma  expresión  según  diferentes  criterios, 

(°/0i34)  log(x~r)-log(x*y)  +  a*log(x/y) ; 


x3  —  4  x2  +  4x  —  2 

X2  +  X  +  1 


( %o34)  —  log  (xy)  +  a  log  +  r  l°g  x 

En  principio,  ha  hecho  la  transformación  logxr  — *  r  log  x .  Esto  se  controla  con  la 
variable  global  logexpand,  cuyo  valor  por  defecto  es  true,  lo  cual  sólo  permite  la  reducción 
recién  indicada.  Si  también  queremos  que  nos  expanda  los  logaritmos  de  productos  y 
divisiones,  le  cambiaremos  el  valor  a  esta  variable  global,  tal  como  se  indica  a  continuación, 
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(°/0i35)  logexpand:  super$ 

(0/.i36)  log(x~r)-log(x*y)  +  a*log(x/y) ; 

( %o36)  —  log  y  +  a  (log  x  —  log  y)  +  r  log  x  —  log  x 

(°/0i37)  /*  pedimos  que  nos  simplifique  la  expresión  anterior  */ 
ratsimp(7o) ; 

( %o37)  (—a  —  1)  log  y  +  (r  +  a  —  log  x 

En  caso  de  no  querer  que  nos  haga  transformación  alguna  sobre  la  expresión  de  entrada, 
simplemente  haremos 

(°/0i38)  logexpand:  false$ 

(0/.i39)  log(x~r)-log(x*y)  +  a*log(x/y) ; 

( %o39)  —  log  (xy)  +  a  log  +  log  xr 

Devolvemos  ahora  la  variable  logexpand  a  su  estado  por  defecto, 

(°/0i40)  logexpand:  true$ 

La  función  logcontract  es  la  que  nos  va  a  permitir  compactar  expresiones  logarítmi¬ 
cas, 

(°/0i41)  logcontract  (2* (a*log(x)  +  2*a*log(y) ) ) ; 

( %o41)  o  log  ( x 2  y4) 

Por  defecto,  Maxima  contrae  los  coeficientes  enteros;  para  incluir  en  este  caso  la  variable 
a  dentro  del  logaritmo  le  asignaremos  la  propiedad  de  ser  un  número  entero, 

(°/0i42)  declare  (a,  integer)$ 

(°/0i43)  ,,0/0i41;  /*  fuerza  reevaluación  de  expresión  nominal  */ 

( %o43)  log  (x2ay4a) 

Esta  es  una  forma  de  hacerlo,  pero  Maxima  permite  un  mayor  refinamiento  sobre  este 
particular  haciendo  uso  de  la  variable  global  logconcoef  fp. 

Toda  esta  casuística  de  diferentes  formas  de  representar  una  misma  expresión  aumen¬ 
ta  considerablemente  con  aquéllas  en  las  que  intervienen  las  funciones  trigonométricas 
e  hiperbólicas.  Empecemos  con  la  función  trigexpand,  cuyo  comportamiento  se  con¬ 
trola  con  las  variables  globales  trigexpand  (que  le  es  homónima),  trigexpandplus  y 
trigexpandtimes,  cuyos  valores  por  defectos  son,  respectivamente,  false,  true  y  true. 
Puestos  a  experimentar,  definamos  primero  una  expresión  trigonométrica  e  interpretamos 
algunos  resultados, 
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(°/0i44)  expr:  x+sin(3*x+y)/sin(x) ; 


( %o44) 


sin(¡,  +  3x)+x 
sin  x 


(°/0i45)  trigexpand(expr) ;  /*  aquí  trigexpand  vale  false  */ 


(%o45) 


eos  (3  x)  sin  y  +  sin  (3  x)  eos  y 

- : - 1-  x 

smx 


Vemos  que  sólo  se  desarrolló  la  suma  del  numerador;  ahora  cambiamos  el  valor  lógico 
de  trigexpand, 


(°/0i46)  trigexpand(expr) ,  trigexpand=true ; 


( %o46) 


(eos3  x  —  3  eos  x  sin2  x)  sin  y  +  (3  eos2  x  sin  x  —  sin3  x)  eos  y 

S1D.X 


+  x 


Cuando  la  asignación  de  la  variable  global  se  hace  como  se  acaba  de  indicar,  sólo 
tiene  efecto  temporal  durante  esa  ejecución,  sin  haberse  alterado  a  nivel  global.  Podemos 
ordenar  a  Maxima  que  simplifique  la  expresión  anterior, 


(°/0i47)  ratsimp(70); 


eos  x  sin  x  —  eos0  x 


( %o47) 

Podemos  inhibir  el  desarrollo  de  la  suma, 

(°/0i48)  trigexpand (expr) ,  trigexpandplus=f  alse ; 


3  x  —  3  eos2  x  s 


y  —  x  sin  x 


( %o48) 


sinfa  +  3  *)+x 
sin  x 


Otra  variable  global  de  interés,  que  no  tiene  nada  que  ver  con  la  función  trigexpand, 
es  halfangles,  con  valor  false  por  defecto,  que  controla  si  se  reducen  los  argumentos 
trigonométricos  con  denominador  dos, 


(°/0i49)  sin(x/2) ; 


( %o49) 


sin(!) 


(°/0i50)  sin(x/2) ,  halfangles=true; 
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( %o50) 


a/1  —  cosa: 

V2 


La  función  trigsimp  fuerza  el  uso  de  las  identidades  fundamentales  sin2  x  +  eos2  x  =  1 
y  cosh2  x  —  sinh2  x  =  1  para  simplificar  expresiones, 

(°/0i51)  5*sin(x)~2  +  4*cos(x) ~2; 


( %o51)  5  sin2  x  +  4  eos2  x 

(0/.i52)  trigsimp  (°/„)  ; 


( %o52)  sin2  x  +  4 

Como  un  último  ejemplo,  veamos  cómo  reducir  productos  y  potencias  de  funciones 
trigonométricas  a  combinaciones  lineales, 

(°/0i53)  -sin(x)  ~2+3*cos  (x)  ~2*tan(x)  ; 

( %o53)  3  eos2  x  tan  x  —  sin2  x 

(°/0i54)  trigreduceC/J  ; 


( %o54) 


3  sin  (2  x)  t  eos  (2  x)  1 

2  +  2  2 


Y  si  queremos  transformar  expresiones  trigonométricas  en  complejas, 
(°/0i55)  exponentialize(0/0) ; 


( %o55) 


4 


3 i  (e2ix  -  e~2ix)  1 
4  2 


Otras  funciones  útiles  en  este  contexto  y  que  el  usuario  podrá  consultar  en  la  docu¬ 
mentación  son:  triginverses,  trigsign  y  trigrat. 
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4.2.  Ecuaciones 

Almacenar  una  ecuación  en  una  variable  es  tan  simple  como  hacer 
C/.il)  ec:3*x=l  +  x; 

(%ol)  3x  =  x  +  l 

A  partir  de  ahí,  aquellas  operaciones  en  las  que  intervenga  la  variable  serán  realizadas 
a  ambos  miembros  de  la  igualdad;  restamos  x  en  los  dos  lados  y  a  continuación  dividimos 
lo  que  queda  entre  2, 

C/.Í2)  %-x; 


( %o2)  2x  1 

(°/.i3)  %/2; 

( %o3)  *=1 

obteniendo  de  esta  manera  la  solución  de  la  ecuación  como  si  hubiésemos  operado  ma¬ 
nualmente. 

Ni  qué  decir  tiene  que  la  ecuación  anterior  se  pudo  haber  resuelto  de  un  modo  más 
inmediato, 

(%i4)  solve(ec); 


(%o4) 

La  instrucción  solve  puede  admitir  como  segundo  argumento  la  incógnita  que  se  pre¬ 
tende  calcular,  lo  que  resultará  de  utilidad  cuando  en  la  ecuación  aparezcan  constantes 
literales, 

(°/0i5)  solve  ( (2-a)  /x-3=b*x+l/x ,  x) ; 

i/(4  —  4a)  6  +  9  +  3  ^(4  -  4o)  b  +  9  -  3 

2b  2b 

Las  soluciones  de  las  ecuaciones  serán  probablemente  utilizadas  en  cálculos  posteriores, 
por  lo  que  nos  interesará  poder  extraerlas  de  la  lista  anterior;  a  continuación  tomamos  el 
primer  resultado  calculado  por  Maxima  mediante  la  función  part  y  después  asignamos  a 
la  variable  sol  el  resultado  numérico, 


(0/.i6)  part  C/0, 1) ; 
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( %o6) 

(°/0i7)  sol:  rhs(°/0); 


V(4-4o)  6  +  9  +  3 
2  b 


%  _  \/(4-4a)  b  +  9  +  3 

2  b 

La  función  rhs  devuelve  el  miembro  derecho  de  la  igualdad,  mientras  que  lhs  haría  lo 
propio  con  el  miembro  izquierdo. 

Es  posible  resolver  ecuaciones  polinómicas  de  grado  <  4,  pero  desgraciadamente,  como 
es  de  esperar,  Maxima  no  dispone  de  un  método  general  algebraico  que  permita  resolver 
ecuaciones  polinómicas  de  grado  mayor  que  cuatro, 

(°/0i8)  solve (x~5  -  3*x~4  +  2*x~3  -2*x~2  -  x  +  4  =  0) ; 


( %o8)  [0  =  x5 -3x4  +  2x3 -2x2 -x  +  4] 

por  lo  que  solve  devolverá  la  misma  ecuación  sin  resolver. 

También  solve  se  puede  utilizar  para  la  resolución  de  sistemas,  en  cuyo  caso  las  ecua¬ 
ciones  deben  ir  agrupadas  en  una  lista,  así  como  las  incógnitas;  nos  planteamos  la  resolución 
del  siguiente  sistema  no  lineal  de  incógnitas  x  e  y, 

í  3  *x2-y2  =  6 
\  x  =  y  +  a 

(°/0i9)  solve (  [3*x~2-y~2=6,x=y+a]  ,  [x,y]) ; 


( %o9) 


[[x  =  =  _V3Vm+3  «' 

^  VaA+i—a  ^  y  _•  ■y/3-y/aTf4— 3aj  j' 


Una  alternativa  al  uso  de  solve  es  la  función  algsys.  Veamos  cómo  algsys  trata  la 
resolución  de  la  ecuación  polinómica  anterior  °/0o8, 

(°/0ilO)  algsys  ([x~5  -  3*x~4  +  2*x~3  -2*x~2  -  x  +  4  =  0]  ,  [x] ) ; 


[[x  =  2.478283086356668] ,  [x  =  .1150057557117295  -  1.27155810694299  i] , 

( %ol0)  [x  =  1.27155810694299  i  +  .1150057557117295] ,  [x  =  -.8598396689940523] , 

[x  =  1.151545166402536]] 

Como  se  ve,  al  no  ser  capaz  de  resolverla  algebraicamente,  nos  brinda  la  oportunidad 
de  conocer  una  aproximación  numérica  de  la  solución.  La  función  algsys  reconoce  la 
variable  global  realonly,  que  cuando  toma  el  valor  true,  hará  que  se  ignoren  las  soluciones 
complejas, 
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(lili)  realonly :true$ 

(°/0il2)  ,,0/oilO;  /*  recalcula  la  entrada  °/0ilO  */ 

(%ol2) 

[[x  =  2.478283086356668] ,  [x  =  1.151545166402536] ,  [x  =  -.8598396689940523]] 
(°/0il3)  realonly :false$  /*  le  devolvemos  el  valor  por  defecto  */ 

Un  sistema  no  lineal  con  coeficientes  paramétricos  y  complejos 

Su  —  av  =  t 

18  =  3»+« 

|  =  1 

(°/0il4)  algsys(  [3*u-a*v=t,  (2+yoi)/(u+t)=3*v+u,t/u=l]  ,  [u,v,t]); 


(%ol4) 


Veamos  cómo  trata  algsys  las  ecuaciones  indeterminadas,  devolviendo  la  solución  en 
términos  paramétricos, 

(°/0il5)  algsys ( [3*x~2-5*y=x]  ,  [x,y]); 


v= 

u  y/2  ,U  y/2y/Ey/a+6,‘ 


_  \fa 
\/2  yj  a+6  5 


_  \!  a+6  +  a+6  _  2  y/í+2  .  _  y/i+2y/a 

y/2  ’  y/2y/ay/a+6  ’  y/2  y/ a+6 


(%ol5) 

(0/„il6)  */.,  %rlíi; 


'  %ruy=- 


2  -  %n 


(%ol6)  [[.T  -  |]] 

Maxima  nombra  los  parámetros  siguiendo  el  esquema  %rn,  siendo  n  un  número  entero 
positivo.  En  la  entrada  °/0il6  pedimos  que  en  °/„ol5  se  sustituya  el  parámetro  por  la  unidad. 

En  ocasiones  conviene  eliminar  las  variables  de  los  sistemas  de  forma  controlada;  en 
tales  casos  la  función  a  utilizar  será  elimínate, 

(°/0il7)  elimínate (  [3*u-a*v=t,  (2+°/0i)/(u+t)=3*v+u,t/u=l]  ,  [u] ) ; 


(%ol7) 


[2í  —  a  v,  —  (a2  +  9  o)  v2  —  (5  a  +  36)  t  v  —  4 12  +  9  i  +  18] 
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Cuando  de  lo  que  se  trata  es  de  resolver  un  sistema  de  ecuaciones  lineales,  la  mejor 
opción  es  linsolve,  cuya  sintaxis  es  similar  a  la  de  las  funciones  anteriores.  En  el  siguiente 
ejemplo,  el  objetivo  es 

2x  -  4  y  +  2z 

\x  +  2  y  +  9z 

-4x  +  y/2y  +  2 

(°/0il8)  linsolve  ( 

[2*x-4*y+2*z=-2, 

1/3*  x+2*y+9*z=x+y, 

-4  *  x  +  sqrt (2)  *y+z=3*y], 

[x,y,z] ) ; 

3021  y/2  -  12405  _  1537  y/2 +  16642  _ _  53  y/2  -  27681 

48457  ,V  ~  48457  48457  J 

Cuando  la  matriz  de  coeficientes  del  sistema  es  dispersa,  la  función  f  ast_linsolve 
será  preferible,  ya  que  aprovechará  tal  circunstancia  para  encontrar  las  soluciones  de  forma 
más  rápida. 

Si  los  coeficientes  del  sistema  se  tienen  en  formato  matricial,  quizás  sea  más  apropiado 
el  uso  de  la  función  linsolve_by_lu,  tal  como  se  indica  en  la  sección  dedicada  a  las 
matrices. 

No  todo  es  resoluble  simbólicamente.  Existen  en  Maxima  varios  procedimientos  cuya 
naturaleza  es  estrictamente  numérica.  Uno  de  ellos  es  realroots,  especializado  en  el  cálcu¬ 
lo  de  aproximaciones  racionales  de  las  raíces  reales  de  ecuaciones  polinómicas;  el  segundo 
parámetro,  que  es  opcional,  indica  la  cota  de  error. 

(°/0il9)  realroots (x~8+x~3+x+l ,  5e-6) ; 


(%ol8) 


(%ol9) 


371267' 
524288 _ 


En  cambio,  allroots  obtiene  aproximaciones  en  formato  decimal  de  coma  flotante  de 
todas  las  raíces  de  las  ecuaciones  polinómicas,  tanto  reales  como  complejas, 


(°/0i20)  allroots  (x~8+x~3+x+l) ; 


[x  =  .9098297401801199  i  +  .2989522918873167,  , 
x  =  .2989522918873167  -  .9098297401801199  i  , 

(V  20)  x=  --7081337759784606>  *  =  -9807253637807569 i  -  .4581925662678885  , 

^  /o°  ’  x-  —.9807253637807569 *  -  .4581925662678885,  , 

x  =  .5359278244124014  i  +  1.013307162369803  , 

x  =  1.013307162369803  -  .5359278244124014  i, x  =  -1.000000000000001] 

Más  allá  de  las  ecuaciones  algebraicas,  f  ind_root  utiliza  el  método  de  bipartición  para 
resolver  ecuaciones  en  su  más  amplio  sentido, 
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(°/0i21)  f(x):=144  *  sin(x)  +  12*sqrt (3) *°/„pi  -  36*x~2  -  12*°/„pi*x$ 

(°/0i22)  find_root(f (z) ,z,l,2) ; 

(%o22)  1.899889962840263 

El  uso  del  método  de  Newton  requiere  cargar  en  memoria  el  módulo  correspondiente. 
Veamos  como  ejemplo  la  ecuación 


2uu  -  5  =  u 

(°/0i23)  load(mnewton)$  /*  carga  el  paquete  */ 

(°/0i24)  mnewton(  [2*u~u-5]  ,  [u]  ,  [1] )  ; 

0  errors,  0  warnings 

( %ol4)  [[u  =  1.70927556786144]] 

y  el  sistema 

í  x  +  3  log(r)  -  y2  =  0 
\  2x 2  —  xy  —  5x  +  1  =  0 

(°/0i25)  mnewton(  [x+3*log(x)-y~2,  2*x~2-x*y-5*x+l]  ,  [x,  y],  [5,  5]); 

( %o25)  [[x  =  3.756834008012769,  y  =  2.779849592817897]] 

(°/0i26)  mnewton(  [x+3*log(x)-y"'2,  2*x~2-x*y-5*x+l]  ,  [x,  y],  [1,  -2]); 

(%o26)  [[x  =  1.373478353409809,  y  =  -1.524964836379522]] 

En  los  anteriores  ejemplos,  el  primer  argumento  es  una  lista  con  la  ecuación  o  ecuaciones 
a  resolver,  las  cuales  se  suponen  igualadas  a  cero;  el  segundo  argumento  es  la  lista  de 
variables  y  el  último  el  valor  inicial  a  partir  del  cual  se  generará  el  algoritmo  y  cuya 
elección  determinará  las  diferentes  soluciones  del  problema. 

4.3.  Matrices 

La  definición  de  una  matriz  es  extremadamente  simple  en  Maxima, 

(7„il)  mi:  matrix(  [3,4,0]  ,  [6 ,0,-2]  ,  [0,6, a] ) ; 

/3  4  0  \ 

( %ol)  6  0  -2 

\0  6  a  ) 

También  es  posible  definir  una  matriz  de  forma  interactiva  tal  y  como  muestra  el 
siguiente  ejemplo, 
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(°/0i2)  entermatrix(2 ,3) ; 
Row  1  Column  1 : 

4/7; 

Row  1  Column  2: 

0; 

Row  1  Column  3: 

°/.pi; 

Row  2  Column  1 : 
sqrt (2) ; 

Row  2  Column  2: 
log(3) ; 

Row  2  Column  3: 

-9; 

Matrix  entered. 


(%°2)  (\/2  log  3  -9) 

Existe  un  tercer  método  para  construir  matrices  que  es  útil  cuando  el  elemento  (i,  j)- 
ésimo  de  la  misma  es  función  de  su  posición  dentro  de  la  matriz.  A  continuación,  se  fija 
en  primer  lugar  la  regla  que  permite  definir  un  elemento  cualquiera  y  luego  en  base  a  ella 
se  construye  una  matriz  de  dimensiones  2x5 

(0/„i3)  a[i,j]  :=i+j$ 

(°/0i4)  genmatrix(a,2,5) ; 


(%o4) 


2  3  4  5  6 
V3  4  5  6  7/ 


Obsérvese  que  el  símbolo  de  asignación  para  el  elemento  genérico  es  :=. 

Podemos  acceder  a  los  diferentes  elementos  de  la  matriz  haciendo  referencia  a  sus 
subíndices,  indicando  primero  la  fila  y  después  la  columna: 

(7.15)  mi  [3,1]; 


( %o5)  0 

Se  puede  extraer  una  submatriz  con  la  función  submatrix,  teniendo  en  cuenta  que  los 
enteros  que  preceden  al  nombre  de  la  matriz  original  son  las  filas  a  eliminar  y  los  que  se 
colocan  detrás  indican  las  columnas  que  no  interesan;  en  el  siguiente  ejemplo,  queremos  la 
submatriz  que  nos  queda  de  mi  después  de  extraer  la  primera  fila  y  la  segunda  columna, 


(°/0i6)  submatrix (1  ,ml  ,2) ; 
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<*->  (o  ;2) 

Otro  ejemplo  es  el  siguiente, 

(°/0i7)  submatrix(l,2,ml,3) ; 

( %o7)  (0  6) 

en  el  que  eliminamos  las  dos  primeras  filas  y  la  última  columna,  ¿se  pilla  el  truco? 

Al  igual  que  se  extraen  submatrices,  es  posible  añadir  filas  y  columnas  a  una  matriz 
dada;  por  ejemplo, 

(°/0i8)  addrow(ml ,  [1,1,1]  ,  [2,2,2]  ) ; 


(3  4 
6  0 
0  6 
1  1 
\2  2 


¡J 


C/.i9)  addcoia,  [7,7,7,7,7]  )  ; 


( %o9) 


/3  4  0  7\ 

6  0-27 
0  6  a  7 
1117 
\2  2  2  7/ 


La  matriz  identidad  es  más  fácil  construirla  mediante  la  función  ident, 


(°/0ilO)  ident  (3); 


(1  0  0\ 

( %ol0)  0  10 

\0  0  l) 

y  la  matriz  con  todos  sus  elementos  iguales  a  cero, 
C7.il!)  zeromatrix(2,4) ; 


/ (w  f o  o  o  0\ 

(%0ll)  (o  0  0  o) 

También,  una  matriz  diagonal  con  todos  los  elementos  de  la  diagonal  principal  iguales 
puede  construirse  con  una  llamada  a  la  función  diagmatrix, 
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(°/0il2)  diagmatrix(4,°/„e) ; 


(%ol2) 


fe 

0 

0 

0' 

0 

e 

0 

0 

0 

0 

e 

0 

^0 

0 

0 

e, 

En  todo  caso,  debe  tenerse  cuidado  en  que  si  la  matriz  no  se  construye  de  forma 
apropiada,  Maxima  no  la  reconoce  como  tal.  Para  saber  si  una  expresión  es  reconocida 
como  una  matriz  se  utiliza  la  función  matrixp;  la  siguiente  secuencia  permite  aclarar  lo 
que  se  pretende  decir, 


(°/0il3)  matrix(  [[1,2,3]  ,  [4,5,6]]  ) ;  /*  construcción  correcta  */ 


(%ol3)  ([1,2,3]  [4,5,6]) 

(°/0il4)  matrixp(0/0)  ;  /*  es  la  anterior  realmente  una  matriz?  */ 

(%ol4)  true 

C/.Í15)  [[7,8,9]  ,  [0,1,2]]  ;  /*  otra  matriz  */ 


(%ol5)  [[7, 8, 9],  [0,1, 2]] 

(°/0il6)  matrixp (%) ;  /*  será  una  matriz?  */ 


(%ol6)  false 

Casos  particulares  de  submatrices  son  las  filas  y  las  columnas;  los  ejemplos  se  explican 
por  sí  solos: 

(0/„il7)  col  (mi,  3); 

(° 

(%ol7)  -2 

(°/0il8)  row(ml,2) ; 

(%ol8)  (6  0  -2) 

Con  las  matrices  se  pueden  realizar  múltiples  operaciones.  Empezamos  por  el  cálculo 
de  la  potencia  de  una  matriz: 

C/.Í19)  ml~~2; 
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/33  12  -8  \ 

(%ol9)  18  12  -2  a 

\36  6a  a2  —  12/ 

Nótese  que  se  utiliza  dos  veces  el  símbolo  ~  antes  del  exponente;  en  caso  de  escribir¬ 
lo  una  sola  vez  se  calcularían  las  potencias  de  cada  uno  de  los  elementos  de  la  matriz 
independientemente,  como  se  indica  en  el  siguiente  ejemplo, 

(°/0i20)  m2 :ml~2; 


/  9  16  0\ 

( %o38)  36  0  4 

\  0  36  a2) 

Para  la  suma,  resta  y  producto  matriciales  se  utilizan  los  operadores  +,  -  y  . ,  respec¬ 
tivamente, 

(°/0i21)  ml+m2; 


(%o21) 


C/.Í22)  ml-m2; 


'12  20  0  \ 
42  0  2 

v  0  42  a2  +  a) 


( %o22) 

(0/„i23)  ml.m2; 


/  -6  -12  0  \ 

-30  0  -6 

V  0  -30  a  -a2/ 


/ 171  48  16  \ 

(%o23)  54  24  -2  a2 

\216  36  a  a3 +  24/ 


Sin  embargo,  tanto  el  producto  elemento  a  elemento  de  dos  matrices,  como  la  multi¬ 
plicación  por  un  escalar  se  realizan  mediante  el  operador  *,  como  indican  los  siguientes 
dos  ejemplos, 

(°/0i24)  ml*m2; 

/  27  64  0  \ 

( %o24)  216  0  -8 

V  0  216  a3) 


(0/„i25)  4*ml; 
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/12  16  0  \ 

( %o25)  24  0  -8 

\0  24  4  a) 

Otros  cálculos  frecuentes  con  matrices  son  la  transposición,  el  determinante,  la  inver¬ 
sión,  el  polinomio  característico,  así  como  los  valores  y  vectores  propios;  para  todos  ellos 
hay  funciones  en  Maxima: 

(°/0i26)  transpose  (mi) ;  /*la  transpuesta*/ 

/3  6  0\ 

( %o26)  4  0  6 

\0  -2  a) 

(°/0i27)  determinant(ml) ;  /*el  determinante*/ 


(%o27)  36-24o 

(°/0i28)  invert(ml);  /*la  inversa*/ 


(%o28) 


(12  4  g 

m~ña  363a24a 

363h24a  36~ 18° 

36-24  a  36-24  a 


(°/0i29)  invert(ml)  ,detout;  /*la  inversa,  con  el  determinante  fuera*/ 


(%o29) 


/  12  -4o  —8  \ 

J  — 6o  3  a  6  I 
\  36  -18  -24/ 

36-24o 


(°/0i30)  charpoly(ml  ,x) ;  /*pol.  caract.  con  variable  x*/ 


( %o30)  (3  —  x )  (12  —  (a  —  x)  x)  —  24  (a  —  x) 

(°/0i31)  expandí0/) ;  /*pol.  caract.  expandido*/ 

( %o31)  — x 3  +  ax2  +  3x2  —  3ax  +  12a;  —  24a  +  36 

Vamos  a  suponer  ahora  que  a  vale  cero  y  calculemos  los  valores  propios  de  la  matriz, 
(°/0i32)  mi ,  a=0 ; 


(%o32) 


'3  4  O  ' 
6  0-2 
.0  6  0  y 
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(°/0i33)  eigenvalues(°/0)  ; 


( %o33) 


v/l5*  +  3  VÍ5¿-3  " 
- 


,[1,1,1] 


El  resultado  que  se  obtiene  es  una  lista  formada  por  dos  sublistas,  en  la  primera  se 
encuentran  los  valores  propios,  que  en  este  caso  son  Ai  =  —  |  —  A2  =  —  §  +  - i  y 

A3  =  6,  mientras  que  en  la  segunda  sublista  se  nos  indican  las  multiplicidades  de  cada  una 
de  las  A 

Para  el  cálculo  de  los  vectores  propios, 

(°/0i34)  eigenvectors(°/„o32) ; 


( %o34) 


1  _ vd5 1+3  3jy%i/ 5 1  — "21  ]  I",  y/l¿t±§'  3-VbVÁ  i+21 

’  ’W'  ’  rffir  ®h¡  ■# 


Lo  que  se  obtiene  es,  en  primer  lugar,  los  valores  propios  junto  con  sus  multiplicida¬ 
des,  el  mismo  resultado  que  se  obtuvo  con  la  función  eigenvalues,  y  a  continuación  los 
vectores  propios  de  la  matriz  asociados  a  cada  uno  de  los  valores  propios.  A  veces  intere¬ 
sa  que  los  vectores  sean  unitarios,  de  norma  1,  para  lo  que  será  de  utilidad  la  función 
uniteigenvectors,  que  se  encuentra  definida  en  el  paquete  eigen.lisp,  lo  que  significa 
que  antes  de  hacer  uso  de  ella  habrá  que  ejecutar  la  orden  load(eigen) .  También  podemos 
solicitar  los  vectores  propios  unitarios  por  medio  de  la  función  uniteigenvectors, 


(°/0i35)  uniteigenvectors  (°/0o32) ; 


( %o35) 


'[[-J^|.^ta|6]>[ljljl]]> 

"_¿2_  V2VlEi+9V2  3\/2  y/3  \/5¿— 21  \/2 1 
,\/23’  8^23  ’  8V23  J  ’ 

'  V2  V2V15Í-9V2  3V2V3V5Í+21 V2]  [2V2  3V2  3y/2  11 

_\/23  ’  8  \/23  ’  8  \/23  J  ’  L^Í7’  2vT7’  2VT7JJ 


El  rango,  el  menor  de  un  elemento  y  una  base  del  espacio  nulo  de  cierta  matriz  pueden 
calcularse,  respectivamente,  con  las  funciones  rank,  minor  y  nullspace, 


(°/0i36)  rank(°/0o32) ;  /*  el  rango  de  la  matriz*/ 


( %o36)  3 

(°/0i37)  minor(°/0o32,2, 1) ;  /*  el  menor  de  un  elemento  */ 


(%o37) 


A  0' 
6  0 
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(°/0i38)  nullspace(°/0o32) ;  /*  base  del  núcleo  */ 

0  errors,  0  warnings 

( %o38)  span  () 

que  es  la  respuesta  que  se  obtiene  cuando  el  espacio  nulo  está  formado  por  un  único 
elemento. 

De  forma  similar  a  como  la  instrucción  map  aplica  una  función  a  todos  los  elementos 
de  una  lista,  matrixmap  hace  lo  propio  con  los  elementos  de  una  matriz, 

(°/0i39)  matrixmap(sin,70o32)  ; 


( %o39) 


Avancemos  un  poco  más  en  otros  aspectos  del  álgebra  lineal.  Empezamos  con  el  cálculo 
de  la  matriz  de  Jordán  J  de  cierta  matriz  dada  A,  esto  es,  la  matriz  que  verifica  A  =  SJS -1, 
para  cierta  S.  Para  ello  es  necesario  cargar  previamente  el  paquete  diag,  se  introduce  la 
matriz  A  y  la  función  j  ordan  se  encargará  de  obtener  lo  que  buscamos  en  un  formato  que 
luego  utilizamos  para  obtener  las  matrices  J  y  S.  Finalmente,  comprobamos  los  resultados: 


(°/0i40)  A:  matrix  ( [2 , 4 , -6 , 0]  ,  [4, 6,-3, -4]  ,  [0,0, 4,0]  ,  [0,4, -6, 2])  ; 


/2  4  -6 

( %o40) 

4  6-3 

0  0  4 

-4 

0 

\0  4  -6 

2  ) 

(°/0i41)  load  (diag)  $ 

jl:  j  ordan  (A); 

( %o42) 


[[6,1],  [2, 2],  [4,1]] 


(°/0i43)  J:  dispJordan(jl) ; 


( %o43) 


(6 

0 

0 

0' 

0 

2 

1 

0 

0 

0 

2 

0 

^0 

0 

0 

4, 

(°/0i44)  S:  ModeMatrix  (A,jl) 
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( %o44) 

(°/„i45)  is(A  =  S.J.  (S~~-l) )  ; 


(1  4  O  0\ 
10  11 
O  O  O  1 
14  10/ 


( %o45)  true 

Cuando  los  coeficientes  de  un  sistema  lineal  se  tengan  en  forma  matricial  y  los  términos 
independientes  en  un  vector,  el  sistema  se  podrá  resolver  con  la  función  linsolve_by_lu, 

(°/0i46)  A  :  matrix(  [sqrt(2)  ,4]  ,  [a, 4/5] ) ; 


( %o46) 

(°/0i47)  B  :  [1/3, a]; 

( %o47) 

Resolvemos  ahora  AX  =  B  y 

(°/0i48)  linsolve_by_lu(A,B)  ; 


luego  simplificamos  algo  el  resultado 


( %o48) 


(°/0i49)  ratsimp(%)  ; 


l_4(g-372) 
3  £-272a 

y/2 

a~W5 

f-2\/2a 


(  %o49) 


15V2a-V2 
15  \/2a— 6 
(15  a/2-5  )o 
”  60  a— 12\/2 


,  false 


El  resultado  que  se  obtiene  es  una  lista  que  contiene  la  solución  y  el  símbolo  false. 
Cuando  la  resolución  del  sistema  sea  de  índole  numérica,  esta  segunda  componente  será  sus¬ 
tituida  por  el  número  de  condición,  una  medida  de  la  bondad  del  resultado  numérico. 
Veamos  de  aplicar  esto  a  cierta  ecuación  matricial  del  estilo  MX  =  MM 


(°/0i50)  M  :  matrix(  [sqrt(3)  ,2/5]  ,  [3, sin (2)] ) ; 
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( %o50) 


V5  I  ^ 

3  sin  2J 


(°/0i51)  linsolve_by_lu(M,M.M,  ’floatfield)  ; 


(%o51) 


Vi. 732050807568877  ,4000000000000005\ 
_  V3.000000000000001  .9092974268256803 ) 


,49.33731560796254 


Veamos  ahora  cómo  obtener  algunas  matrices  especiales:  la  hessiana  de  una  función, 
la  de  Hilbert,  la  de  Toeplitz,  la  de  Vandermone  y  la  de  Hankel, 


(°/0i52)  hessian(exp(x~2*y+a*z) ,  [x,y,z] ) ; 


(%o52) 


(Ax2  y 2  eaz+x2y  +  2yeaz+x2y 

2x3yeaz+x2y  +  2xeaz+x2y 

2  axyeaz+x2y 


2  x3yeaz+x2y 
xA  ea 
ax2  ec 


2xea' 

x2y 


y 


2  axyeaz+x2y\ 
ax2  eazWy 
a2eaz+x2y  J 


(°/0i53)  hilbert_matrix(5) ; 


( %o53) 


/I 

1 

2 

? 

t  A 

? 

\ 

f 

í  ? 

i 

í 

5 

5  1 

y 

! 

6 

f 

7 

I  iJ 

8  9/ 

(°/0i54)  toeplitz  ( [1 , 2 , 3 , 4]  ,  [t ,  x ,  y ,  z] ) ; 


( %o54) 


(1  x  y  z\ 
2  1  x  y 

3  2  1  x 
4  3  2  1/ 


(°/0i55)  vandermonde_matrix(  [u,v,x,y,z] ) ; 
0  errors,  0  warnings 


/I  u  u 2  u3  uA\ 

1  V  V2  V3  VA 

1  X  X2  X3  XA 

i  y  y2  y3  y 4 

\i  z  V  *3  zA) 


( %o55) 
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(°/0i56)  hankel  ( [v,x,y,z]  ,  [p,q,r,s]  )  ; 


( %o56) 


¡v  x  y  z\ 
x  y  z  q\ 
y  z  q  rj 
\z  q  r  s) 


El  lector  interesado  puede  obtener  información  sobre  otras  funciones  matriciales  acce¬ 
diendo  a  la  documentación  sobre  el  paquete  linearalgebra. 


4.4.  Patrones  y  reglas 

En  un  programa  de  cálculo  simbólico,  éste  no  sólo  debe  tener  información  sobre  una 
función  u  operación  a  partir  de  su  definición,  sino  que  también  habrá  propiedades  y  reglas 
de  transformación  de  expresiones  de  las  que  un  programa  como  Maxima  debe  tener  noticia. 

Nos  planteamos  la  siguiente  situación:  necesitamos  trabajar  con  una  función  G(x,  y) 
que,  independientemente  de  que  esté  definida  o  no,  sabemos  que  es  igual  a  la  expresión 
en  todo  su  dominio,  siendo  H(x,  y)  otra  función;  queremos  que  la  primera  expresión 
sea  sustituida  por  la  segunda  y  además  queremos  tener  bajo  control  estas  sustituciones. 
Todo  ello  se  consigue  trabajando  con  patrones  y  reglas. 

Antes  de  definir  la  regla  de  sustitución  es  necesario  saber  a  qué  patrones  será  aplicable, 
para  lo  cual  admitiremos  que  los  argumentos  de  G(x,  y)  pueden  tener  cualquier  forma: 

(7„il)  matchdeclare  ([x,y],  true) ; 


( %ol)  done 

En  este  caso,  las  variables  patrón  serán  x  e  y,  siendo  el  segundo  argumento  una  función 
de  predicado8  que  devolverá  true  si  el  patrón  se  cumple;  en  el  caso  presente,  el  patrón  es 
universal  y  admite  cualquier  formato  para  estas  variables.  Si  quisiésemos  que  la  regla  se 
aplicase  sólo  a  números  enteros,  se  debería  escribir  matchdeclare  ( [x,y]  ,  integerp) ;  si 
quisiésemos  que  la  regla  se  aplicase  siempre  que  x  sea  un  número,  sin  imponer  restricciones 
a  y,  escribiríamos  matchdeclare  (x,  numberp,  y,  true).  Como  se  ve,  los  argumentos 
impares  son  variables  o  listas  de  variables  y  los  pares  las  condiciones  de  los  patrones. 

Se  define  ahora  la  regla  de  sustitución  indicada  más  arriba,  a  la  que  llamaremos  reglal, 

(°/0i2)  defrule  (reglal,  G(x,y) ,  H(x,y)/x); 


( %o2)  reglai  :  G  (x,  y)  — >  ^  ^ 

Aplicamos  ahora  la  regla  a  las  expresiones  G(f,  2  +  k)  y  G(G(4, 6),  2  +  k), 

8  Se  denominan  así  todas  aquellas  funciones  que  devuelven  como  resultado  de  su  evaluación  true  o  false. 
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(°/0i3)  applyl(G(f ,2+k)  ,reglal) ; 


( %o3) 


H(f,k  +  2) 

/ 


(°/0i4)  applyl(G(G(4,6)  ,2+k)  ,reglal)  ; 


( %o4) 


444(^1^,  k  +  2) 
H(  4,6) 


Como  se  ve  en  este  último  ejemplo,  la  regla  se  aplica  a  todas  las  subexpresiones  que 
contengan  a  la  función  G.  Además,  applyl  aplica  la  regla  desde  fuera  hacia  adentro. 
La  variable  global  maxapplydepth  indica  a  Maxima  hasta  qué  nivel  puede  bajar  en  la 
expresión  para  aplicar  la  regla;  su  valor  por  defecto  es  10000,  pero  se  lo  podemos  cambiar, 


(°/0i5)  maxapplydepth; 


( %o5)  10000 

(°/0i6)  maxapplydepth :  1 ; 


( %o6)  1 

(°/0i7)  applyl(G(G(4,6)  ,2+k)  ,reglal)  ; 


( %o7) 


H  {G  (4, 6) ,  k  +  2) 
G  (4, 6) 


Quizás  sea  nuestro  deseo  realizar  la  sustitución  desde  dentro  hacia  fuera,  controlando 
también  a  qué  niveles  se  aplica;  applybl  y  maxapplyheight  son  las  claves  ahora. 


(°/0i8)  maxapplyheight :  1 ; 


( %o8)  1 

(°/0i9)  applybl (G(G (4, 6)  ,2+k)  ,reglal)  ; 


(%o9) 


Obsérvese  que  hemos  estado  controlando  el  comportamiento  de  las  funciones  G  y  H 
sin  haberlas  definido  explícitamente,  pero  nada  nos  impide  hacerlo, 
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(°/„ilO)  H(u, v)  :  =  u~v+l; 

v 

(°/0olO)  H(u,  v)  :=  u  +  1 

(°/.ill)  applybl  (G (G (4 , 6)  ,2+k)  ,reglal)  ; 


Continuemos  esta  exposición  con  un  ejemplo  algo  más  sofisticado.  Supongamos  que 
cierta  función  F  verifica  la  igualdad 

F(x i  +  X2  +  ■  ■  .  +  xn)  =  F(x i)  +  F(x 2)  +  . . .  +  F(xn)  +  X1X2  ■ . .  xn 

y  que  queremos  definir  una  regla  que  realice  esta  transformación.  Puesto  que  la  regla  se 
aplicará  sólo  cuando  el  argumento  de  F  sea  un  sumando,  necesitamos  una  función  de 
predicado  que  defina  este  patrón, 

(°/0il2)  esunasuma(expr)  :=  not  atom(expr)  and  op(expr)  =  "  +  "  $ 

(°/0il3)  matchdeclare (z , esunasuma) $ 

En  la  definición  de  la  nueva  regla,  le  indicamos  a  Maxima  que  cada  vez  que  se  encuentre 
con  la  función  F  y  que  ésta  tenga  como  argumento  una  suma,  la  transforme  según  se  indica: 
map(F,  args(z))  aplica  la  función  a  cada  uno  de  los  sumandos  de  z,  sumando  después 
estos  resultados  con  la  función  apply,  finalmente  a  este  resultado  se  le  suma  el  producto 
de  los  sumandos  de  z, 

(°/0il4)  def rule  (regla2 , 

F (z) , 

apply ("+" ,map(F,  args(z)))  +  applyC1*",  args(z))); 


( %ol4)  regld2  :  F  (z)  — ►  apply  (+,  map  (F,  args  (z)))  +  apply  (*,  args  (z)) 

Veamos  unos  cuantos  resultados  de  la  aplicación  de  esta  regla, 

(*/,il5)  applyl  (F(a+b) ,  regla2) ; 

(%ol5)  F(b)  +  ab  +  F(a) 

(°/0il6)  applyl  (F(a+b) ,  regla2)  -  applyl  (F (a+c) ,  regla2)  ; 

(%ol6)  ~F(c)  -ac  +  F(b)  +  ab 

(°/0il7)  applyl  (F(a+b+c) ,  regla2) ; 
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(%ol7)  F(c)  +  abc  +  F(b)  +  F(a) 

(0/.il8)  applyl  (F(4)  ,  regla2)  ; 

( %ol8)  F  (4) 

(°/0il9)  applyl  (F (4+5) ,  regla2) ; 

( %ol9)  F  (9) 

(°/0i20)  simp:false$  /*  inhibe  la  simplificación  de  Maxima  */ 

(°/0i21)  applyl  (F  (4+5) ,  regla2) ; 

(%o21)  45  +  (E(4)  +  F(5)) 

(°/0i22)  simp:true$  /*  restaura  la  simplificación  de  Maxima  */ 

C/.Í23)  %o21 ; 

( %o23)  F  (5)  +  F  (4)  +  20 

En  los  ejemplos  recién  vistos,  hemos  tenido  que  indicar  expresamente  a  Maxima  en 
qué  momento  debe  aplicar  una  regla.  Otro  mecanismo  de  definición  de  reglas  es  el  aportado 
por  tellsimp  y  tellsimpafter;  el  primero  define  reglas  a  aplicar  antes  de  que  Maxima 
aplique  las  suyas  propias,  y  el  segundo  para  reglas  que  se  aplican  después  de  las  habituales 
de  Maxima. 

Otra  función  útil  en  estos  contextos  es  declare,  con  la  que  se  pueden  declarar  pro¬ 
piedades  algebraicas  a  nuevos  operadores.  A  continuación  se  declara  una  operación  de 
nombre  o  como  infija  (esto  es,  que  sus  operandos  se  colocan  a  ambos  lados  del  operador); 
obsérvese  que  sólo  cuando  la  operación  se  declara  como  conmutativa  Maxima  considera 
que  a  o  b  es  lo  mismo  que  boa,  además,  no  ha  sido  necesario  definir  qué  es  lo  que  hace 
la  nueva  operación  con  sus  argumentos, 

(°/0i24)  inf  ix  ( "  o " ) ; 

( %o24)  o 

(°/0i25)  is(a  o  b  =  b  o  a); 

( %o25)  false 

(°/0i26)  declare  ("o",  commutative)$ 

(°/0i27)  is(a  o  b  =  b  o  a) ; 


(%o27) 


true 
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5.  Cálculo 

5.1.  Funciones  matemáticas 

En  Maxima  están  definidas  un  gran  número  de  funciones,  algunas  de  las  cuales  se 
presentan  en  la  Figura  6.  A  continuación  se  desarrollan  algunos  ejemplos  sobre  su  uso. 
Empecemos  por  algunos  ejemplos, 

(°/.il)  log(%e) ; 

(%ol)  1 

(°/0i2)  abs(-3~-x) ; 

(%°2)  ^ 

(°/0i3)  signum(-3~-x)  ; 

( %o3)  -1 

La  función  genf  act  (m ,  n ,  p)  es  el  factorial  generalizado,  de  forma  que  geni  act  (m ,  m ,  1 ) 
coincide  con  m!, 

(°/0i4)  genf  act  (5, 5, 1) -5!  ; 

( %o4)  0 

y  genfact(m,m/2,2)  es  igual  a  mü, 

(0/.i5)  genf  act  (5 , 5/2 , 2)  -5 !  !  ; 

( %o5)  0 

Maxima  siempre  devuelve  resultados  exactos,  que  nosotros  podemos  solicitar  en  for¬ 

mato  decimal, 

(°/0i6)  asin(l)  ; 

(%o6)  | 

(0/.i7)  floatC/O; 

( %o7)  1.570796326794897 

Pero  si  damos  el  argumento  en  formato  decimal,  Maxima  también  devuelve  el  resultado 
en  este  mismo  formato,  sin  necesidad  de  hacer  uso  de  f  loat, 
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abs(x) 

abs(x) 

min(xl ,x2,  .  .  . ) 

mín(xi,  X2,  ■  ■  •) 

max(xl ,x2 ,  .  .  . ) 

máx(xi,  X2,  ■  ■  •) 

signum(x) 

r  — i  si  x  <  o 

signo(x)  —  <  0  si  x  =  0 

{  1  six>0 

x ! 

x\ 

x!  ! 

x\\ 

binomial(m,n) 

/mN  _  m(m— (n— 1)J 
\n)  n\ 

genf act (m,n,p) 

m(m  —  p)(m  —  2 p)  . . .  [m  —  (n  —  1 

sqrt (x) 

\fx 

exp(x) 

ex 

log(x) 

ln(x) 

sin(x) 

sin(x) 

cos(x) 

cos(x) 

tan(x) 

tan(x) 

csc(x) 

csc(x) 

sec(x) 

sec(x) 

cot(x) 

cot(x) 

asin(x) 

arcsin(x) 

acos(x) 

are  cos(x) 

atan(x) 

arctan(x) 

atan2(x,y) 

arctan(|)  G  (— 7r,  7 r) 

sinh(x) 

sinh(x)  =  ^(ex  —  e~x ) 

cosh(x) 

cosh(x)  =  7¿(ex  +  e~x) 

tanh(x) 

tanh(x)  = 

asinh(x) 

arcsinh(x) 

acosh(x) 

arccosh(x) 

atanh(x) 

arctanh(x) 

gamma (x) 

T(x)  =  e~uux  ídu,  Vx  >  0 

beta(x,y) 

p(x,y)  =  ',S+íy 

erf (x) 

erf(x)  =  Jq  du 

Figura  6:  Algunas  funciones  de  Maxima. 
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(°/0i8)  asin(l.O); 


( %o8)  1.570796326794897 

Recordemos  que  el  formato  decimal  lo  podemos  pedir  con  precisión  arbitraria, 

(°/0i9)  fpprec:50$  bf  loat(°/0o8) ; 

( %o9)  1.5707963267948966192313216916397514420985846996876b  x  10o 

La  función  de  error  está  relacionada  con  la  función  de  distribución  de  la  variable 

aleatoria  normal  X  ~  A/"(0, 1)  de  la  forma 

$(x)  =  Pr(X  <  x)  =  ^  +  ^erf  (j^j  , 

por  lo  que  la  probabilidad  de  que  esta  variable  tome  un  valor  menor  que  1.5  es 
(°/0ilO)  0 . 5+0 . 5*erf  (1 . 5/ sqrt  (2) )  ,numer; 

(%ol0)  0.9331927987311419 

Una  forma  más  elegante  de  hacer  lo  anterior  es  definir  nuestra  propia  función  de 
distribución  a  partir  de  la  de  error,  para  lo  que  se  hace  uso  del  símbolo  :  =, 

(Xill)  F(x) :=l/2+erf (x/sqrt(2))/2  $ 

(°/0il2)  F(1 . 5)  ,numer; 

(%ol2)  0.9331927987311419 

No  terminan  aquí  las  funciones  de  Maxima;  junto  a  las  ya  expuestas  habría  que  incluir 
las  funciones  de  Airy,  elípticas  y  de  Bessel,  sobre  las  que  se  podrá  obtener  más  información 
ejecutando  la  instrucción  ??  y  utilizando  como  argumento  airy,  elliptic  o  bessel,  según 
el  caso.  Por  ejemplo, 

C/.Í13)  ??  airy 

0:  airy_ai  (Funciones  y  variables  para  las  funciones  especiales) 

1:  airy_bi  (Funciones  y  variables  para  las  funciones  especiales) 

2:  airy_dai  (Funciones  y  variables  para  las  funciones  especiales) 

3:  airy_dbi  (Funciones  y  variables  para  las  funciones  especiales) 

Enter  space-separated  numbers,  ‘all’  or  ‘none’ :  0 

—  Función:  airy_ai  (<x>) 

Función  Ai  de  Airy,  tal  como  la  definen  Abramowitz  y  Stegun, 

Handbook  of  Mathematical  Functions,  Sección  10.4. 
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La  ecuación  de  Airy  ‘diff  (y(x) ,  x,  2)  -  x  y(x)  =  O’  tiene  dos 
soluciones  linealmente  independientes,  ‘y  =  Ai(x)’  y  ‘y  =  Bi(x)’. 

La  derivada  ‘diff  (airy_ai(x),  x) 5  es  ‘airy_dai(x) ’ . 

Si  el  argumento  ‘x’  es  un  número  decimal  real  o  complejo,  se 
devolverá  el  valor  numérico  de  ‘airy_ai’  siempre  que  sea  posible. 

Véanse  ‘airy_bi’,  ‘airy_dai’  y  ‘airy_dbi’ . 

Maxima  reconoce  los  dominios  en  el  plano  complejo  de  las  funciones;  los  siguientes 
ejemplos  lo  demuestran: 

(0/,il4)  sin(%i) ; 

( %ol4)  i  sinh  1 

C/.Í15)  log(-3 . 0) ; 

(%ol5)  3.141592653589793  ¿  +  1.09861228866811 

(°/0il6)  asin(2 . 0) ; 

(%ol6)  1.570796326794897  -  1.316957896924817  i 

5.2.  Límites 

Sin  más  preámbulos,  veamos  algunos  ejemplos  de  cómo  calcular  límites  con  la  asistencia 
de  Maxima.  En  primer  lugar  vemos  que  es  posible  hacer  que  la  variable  se  aproxime  al 
infinito  ( x  — ►  oo)  haciendo  uso  del  símbolo  inf,  o  que  se  aproxime  al  menos  infinito 
(x  — ►  — oo)  haciendo  uso  de  minf , 

(•/.il)  limit  (1/ sqrt  (x)  ,x, inf ) ; 

(%ol)  0 

(°/0i2)  limit ((exp(x)-exp(-x))/(exp(x)+exp(-x))  ,x,minf) ; 


( %o2)  -1 

que  nos  permite  calcular  lím^oo  y  lím^^-oo  I ,  respectivamente. 

Los  siguientes  ejemplos  muestran  límites  en  los  que  la  variable  x  se  aproxima  a  puntos 
de  discontinuidad, 


(°/„i3)  limit ( (x~2-x/3-2/3) / (5*x~2-ll*x+6)  ,x,l) ; 
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( %o3) 

(°/0i4)  limit(l/(x-l)~2,x,l) ; 


(%o4) 


oo 


donde  hemos  obtenido  los  resultados 


x  —  £  —  ñ 


5 

3 


y 


1  ( x  —  l)2 


=  oo. 


Sin  embargo,  ciertos  límites  no  se  pueden  resolver  sin  aportar  información  adicional, 
tal  es  el  caso  de  \ímx^\  para  el  que  podemos  hacer 


(°/0i5)  limit  (l/(x-l)  ,x,  1) ; 


( %o5)  und 

donde  Maxima  nos  responde  con  el  símbolo  und  de  undeftned,  o  indefinido.  En  tales  si¬ 
tuaciones  podemos  indicarle  al  asistente  que  la  variable  x  se  aproxima  a  1  por  la  derecha 
(x  —>  1+)  o  por  la  izquierda  (x  — >  1_), 

(°/0i6)  limit  (1/ (x-1)  ,x,  1  .plus) ; 


( %o6)  oo 

(°/0i7)  limit  (1/ (x-1)  ,x,  1  ,minus) ; 


( %o7)  — oo 

5.3.  Derivadas 

Maxima  controla  el  cálculo  de  derivadas  mediante  la  instrucción  diff .  A  continuación 
se  presentan  algunos  ejemplos  sobre  su  uso, 

(°/0il)  diff  (x~log(a*x)  ,x) ;  /*  primera  derivada  */ 


(%ol) 


log^)  ( logM 


logx\ 

X  ) 


(°/0i2)  diff  (x~log(a*x)  ,x,2)  ;  /^derivada  segunda*/ 
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( %o2)  *««*>  (%M  +  M!  +  jí„gM  /íog^iog*  2\ 

\  X  X  )  \  x¿  x¿  x¿ ) 

(°/0i3)  factor (%) ;  /*  ayudamos  a  Maxima  a  mejorar  el  resultado  */ 

( %o3)  xl°s{ax)-2  ^Qg2  _|_  2  log  x  log  (a  x)  —  log  (a x)  +  log2  x  —  log  x  +  2) 

donde  se  han  calculado  la  primera  y  segunda  derivadas  de  la  función  y  =  x;ln(axh  Nótese 
que  en  la  entrada  (0/„i3)  se  le  ha  pedido  al  programa  que  factorizase  la  salida  (0/„o2) . 

Pedimos  ahora  a  Maxima  que  nos  calcule  el  siguiente  resultado  que  implica  derivadas 
parciales, 

$io 

VxHyHz2  (e*sinÜ/)  tan(¿0)  =  2e*cos(y)sec2(z)  tau(z). 

(°/0i4)  diff (exp(x)*sin(y)*tan(z)  ,x,3,y,5,z,2)  ; 


( %o4)  2  ex  eos  y  sec2  z  tan  z 

Maxima  también  nos  puede  ayudar  a  la  hora  de  aplicar  la  regla  de  la  cadena  en  el 
cálculo  de  derivadas  de  funciones  vectoriales  con  variable  también  vectorial.  Supónganse 
que  cierta  variable  z  depende  de  otras  dos  x  y  y,  las  cuales  a  su  vez  dependen  de  u  y  v. 
Veamos  cómo  se  aplica  la  regla  de  la  cadena  para  obtener  |^,  o 
(°/0i5)  depends  (z ,  [x ,  y]  ,  [x ,  y]  ,  [u ,  v]  )  ; 


(%o5)  [z(x,y)  ,x(u,v)  ,y(u,v)] 

(°/0i6)  diff  (z,v,  1)  ; 


( %o6) 

(°/0i7)  diff  (z ,y,  1  ,v,  1)  ; 


'd_  \ 
dy  /  +  d 


(  — 

\dx 


( %o7) 

(°/0i8)  diff  (z ,u,  1  ,v,  1)  ; 


A  ( (j2 

dvV  \dy2 


A 

d 


(— 

\dxdy 


( %o8) 


+  dh¡L>y  (áz) + 

+  dVdü  x  (A  z ) 


La  función  depends  admite  un  número  par  de  argumentos,  los  cuales  pueden  ser  varia¬ 
bles  o  listas  de  variables.  Las  que  ocupan  los  lugares  impares  son  las  variables  dependientes 
y  las  que  están  en  las  posiciones  pares  las  independientes;  denominación  ésta  que  tiene 
aquí  un  sentido  relativo,  ya  que  una  misma  variable,  como  la  x  de  la  entrada  %i5,  puede 
actuar  como  independiente  respecto  de  z  y  como  dependiente  respecto  de  u  y  v. 

En  cualquier  momento  podemos  solicitarle  a  Maxima  que  nos  recuerde  el  cuadro  de 
dependencias, 
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(°/0i9)  dependencies ; 

( %o9)  [z  ( x ,  y)  ,x(u,v)  ,y  ( u ,  -y)] 

También  podemos  eliminar  algunas  de  la  dependencias  previamente  especificadas, 
(°/0ilO)  remove (x,dependency) ; 

( %ol0)  done 

(lili)  dependencies; 

(%oll)  [z(x,y)  ,y(u,v)] 

(°/0il2)  diff  (z ,y,  1  ,v,  1) ; 


( %o!2) 


Veamos  cómo  deriva  Maxima  funciones  definidas  implícitamente.  En  el  siguiente  ejem¬ 
plo,  para  evitar  que  y  sea  considerada  una  constante,  le  declararemos  una  dependencia 
respecto  de  x, 


(°/0il3)  depends(y,x)$ 


(°/„il4)  diff  (x~2+y~3=2*x*y,x) ; 


( %oi4)  3»2  (Jxy)+2x=ÍX{Áy)+2y 

Cuando  se  solicita  el  cálculo  de  una  derivada  sin  especificar  la  variable  respecto  de  la 
cual  se  deriva,  Maxima  utilizará  el  símbolo  del  para  representar  las  diferenciales, 

(°/„i!5)  diff  (x~2) ; 


(%ol5)  2  xdel(x) 

lo  que  se  interpretará  como  2 xdx.  Si  en  la  expresión  a  derivar  hay  más  de  una  variable, 
habrá  diferenciales  para  todas, 


(°/„i!6)  diff  (x~2+y~3=2*x*y)  ; 
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(%ol6) 

3  y 2  del  (y)  +  ^3  y2  +  2  del  (x)  =  2  x  del  (y)  +  ^2  x  y^j  +  2yj  del  (x) 

Recuérdese  que  durante  este  cálculo  estaba  todavía  activa  la  dependencia  declarada 
en  la  entrada  (°/0il3). 

Finalmente,  para  acabar  esta  sección,  hagamos  referencia  al  desarrollo  de  Taylor  de 
tercer  grado  de  la  función 

xlnx 


en  el  entorno  de  x  =  1, 

(°/0il7)  taylor( (x*log(x) )/(x~2-l)  ,x,l,3) ; 

1  (x  —  l)2  (x  —  1) 


( %ol7) 

(°/0il8)  expand(°/0) ; 


12 


12 


3 

-  +  ■■■ 


A  continuación  un  ejemplo  de  desarrollo  multivariante  de  la  función  y  =  exp(x2  sin(x'y)) 
alrededor  del  punto  (2,0)  hasta  grado  2  respecto  de  cada  variable, 

(°/0il9)  taylor  (exp(x~2*sin(x*y) )  ,  [x,2,2]  ,  [y, 0,2]  )  ; 

(%ol9) 

1  +  8í/  +  32í/2  +  --  -  +  (12  j/  +  96  y2  +  •  •  • )  (x  -  2)  +  (6y  +  120  í/2  +  •  •  • )  (x  -  2)2  +  •  •  • 
(°/0i20)  expand(°/0) ; 

( %o20)  120  x2  y2  —  384  x  y2  +  320  y2  +  6 x2  y  —  12  x  y  +  8y  +  1 

En  ocasiones  resulta  necesario  operar  con  funciones  cuyas  derivadas  son  desconocidas, 
quizás  porque  la  propia  función  lo  es.  Esta  situación  puede  llevar  a  que  Maxima  devuelva 
expresiones  realmente  complicadas  de  manipular.  Un  ejemplo  es  el  siguiente,  en  el  que 
trabajamos  con  una  función  /  arbitraria 

(°/0i21)  taylor (f  (x  +  x~2)  ,x,  1 , 1) ; 


(%o2l)  /(2)+(^;/(x2  +  x;)|  J  (x-  1)H - 

Podemos  facilitar  las  cosas  si  le  definimos  una  función  derivada  a  /,  a  la  que  llamaremos 

df, 

(0/.i22)  gradef  (f  (x)  ,  df  (x) )  $ 

(°/0i23)  taylor (f  (x+x~2)  ,x,  1 , 1) ; 


(%o23)  /  (2)  +  3  df  (2)  (x  —  1)  +  •  •  • 

El  paquete  pdif  f ,  que  se  encuentra  en  la  carpeta  share/contrib/,  aporta  una  solución 
alternativa  a  este  problema. 
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5.4.  Integrales 

La  función  de  Maxima  que  controla  el  cálculo  de  integrales  es  intégrate,  tanto  para 
las  definidas  como  indefinidas;  empecemos  por  estas  últimas, 

C/.il)  integrate(cos(x)~3/sin(x)~4,x) ; 


(%ol) 

(°/.i2) 


3  sin2  x  —  1 
3  sin3  x 

intégrate (a [3] *x~3+a[2] *x~2+a[l] *x+a[0] ,x) ; 


( %o2) 


a 3  x  02  x°  a  i  x 

4  +  3  +  2 

que  nos  devuelve,  respectivamente,  las  integrales 


/eos3  x 
.  4  dx 
sm  x 

J  (ü3X3  +  Ü2X2  +  CL\X  +  CLo)dx. 


Además,  este  último  ejemplo  nos  ofrece  la  oportunidad  de 
con  subíndices. 

Ahora  un  par  de  ejemplos  sobre  la  integral  definida, 


cómo  escribir  coeficientes 


(°/0i3)  intégrate (2*x/ ((x-l)*(x+2))  ,x,3,5)  ; 


( %o3) 

(°/0i4)  yo,numer; 


( 2  log7  2  log 5  log4  log2\ 
As  3  +  3  3  ) 

/^aproximación  decimal*/ 


( %o4)  0.9107277692015807 

(°/0i5)  intégrate (asin(x)  ,x,0,u) ; 

Is  u  positive,  negative,  or  zero? 

positive; 


(%o5) 


u  aresin  u  +  \ 


2x 

(x  -  \){x  +  2) 


dx  «  0.91072776920158 
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/  arcsin(x)o?x  =  u  arcsin(u)  +  \/l  —  u2  —  1,  Vu  >  0. 

Jo 

Nótese  en  este  último  ejemplo  cómo  antes  de  dar  el  resultado  Maxima  pregunta  si  u  es 
positivo,  negativo  o  nulo;  tras  contestarle  escribiendo  positive;  (punto  y  coma  incluido) 
obtenemos  finalmente  el  resultado.  En  previsión  de  situaciones  como  esta,  podemos  darle 
al  sistema  toda  la  información  relevante  sobre  los  parámetros  utilizados  antes  de  pedirle 
el  cálculo  de  la  integral, 

(°/0i6)  assume(u>0) ; 

( %o6)  [u  >  0] 

(°/0i7)  intégrate  (asín (x)  ,x,0,u) ; 

( %o7)  u  arcsin  u  +  \fl  —  u2  —  1 

Cuando  Maxima  no  puede  resolver  la  integral,  siempre  queda  el  recurso  de  los  métodos 
numéricos.  El  paquete  quadpack,  escrito  inicialmente  en  Fortran  y  portado  a  Lisp  para 
Maxima,  es  el  encargado  de  estos  menesteres;  dispone  de  varias  funciones,  pero  nos  deten¬ 
dremos  tan  sólo  en  dos  de  ellas,  siendo  la  primera  la  utilizada  para  integrales  definidas  en 
intervalos  acotados, 

(°/0i8)  /*  El  integrador  simbólico  no  puede  con  esta  integral  */ 
intégrate (exp (sin (x)) ,x,2,7) ; 


(°/0i9)  /*  Resolvemos  numéricamente  */ 
quad_qag(exp(sin(x)) ,x,2,7,3) ; 


( %o9)  [4.747336298073747, 5.27060206376023  x  1(T14, 31, 0] 

La  función  quad_qag  tiene  un  argumento  extra,  que  debe  ser  un  número  entero  entre 
1  y  6,  el  cual  hace  referencia  al  algoritmo  que  la  función  debe  utilizar  para  la  cuadratura; 
la  regla  heurística  a  seguir  por  el  usuario  es  dar  un  número  tanto  más  alto  cuanto  más 
oscile  la  función  en  el  intervalo  de  integración.  El  resultado  que  obtenemos  es  una  lista  con 
cuatro  elementos:  el  valor  aproximado  de  la  integral,  la  estimación  del  error,  el  número  de 
veces  que  se  tuvo  que  evaluar  el  integrando  y,  finalmente,  un  código  de  error  que  será  cero 
si  no  surgieron  problemas. 

La  otra  función  de  integración  numérica  a  la  que  hacemos  referencia  es  quad_qagi,  a 
utilizar  en  intervalos  no  acotados.  En  el  siguiente  ejemplo  se  pretende  calcular  la  proba¬ 
bilidad  de  que  una  variable  aleatoria  x2  de  n  =  4  grados  de  libertad,  sea  mayor  que  la 
unidad  (Pr(x4  >  1)), 
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(°/.ilO)  n:4$ 

(Xill)  intégrate (x~ (n/2-l)*exp(-y/2)/2~ (n/2) *gamma(n/2) ,x, 1 ,  inf ) ; 
Integral  is  divergent 

—  an  error.  To  debug  this  try  debugmode (true) ; 

(°/0il2)  quad_qagi(x~  (n/2-l)*exp(-x/2)/2~  (n/2) *gamma(n/2)  ,x,  1  ,  inf ) ; 


( %ol2)  [.9097959895689501, 1.913452127046495  x  10"10, 165, 0] 

El  integrador  simbólico  falla  emitiendo  un  mensaje  sobre  la  divergencia  de  la  integral. 

La  función  quad_qagi  ha  necesitado  165  evaluaciones  del  integrando  para  alcanzar  una 
estimación  numérica  de  la  integral,  la  cual  se  corresponde  aceptablemente  con  la  estimación 
que  hacen  los  algoritmos  del  paquete  de  distribuciones  de  probabilidad  (ver  Sección  6.1). 

Otras  funciones  de  cuadratura  numérica  son  quad_qags,  quad_qawc,  quad_qawf ,  quad_qawo 
y  quad_qaws,  cuyas  peculiaridades  podrá  consultar  el  lector  interesado  en  el  manual  de 
referencia. 

La  transformada  de  Laplace  de  una  función  f(x)  se  define  como  la  integral 
L(p)  =  [  f(x)e~pxdx, 

Jo 

siendo  p  un  número  complejo.  Así,  la  transformada  de  Laplace  de  f(x)  =  ke~kx  es 
(°/0il3)  laplace (k*exp(-k*x)  ,x,p) ; 


(%ol3) 


k 

p  +  k 


y  calculando  la  transformada  inversa  volvemos  al  punto  de  partida, 


(0/.il4)  ilt C/0,p,x) ; 


(%ol4)  ke~kx 

La  transformada  de  Fourier  de  una  función  se  reduce  a  la  de  Laplace  cuando  el  argu¬ 
mento  p  toma  el  valor  —  it,  siendo  i  la  unidad  imaginaria  y  t  €  M, 

F(t)  =  í  f(x)eltxdx. 

Jo 

De  esta  manera,  la  transformada  de  Fourier  de  f(x)  =  ke~kx  es 
(°/0il5)  laplace (k*exp(-k*x)  ,x,-°/0i*t) ; 


(%ol5) 


k 

k  —  it 


Nótese  que  si  x  >  0,  la  f(x)  anterior  es  precisamente  la  función  de  densidad  de  una 
variable  aleatoria  exponencial  de  parámetro  k,  por  lo  que  este  último  resultado  coincide 
precisamente  con  la  función  característica  de  esta  misma  distribución. 
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5.5.  Ecuaciones  diferenciales 

Con  Maxima  se  pueden  resolver  analíticamente  algunas  ecuaciones  diferenciales  ordi¬ 
narias  de  primer  y  segundo  orden  mediante  la  instrucción  ode2. 

Una  ecuación  diferencial  de  primer  orden  tiene  la  forma  general  F(x,  y,  y')  =  0,  donde 
y'  =  %■  Para  expresar  una  de  estas  ecuaciones  se  hace  uso  de  diff , 

C/.il)  /*  ecuación  de  variables  separadas  */ 
ec : (x-1) *y~3+(y-l) *x~3* ’diff (y ,x)=0; 


(%ol)  x3  {y-1)  (j^y'j +{x-l)  y3  =  0 

siendo  obligatorio  el  uso  de  la  comilla  simple  ( ’ )  antes  de  diff  al  objeto  de  evitar  el  cálculo 
de  la  derivada,  que  por  otro  lado  daría  cero  al  no  haberse  declarado  la  variable  y  como 
dependiente  de  x.  Para  la  resolución  de  esta  ecuación  tan  solo  habrá  que  hacer 

(°/0i2)  ode2  (ec ,  y ,  x)  ; 


( %o2) 


2y-l 
2  y2 


2x  —  1 
2  x2 


donde  7.C  representa  una  constante,  que  se  ajustará  de  acuerdo  a  la  condición  inicial  que 
se  le  imponga  a  la  ecuación.  Supóngase  que  se  sabe  que  cuando  x  =  2,  debe  verificarse  que 
y  =  —3,  lo  cual  haremos  saber  a  Maxima  a  través  de  la  función  icl, 


(°/0i3)  icl(%o2,x=2,y=-3) ; 


(%o3) 


2y-l  _  x2  +  72x-36 
2y-  ~  72 

Veamos  ejemplos  de  otros  tipos  de  ecuaciones  diferenciales  que  puede  resolver  Maxima, 


(°/0i4)  /*  ecuación  homogénea  */ 

ode2(x~3+y~3+3*x*y~2*’diff (y,x) ,y,x) ; 


( %o4) 


4  xy6 


En  este  caso,  cuando  no  se  incluye  el  símbolo  de  igualdad,  se  da  por  hecho  que  la  expresión 
es  igual  a  cero. 


(°/0i5)  /*  reducible  a  homogénea  */ 

ode2(’diff (y,x)=(x+y-l)/ (x-y-1) ,y,x) ; 


log  (r/2  +  x2  —  2  x  +  l)  +  2  arctan 


(%o5) 


4 


=  %c 
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(°/0i6)  /*  ecuación  exacta  */ 

ode2 ( (4*x~3+8*y) + (8*x-4*y~3) * ’ dif f (y , x) ,y,x) ; 


— y 4  +  8  x  y  +  x4  = 


(°/0i7)  /*  Bernoulli  */ 

ode2 ( ’ dif f (y , x) -y+sqrt (y) , y , x) ; 


( %o7)  2  log  (y/y  —  1)  =  x  +  %c 

(0/.i8)  solve(°/„,y); 


( %o8)  | y  =  ex+  %c  +  2  ei+^  +  1  j 

En  este  último  caso,  optamos  por  obtener  la  solución  en  su  forma  explícita. 

Una  ecuación  diferencial  ordinaria  de  segundo  orden  tiene  la  forma  general  F(x,  y,  y',  y ") 
0,  siendo  y"  la  segunda  derivada  de  y  respecto  de  x.  Como  ejemplo, 

(°/„i9)  >  dif  f  (y ,  x)  =x+  ’  dif  f  (y ,  x ,  2) ; 


( %o9) 

(°/0ilO)  ode2(7OJy,x)  ; 


d 

TxV 


o 

dx 2 


y  +  x 


( %ol0)  y  =  %kt  ex  +  +  22X  +  2  +  %k2 

Maxima  nos  devuelve  un  resultado  que  depende  de  dos  parámetros,  °/0kl  y  °/0k2,  que  para 
ajustarlos  necesitaremos  proporcionar  ciertas  condiciones  iniciales;  si  sabemos  que  cuando 
x  =  1  entonces  y  =  —  1  y  2//  =  |  =2,  haremos  uso  de  la  instrucción  ic2, 

(lili)  ic2(°/0,x=l  ,y=-l  ,diff  (y,x)=2)  ; 


(%oll) 


y  = 


x2  +  2x  +  2 
2 


7 

2 


En  el  caso  de  las  ecuaciones  de  segundo  orden,  también  es  posible  ajustar  los  paráme¬ 
tros  de  la  solución  especificando  condiciones  de  contorno,  esto  es,  fijando  dos  puntos  del 
plano  por  los  que  pase  la  solución;  así,  si  la  solución  obtenida  en  C/„olO)  debe  pasar  por 
los  puntos  (—1,3)  y  (2,  |),  hacemos 

(°/„i!2)  bc2  C/0olO ,  x=-l ,  y=3 ,  x=2 ,  y=5/3)  ; 
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(%ol2) 


35  ex+1  x2  +  2x  +  2  15e3  +  20 

6e3-6  +  2  +  6  e3  —  6 


Nótese  que  este  cálculo  se  le  solicita  a  Maxima  con  bc2. 

La  resolución  de  sistemas  de  ecuaciones  diferenciales  se  hace  con  llamadas  a  la  función 
desolve.  En  este  contexto  es  preciso  tener  en  cuenta  que  se  debe  utilizar  notación  funcional 
dentro  de  la  expresión  diff ;  un  ejemplo  aclarará  este  punto,  resolviendo  el  sistema 


í  W  •  3f(x)  —  2g(x) 
1  W  =  2/(x)-29(z) 
(°/0il3)  desolve(  [’diff  (f  (x)  ,x)=3*f  (x)-2*g(x)  , 
’diff (g(x) ,x)=2*f (x)-2*g(x)] , 

[f (x) ,g(x)]) ; 


(%o!3) 


(2g(0)-/Cg)> f~f.  (2~¿(Ó)--4¿tO)) 

g(x)  /?»))«■*  L  ím-ifmt 


Como  se  ve,  las  referecias  a  las  funciones  deben  incluir  la  variable  independiente  y 
las  ecuaciones  estarán  acotadas  entre  corchetes,  así  como  los  nombres  de  las  funciones. 
Observamos  en  la  respuesta  que  nos  da  Maxima  la  presencia  de  f  (0)  y  g(0),  lo  cual  es 
debido  a  que  se  desconocen  las  condiciones  de  contorno  del  sistema. 

En  este  último  ejemplo,  supongamos  que  queremos  resolver  el  sistema  de  ecuaciones 
diferenciales 

í  NfaT  =  f(x)  +  g(x)  +  3h(x) 

*£?  =  9(x)~2h(x) 

{  1?  =  f(*)  +  h(x ) 

bajo  las  condiciones  /(O)  =  —1,  g( 0)  =  3  y  /(O)  =  1.  En  primer  lugar  introduciremos 
estas  condiciones  con  la  función  at valué,  para  posteriormente  solicitar  la  resolución  del 
sistema, 


(°/0il4)  atvalue(f  (x)  ,x=0,-l)$ 

(°/0il5)  atvalue(g(x)  ,x=0,3)$ 

(°/0il6)  atvalue(h(x)  ,x=0,l)$ 

(°/0il7)  desolve ( [’diff  (f  (x)  ,x)=f  (x)+g(x)+3*h(x) , 

’diff (g (x) ,x)=g(x)-2*h(x) , 

’diff (h(x) ,x)=f (x)+h(x)] ,  [f (x) ,g(x) ,h(x)] ) ; 

( %ol7) 

[/  (x)  =  xe2x  +  e2x  —  2  e~x,g(x)  =  —2xe2x  +  2e2x  +  e~x,  h  (x)  =  xe2x  +  e~x ] 

La  función  desolve  también  nos  va  a  permitir  resolver  ecuaciones  de  orden  mayor  que 
dos.  En  el  siguiente  ejemplo,  abordamos  la  resolución  de  la  ecuación 

d3y(x)  d2y(x)  dy(x)  _ 

dx3  dx2  dx 

bajo  las  condiciones  que  se  describen  a  continuación, 
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(°/0il8)  atvalue ( ’diff  (y(x)  ,  x,2)  ,x=0,v)$ 

(°/0il9)  atvalue ( ’ dif f  (y (x)  ,x)  ,x=0,u)$ 

(°/0i20)  atvalue  (y (x)  ,  x=0 ,  w)  $ 

Ahora  resolvemos, 

(°/0i21)  desolve( ’diff  (y(x)  , x,3)+’diff  (y(x)  ,x,2)+’diff  (y(x)  ,x)+y(x)=0,  y(x) ) ; 


(w  +  v  +  2  u)  sin  x  (w  —  v)  ( 
:  2  +  2 


El  paquete  plotdf  permite 
dales  de  primer  orden 


bien  de  sistemas 


Los  argumentos  a  pasar  a  la  función  plotdf  son  la  función  F,  en  el  primer  caso,  y 
una  lista  con  las  funciones  F  y  G  en  el  segundo.  Las  variables  serán  siempre  x  e  y.  Como 


ejemplo,  pidamos  a  Maxima  que  genere  el  campo  de  direcciones  de  la  ecuación  diferencial 

t  =  l  +  y+y 2 


(°/0i22)  load  (plotdf  )$ 

(°/0i23)  plotdf  (1  +  y  +  y~2)  ; 


El  gráfico  que  resulta  es  el  de  la  Figura  7  a),  en  el  que  además  se  observan  dos  trayec¬ 
torias  que  se  dibujaron  de  forma  interactiva  al  hacer  clic  sobre  dos  puntos  del  plano. 

La  función  plotdf  admite  varias  opciones,  algunas  de  las  cuales  aprovechamos  en  el 
siguiente  ejemplo.  Supongamos  el  modelo  predador-presa  de  Lotka-Volterra,  dependiente 


de  dos  parámetros  h  y  k, 


í  %  =  2x  +  hxy 
l  ^  =  -x  +  kxy 


El  siguiente  código  permite  generar  el  campo  de  direcciones  correspondiente,  dándoles 
inicialmente  a  h  y  k  los  valores  -1.2  y  0.9,  respectivamente;  además,  aparecerán  sobre  la 
ventana  gráfica  dos  barras  de  deslizamiento  para  alterar  de  forma  interactiva  estos  dos 
parámetros  y  ver  los  cambios  que  se  producen  en  el  campo. 


(°/0i24)  plotdf  ( [2*x+k*x*y,  -y+h*x*y]  , 

[parameters, "k=-l .2,h=0.9"]  , 
[sliders, "k=-2:2,h=-2:2"] ) ; 
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a)  b) 


Figura  7:  Campos  de  direcciones  creados  con  la  función  ’plotdf’:  a)  campo  de  la  ecuación 
^|  =  1  +  y  +  y2',  b)  campo  correspondiente  al  modelo  predador- presa. 


En  la  Figura  7  b)  se  observa  el  gráfico  obtenido  despúes  de  pedir  trayectorias  concretas 
que  pasan  por  varios  puntos  del  plano  (en  el  archivo  gráfico  no  aparecen  las  barras  de 
deslizamiento) . 

Cuando  Maxima  no  es  capaz  de  resolver  la  ecuación  propuesta,  se  podrá  recurrir  al 
método  numérico  de  Runge-Kutta,  el  cual  se  encuentra  programado  en  el  paquete  dif  f  eq. 
Como  primer  ejemplo,  nos  planteamos  la  resolución  de  la  ecuación 

^  =  -2y2  +  exp(— 3í), 

con  la  condición  y( 0)  =  1.  La  función  ode2  es  incapaz  de  resolverla: 

(°/0i25)  ec:  ’diff (y,t)+2*y~2-exp(-3*t)=0; 

(%o25)  Aty  +  2y2-e- 3Í  =  0 

(°/0i26)  ode2(ec,y,t) ; 

( %o26)  false 

Abordamos  ahora  el  problema  con  un  enfoque  numérico,  para  lo  cual  definimos  la 
expresión 

f(t,  y)  =  áJt  =  -2 y1  +  exp(— 3í) 


en  Maxima, 
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(0/„i27)  load(diffeq)$ 

(0/.i28)  f  (t,y)  :=  -2*y~2+exp(-3*t)  $ 
(°/0i29)  res:  rungel (f ,0,5 ,0 . 5, 1) ; 


( %o29) 

[[0.0, 0.5, 1.0, 1.5,  2.0, 2.5, 3.0,  3.5, 4.0, 4.5, 5.0]  , 

[1.0,  .5988014211752297,  .4011473182183033,  .2915932807721147  , 

.2260784415237641, 0.183860316087325,  .1547912058210609,  .1336603954558797, 
.1176289334565761,  .1050518038293819,  .09491944625388439] , 

[-1.0,  -0.49399612385452,  -.2720512734596094,  -.1589442862446483  , 
-.09974417126696171,  -.06705614729331426,  -.04779722499498942,  -.03570266617749457, 
-.02766698775990988,  -.02207039201652747,  -.01801909665196759]] 


(°/0i30)  draw2d( 

points_joined  =  true, 
point_type  =  dot, 
points (res [1] ,res [2] ) , 
terminal  =  eps)$ 

La  función  rungel  necesita  cinco  argumentos:  la  función  derivada  ^|,  los  valores  inicial, 
¿o,  y  final,  ¿i,  de  la  variable  independiente,  la  amplitud  de  los  subintervalos  y  el  valor  que 
toma  y  en  to-  El  resultado  es  una  lista  que  a  su  vez  contiene  tres  listas:  las  abscisas  t. 
las  ordenadas  y  y  las  correspondientes  derivadas.  Al  final  del  ejemplo  anterior,  se  solicita 
la  representación  gráfica  de  la  solución,  cuyo  aspecto  es  el  mostrado  por  la  Figura  8  a). 
(Consúltese  la  sección  7  para  una  descripción  de  la  función  draw2d.) 

Nos  planteamos  ahora  la  resolución  de  la  ecuación  diferencial  de  segundo  orden 


dt 2 


-  o.2(i 


con  2/(0)  =  0.1  y  y'(0)  =  1  para  lo  cual  definimos  en  Maxima  la  función 


g(t,y,y') 


fy 

di2 


=  0.2(1  -  y2)y'  -  y, 


(0/.i31)  g(t ,y,yp)  :=  0.2*(l-y~2)*yp  -  y  $ 

(°/0i32)  res:  runge2(g,0, 100 ,0 . 1 ,0 . 1 , 1)$ 

(°/0i33)  midibujo(i,  j)  :=  draw2d(points_joined  =  true, 

point_type  =  dot, 

points (res [i] ,res  [j] ) , 
terminal  =  eps)$ 

(°/0i34)  midibujo  (1 , 2)  $ 

(°/0i35)  midibujo (2,3)$ 

(°/0i36)  midibujo (2,4)$ 


La  función  runge2  necesita  seis  argumentos:  la  función  derivada  los  valores  inicial, 
ío,  y  final,  ¿i,  de  la  variable  independiente,  la  amplitud  de  los  subintervalos  y  los  valores 
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c)  d) 


Figura  8:  Resolución  numérica  de  ecuaciones  diferenciales  con  ’diffeq’:  a),  solución  de  la 
ecuación  de  primer  orden  ^  =  —2 y2  +  exp(— 3£),  y(0)  =  1;  b),  solución  de  la  ecuación  de 
segundo  orden  í¡^  =  0.2(1  —  y2)  ^  —  y ,  con  las  condiciones  y  {ti)  =  0.1  y  y' (ti)  =  1;  c), 
diagrama  (y, y');  d),  diagrama  (y,  y"). 


que  toman  y  y  su  primera  derivada  en  ¿o-  El  resultado  es  una  lista  que  a  su  vez  contiene 
cuatro  listas:  las  abscisas  t,  las  ordenadas  y,  las  correspondientes  primeras  derivadas  y, 
por  último,  las  segundas  derivadas.  Al  final  de  este  último  ejemplo  se  solicitan  algunos 
gráficos  asociados  a  la  solución,  los  formados  con  los  pares  ( t,y ),  (y,  y')  y  (y,  y"),  que  son 
los  correspondientes  a  los  apartados  b),  c)  y  d)  de  la  Figura  8. 

El  paquete  dynamics  dispone  de  otra  rutina  para  el  método  de  Runge-Kutta,  rk,  que 
permite  la  resolución  de  sistemas  de  ecuaciones  difereciales.  Para  resolver  el  sistema 

í  |  =  4 -x2  — 4y2 

1  t  =  y2-*2  +  l 

cargamos  el  paquete  y  ejecutamos  la  sentencia  correspondiente,  junto  con  el  gráfico  aso¬ 
ciado  que  no  mostramos. 

(°/0i37)  load(dynamics) $ 

(0/„i38)  rk([4-x-'2-4*y~2,y~2-x-'2+l]  ,  [x, y]  ,  [-1 .25,0.75]  ,  [t  ,0,4,0.02]  )$ 

(0/„i39)  draw3d ( 

points_joined  =  true. 
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point_type  =  dot, 

points(%) , 

terminal  =  eps)$ 

Para  más  información  sobre  la  sintaxis  de  esta  función,  ejecútese  ?  rk.  El  paquete 
dynamics  contiene  otras  funciones  relacionadas  con  los  sistemas  dinámicos  y  los  fractales. 

5.6.  Vectores  y  campos  vectoriales 

En  Maxima,  los  vectores  se  introducen  como  simples  listas  (Sección  3.3),  siendo  el  caso 
que  con  ellas  se  pueden  realizar  las  operaciones  de  adición,  producto  por  un  número  y 
producto  escalar  de  vectores, 

C/.il)  [l,2,3]  +  [a,b,c]  ; 

(%ol)  [a+ 1,6  +  2, c+3] 

(°/0i2)  s*[a,b,c]; 

( %o2)  [as,  6 s,  es] 

(°/oi3)  [1,2,3] . [a,b,c] ;  /*  producto  escalar  */ 

( %o3)  3  c  +  2  6  +  a 

El  cálculo  del  módulo  de  un  vector  se  puede  hacer  mediante  la  definición  previa  de  una 
función  al  efecto: 

(°/0i4)  modulo  (v):  = 

if  listp(v) 

then  sqrt(apply("+",v~2)) 

else  error ("Mucho  ojito:  ",  v,  "  no  es  un  vector  !!!!")$ 

(°/0i5)  xx:  [a,b,c,d,e]  $ 

(°/oi6)  yy:  [3, 4, -6, 0,4/5]$ 

(°/0i7)  modulo  (xx-yy)  ; 


(%o7)  ^  ^  +  d2  +  (c  +  6)2  +  (6  —  4)2  +  (a  —  3)2 

Los  operadores  diferenciales  que  son  de  uso  común  en  el  ámbito  de  los  campos  vec¬ 
toriales  están  programados  en  el  paquete  vect,  lo  que  implica  que  debe  ser  cargado  en 
memoria  antes  de  ser  utilizado.  Sigue  a  continuación  una  sesión  de  ejemplo  sobre  cómo 
usarlo. 

Partamos  de  los  campos  escalares  c/)(x,  y ,  z)  =  —x+y2+z2  y  ip(x,  y ,  z)  =  4x+log(y2+2:2) 
y  demostremos  que  sus  superficies  de  nivel  son  ortogonales  probando  que  V</>  •  VV’  =  0, 
siendo  V  el  operador  gradiente, 
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(°/0i8)  /*  Se  carga  el  paquete  */ 
load(vect)$ 

(°/0i9)  /*  Se  definen  los  campos  escalares  */ 
phi :  y~2+z~2-x$  psi : log(y~2+z~2)+4*x$ 

(lili)  grad(phi)  .  grad(psi) ; 

( %oll)  grad  ( z 2  +  y2  —  x)  ■  grad  (log  (z2  +  y2)  +  Ax) 

Como  se  ve,  Maxima  se  limita  a  devolvernos  la  misma  expresión  que  le  introducimos; 
el  estilo  de  trabajo  del  paquete  vect  requiere  el  uso  de  dos  funciones:  express  y  ev,  la 
primera  para  obtener  la  expresión  anterior  en  términos  de  derivadas  y  la  segunda  para 
forzar  el  uso  de  éstas. 

(°/0il2)  express  ('/,)  ; 


£  (z2  +  y2-x)  (á  (lo§  ( 2 2  +  y2)  +  4a0)  + 

(%ol2)  ¿  (; z2  +  y2-x )  (log  (z2  +  y2)  +  4*))  + 

£(z2  +  y2-*)  (£  (l°g  (-Z2  +  y2)  +  4r)) 

(°/0il3)  ev(°/0,diff ) ; 


( %ol3) 

(°/0il4)  ratsimp(7J  ; 


4  y2 


(%ol4)  0 

Al  final,  hemos  tenido  que  ayudar  un  poco  a  Maxima  para  que  terminase  de  reducir 
la  última  expresión. 

Sea  ahora  el  campo  vectorial  definido  por  F  =  xyi+x2zj—ex+y'k  y  pidámosle  a  Maxima 
que  calcule  su  divergencia,  (V  •  F),  rotacional  (V  x  F)y  laplaciano  (V2F) 

(°/0il5)  F:  [x*y,x~2*z,exp(x+y)] $ 

(°/0il6)  div  (F) ;  /*  divergencia  */ 


div  [[xy,x2  z,ey+x]) 


(°/0il7)  express  (°/0)  ; 


A 

dy 


{x2 


+  f  ey~x+  (xy) 
dz  dx 


( %o!7) 
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(°/0il8)  ev  (%,  diff )  ; 


(%ol8)  y 

(°/0il9)  curl  (F) ;  /*  rotacional  */ 

(%ol9)  curl  ([xy,x2  z,ey+x]) 

(°/0i20)  express  (°/0)  ; 


( %o20) 

(°/0i21)  ev  (°/„,  diff); 


(%o21)  [ey+x  -x2,-ey+x,2xz-x\ 

(°/0i22)  laplacian  (F) ;  /*  laplaciano  */ 

( %o22)  laplacian  ( [x  y,  x2  z,  ev+x ] ) 

(°/0i23)  express  (°/0)  ; 

( %o23)  ^2  [xy,x2  z,ey+x]  +  [X1L  z->  ey+x]  +  [xy,x2  z,ey+x] 

(°/0i24)  ev  (% ,  diff); 


(%o24)  [0,2*,  2  e***] 

Nótese  en  todos  los  casos  el  usos  de  la  secuencia  express  -  ev.  Si  el  usuario  encuentra 
incómoda  esta  forma  de  operar,  siempre  podrá  definir  una  función  que  automatice  el 
proceso;  por  ejemplo, 

(°/0i25)  milaplaciano(v)  :=  ev(express (laplacian(v) )  ,  diff)  $ 

(°/0i26)  milaplaciano(F) ; 


(%o26)  [0 ,2z,2ey+x\ 

Por  último,  el  paquete  vect  incluye  también  la  definición  del  producto  vectorial,  al 
cual  le  asigna  el  operador  ", 

(°/.i27)  [a,  b,  c]  ~  [x,  y,  z]  ; 


[a,  b,  c ]  ~  [x,  y,  z\ 


(%o27) 

(°/0i28)  express  (*/,); 
(%o28) 


[bz  —  cy, ex  —  az, ay  —  bx\ 
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6.  Análisis  de  datos 

6.1.  Probabilidad 

El  paquete  distrib  contiene  la  definición  de  las  funciones  de  distribución  de  probabilidad 
más  comunes,  tanto  discretas  (binomial,  de  Poisson,  de  Bernoulli,  geométrica,  uniforme 
discreta,  hipergeométrica  y  binomial  negativa),  como  continuas  (normal,  t  de  Student, 
X2  de  Pearson,  F  de  Snedecor,  exponencial,  lognormal,  gamma,  beta,  uniforme  continua, 
logística,  de  Pareto,  de  Weibull,  de  Rayleigh,  de  Laplace,  de  Cauchy  y  de  Gumbel).  Para 
cada  una  de  ellas  se  puede  calcular  la  probabilidad  acumulada,  la  función  de  densidad,  los 
cuantiles,  medias,  varianzas  y  los  coeficientes  de  asimetría  y  curtosis: 

C7.il)  load(distrib)$ 

(°/0i2)  assume(s>0)$ 


(°/0i3)  cdf_normal(x,mu,s) ; 


( %o3) 


2 


1 


(°/0i4)  pdf_poisson(5, 1/s) ; 


(%o4) 


e  s 


120  s5 


(°/0i5)  quantile_student_t  (0.05, 25)  ; 

(%o5)  -1.708140543186975 

(°/0i6)  mean_weibull(3,67) ; 


( %o6) 

(°/0i7)  var_binomial(34, 1/8)  ; 


67 PQ) 

3 


119 

"32" 


( %o7) 
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(°/0i8)  skewness_rayleigh(l/5) ; 

ZlÍ  _  ÉÉt 

( %o8)  ^ - V 

(°/0i9)  kurtosis_gumbel  (2,3); 

(%o9)  y 

Si  su  argumento  es  un  número  entero  positivo,  la  función  random(n)  genera  un  número 
seudoaleatorio  con  distribución  uniforme  discreta  entre  Oyn-1,  ambos  inclusive;  así,  una 
simulación  del  lanzamiento  de  un  dado  sería 

(°/0ilO)  random(6)+l ; 

( %ol0)  3 

y  una  serie  de  100  lanzamientos  de  una  moneda, 

(°/0ill)  makelist  (random(2) ,  i ,  1 , 100)  ; 

[0,0, 1,0, 1,1, 1,0, 1,1, 1,0, 0,1, 0,1,  0,1,  0,1, 1,1, 1,0, 0,0, 1,0, 1,0, 1,1,0, 
(%oll)  0,1, 1,1, 1,0, 0,0, 1,0, 1,0, 1,0, 1,0,  0,1, 0,0, 0,0, 0,1, 0,0, 1,1,  0,0,  0,0,  0,0, 

1,1, 1,0, 1,1, 1,0, 1,0, 0,0, 0,1, 1,0, 0,1, 0,1, 0,1, 1,0, 0,1, 1,0, 0,1, 0,0,1] 

Cuando  el  argumento  es  un  número  decimal  positivo,  la  variable  aleatoria  que  se  simula 
es  la  uniforme  continua,  dando  como  resultado  un  número  real  perteneciente  al  intervalo 
[0,r), 

(°/0il2)  random(6 . 0)  ; 

( %ol2)  0.373047098775396 

El  algoritmo  generador  de  los  números  seudoaleatorios  es  determinista,  de  manera  que 
partiendo  de  una  misma  semilla  o  valor  inicial,  se  generará  la  misma  secuencia  de  números. 
Para  controlar  el  valor  de  esta  semilla  disponemos  de  las  funciones  make_random_state  y 
set_random_state;  por  ejemplo,  para  definir  una  semilla  que  se  genere  a  partir  del  estado 
actual  del  reloj  del  sistema  haremos 

(0/oil3)  nueva_semilla:  make_random_state(true)$ 

Sin  embargo,  para  que  tal  semilla  se  active  en  el  generador,  debemos  indicarlos  expre¬ 
samente  haciendo 
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(°/0il4)  set_random_state (nueva_semilla) $ 

El  argumento  de  la  función  make_random_state  puede  ser  también  un  número  entero, 
como  se  hace  en  el  ejemplo  de  más  abajo. 

Veamos  un  caso  de  aplicación  de  todo  esto.  Supongamos  que  queremos  simular  dife¬ 
rentes  series  estocásticas,  pero  que  todas  ellas  sean  iguales.  Si  hacemos 

(°/0il5)  makelist (random(6) , i,  1,10)  ; 

(%ol5)  [3, 0,0,  5, 0,5, 1,4, 4, 5] 

(°/0il6)  makelist (random(6) , i,  1,10)  ; 

(%ol6)  [0,4,  0,0,  5, 4, 0,0, 1,5] 

(°/0il7)  makelist  (random(6) ,  i ,  1 , 10)  ; 

(%ol7)  [3, 3,  3, 1,3, 2, 1,5, 2, 4] 

lo  más  probable  es  que  obtengamos  tres  secuencias  distintas,  como  en  el  ejemplo.  Pero  si 
hacemos 

(°/0il8)  semilla:  make_random_state(  123456789)$ 

(°/0il9)  set_random_state (semilla) $  makelist  (random(6) , i ,  1 , 10) ; 

(%ol9)  [4, 4, 0,1, 0,3, 2, 5, 4, 4] 

(°/0i20)  set_random_state (semilla) $  makelist  (random(6) , i ,  1 , 10) ; 

(%o20)  [4, 4,  0,1, 0,3, 2, 5, 4, 4] 

(°/0i21)  set_random_state (semilla) $  makelist  (random(6) , i ,  1 , 10) ; 

(%o2l)  [4, 4,  0,1, 0,3, 2, 5, 4, 4] 

se  verá  que  las  tres  secuencias  son  iguales,  ya  que  antes  de  generar  cada  muestra  aleatoria 
reiniciamos  el  estado  del  generador.  La  función  random  y  las  otras  funciones  relacionadas 
con  ella  discutidas  hasta  aquí  están  disponibles  sin  necesidad  de  cargar  el  paquete  distrib. 

Sin  embargo,  el  paquete  adicional  distrib  también  permite  simular  muchas  otras  va¬ 
riables  aleatorias,  tanto  discretas  como  continuas.  A  modo  de  ejemplo,  pedimos  sendas 
muestra  de  tamaño  5  de  las  variables  aleatorias  binomial  5(5,  |),  Poisson  5(7),  hiper- 
geométrica  55(15,20,7),  exponencial  Exp(  12.5)  y  Weibull  Wei( 7,  ^);  puesto  que  ya  se 
ha  cargado  más  arriba  el  paquete,  no  es  necesario  ejecutar  nuevamente  load  (distrib) . 

(°/0i22)  random_binomial  (5 , 1/3 , 5)  ; 
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(%o22)  [3, 1,2, 3, 2] 

(°/0i23)  random_poisson(7,5)  ; 

( %o23)  [8, 3, 10, 5, 6] 

(°/0i24)  random_hypergeometric  (15,20,7,5) ; 


(%o24)  [4, 2, 4, 3, 2] 

(°/0i25)  random_exp(12.5,5) ; 


( %o25) 


[.05865376074017901,  .2604319923173137,  .07552948674579418, 
.02948508382731128,  .2117111885482312] 


C/.Í26)  random.weibull  (7 , 23/3 , 5)  ; 


[6.35737206358163,7.436207845095266,8.101343432607079, 

'  u  ‘  7.835164678709573, 6.350884234996046] 

Para  más  información  sobre  estas  y  otras  funciones  de  simulación  estocástica,  tecléese 

?  distrib. 

6.2.  Estadística  descriptiva 

Maxima  es  capaz  de  realizar  ciertos  tipos  de  procesamiento  de  datos.  El  paquete 
descriptive,  junto  con  otros,  es  el  que  nos  puede  ayudar  en  este  tipo  de  tareas. 

Durante  la  instalación  de  Maxima  se  almacenan  junto  al  fichero  descriptive  .mac  tres 
muestras  de  datos:  pidigits . data,  wind.data  y  biomed.data,  los  cuales  cargaremos  en 
memoria  con  las  funciones  read_list  y  read_matrix. 

Las  muestras  univariantes  deben  guardarse  en  listas,  tal  como  se  muestra  a  continua¬ 
ción 

(°/.il)  si:  [3, 1,4, 1,5, 9, 2, 6, 5, 3, 5]  ; 


(%ol)  [3, 1,4, 1,5, 9,  2, 6, 5, 3,  5] 

y  muestras  multivariantes  en  matrices,  como  en 

(°/0i2)  s2:matrix(  [13. 17,  9 . 29]  ,  [14 . 71 ,  16.88]  ,  [18.50,  16.88], 
[10.58,  6.63] ,  [13.33,  13.25] ,  [13.21,  8.12]); 
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/13.17 

9.29  \ 

14.71 

16.88 

18.5 

16.88 

10.58 

6.63 

13.33 

13.25 

\13.21 

8.12  / 

En  este  caso,  el  número  de  columnas  es  igual  a  la  dimensión  de  la  variable  aleatoria  y 
el  número  de  filas  al  tamaño  muestral. 

Los  datos  se  pueden  introducir  manualmente,  pero  las  muestras  grandes  se  suelen 
guardar  en  ficheros  de  texto.  Por  ejemplo,  el  fichero  pidigits .  data  contiene  los  100 
primeros  dígitos  del  número  7r: 

3 
1 

4 
1 

5 
9 
2 

6 
5 

3  ... 

Para  cargar  estos  dígitos  en  Maxima, 

(°/0i3)  si  :  read_list  (file_search  ("pidigits .data") )$ 

(°/0i4)  length(sl) ; 


( %o4)  100 

El  fichero  wind.  data  contiene  las  velocidades  medias  del  viento  registradas  diariamente 
en  5  estaciones  meteorológicas  de  Irlanda.  Lo  que  sigue  carga  los  datos, 

(°/0i5)  s2:read_matrix(file_search  ("wind. data") )$ 

(°/0i6)  length(s2) ; 

( %o6)  100 

(°/0i7)  s2  [°/0]  ;  /*  último  registro  */ 

( %o7)  [3.58, 6.0, 4.58, 7.62, 11.25] 

Algunas  muestras  incluyen  datos  no  numéricos.  Como  ejemplo,  el  fichero  biomed. data 
contiene  cuatro  medidas  sanguíneas  tomadas  en  dos  grupos  de  pacientes,  A  y  B,  de  edades 
diferentes, 
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(°/0i8)  s3:read_matrix(file_search  ("biomed.data"))$ 
(°/0i9)  length(s3) ; 


( %o9)  100 

(°/0ilO)  s3[l];  /*  primer  registro  */ 

( %ol0)  [A,  30, 167.0, 89.0, 25.6, 364] 

El  primer  individuo  pertenece  al  grupo  A,  tiene  30  años  de  edad  y  sus  cuatro  medidas 
sanguíneas  fueron  167.0,  89.0,  25.6  y  364. 

El  paquete  descriptive  incluye  dos  funciones  que  permiten  hacer  recuentos  de  datos: 
continuous_freq  y  discrete_f req;  la  primera  de  ellas  agrupa  en  intervalos  los  valores 
de  una  lista  de  datos  y  hace  el  recuento  de  cuántos  hay  dentro  de  cada  una  de  las  clases, 

C7.Í11)  load("descriptive")$ 

(°/0il2)  continuous_f req(sl , 5) ; 

( %ol2)  [[0, 1.8, 3.6, 5.4,  7.2, 9.0] ,  [16,  24, 18, 17,  25]] 

La  primera  lista  contiene  los  límites  de  los  intervalos  de  clase  y  la  segunda  los  recuentos 
correspondientes:  hay  16  dígitos  dentro  del  intervalo  [0, 1.8],  eso  es  ceros  y  unos,  24  dígitos 
en  (1.8, 3.6],  es  decir,  doses  y  treses,  etc.  El  segundo  parámetro  indica  el  número  de  clases 
deseadas,  que  por  defecto  es  10. 

La  función  discrete_freq  hace  el  recuento  de  las  frecuencias  absolutas  en  muestras 
discretas,  tanto  numéricas  como  categóricas.  Su  único  argumento  es  una  lista, 

(°/0il3)  discrete_freq(sl) ; 

( %ol3)  [[0, 1, 2, 3, 4, 5, 6, 7, 8, 9] ,  [8, 8, 12, 12, 10, 8, 9, 8, 12, 13]] 

(°/0il4)  discrete_freq(transpose (col (s3 , 1)  )  [1] ) ; 

(%ol4)  [[A,  B] ,  [35,  65]] 

La  primera  lista  contiene  los  valores  muéstrales  y  la  segunda  sus  frecuencias  absolutas. 

Las  instrucciones  ?  col  y  ?  transpose  ayudarán  a  comprender  la  última  entrada. 

En  cuanto  al  cálculo  de  parámetros  muéstrales,  el  paquete  descriptive  incluye,  tanto 
para  muestras  univariantes  como  multivariantes,  medias  (mean),  varianzas  (var),  desviacio¬ 
nes  típicas  (std),  momentos  centrales  (central_moment)  y  no  centrales  (noncentral_moment), 
coeficientes  de  variación  (cv),  rangos  (range),  medianas  (median),  cuantiles  (quantile), 
coeficientes  de  curtosis  (kurtosis)  y  de  asimetría  (skewness),  y  otros.  En  los  siguientes 
ejemplos,  si  es  una  muestra  univariante  y  s2  multivariante: 
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(°/0il5)  mean(sl) ; 
( %ol5) 

(°/0il6)  f  loat  (%) ; 


471 

100 


(%ol6)  4.71 

(°/0il7)  mean(s2) ;  /*  vector  de  medias  */ 

( %ol7)  [9.9485, 10.1607, 10.8685, 15.7166, 14.8441] 

Cuando  la  muestra  contiene  datos  en  formato  decimal,  los  resultados  serán  también 
en  este  formato,  pero  si  la  muestra  contiene  tan  solo  datos  enteros,  los  resultados  serán 
devueltos  en  formato  racional.  En  caso  de  querer  forzar  que  los  resultados  sean  devueltos 
siempre  en  formato  decimal,  podemos  jugar  con  la  variable  global  numer. 

(°/0il8)  numer:  true$ 

(°/0il9)  var  (si)  ; 

( %ol9)  8.425899999999999 

(°/0i20)  var(s2);  /*  vector  de  varianzas  */ 

( %o20) 

[17.22190675000001, 14.98773651, 15.47572875, 32.17651044000001, 24.42307619000001] 

(°/0i21)  /*  1er  y  3er  cuartiles  */ 

[quantile(sl , 1/4) ,quantile(sl,3/4)] ; 

(%o21)  [2.0,7.25] 

(°/0i22)  range(sl) ; 

( %o22)  9 

(°/0i25)  range(s2);  /*  vector  de  rangos  */ 

( %o25)  [19.67,  20.96, 17.37, 24.38, 22.46] 

(°/0i26)  kurtosis  (si)  ; 


(%o26) 


-1.273247946514421 


6  ANÁLISIS  DE  DATOS 


97 


(°/0i27)  kurtosis (s2) ;  /*  vector  de  coef .  curtosis  */ 

s  [-.2715445622195385,0.119998784429451,  -.4275233490482866, 

^  /o°  ^  -.6405361979019522,  -.4952382132352935] 

Un  aspecto  crucial  en  el  análisis  estadístico  multivariante  es  la  dependencia  estocástica 
entre  variables,  para  lo  cual  Maxima  permite  el  cálculo  de  matrices  de  covarianzas  (cov) 
y  correlaciones  (cor),  así  como  de  otras  medidas  de  variabilidad  global.  La  variable  global 
fpprintprec  es  utilizada  para  limitar  el  número  de  decimales  a  imprimir;  lo  cual  no  afecta 
a  la  precisión  de  los  cálculos. 

(°/0i28)  fpprintprec: 7$  /*  número  de  dígitos  deseados  en  la  salida  */ 

(°/0i29)  cov(s2) ;  /*  matriz  de  covarianzas  */ 


/17.22191 

13.61811 

14.37217 

19.39624 

15.42162N 

13.61811 

14.98774 

13.30448 

15.15834 

14.9711 

( %o29) 

14.37217 

13.30448 

15.47573 

17.32544 

16.18171 

19.39624 

15.15834 

17.32544 

32.17651 

20.44685 

\15.42162 

14.9711 

16.18171 

20.44685 

24.42308/ 

(°/.i30) 

cor (s2) ; 

/*  matriz  de  correlaciones 

*/ 

/  1.0 

.8476339 

.8803515 

.8239624 

.7519506X 

.8476339 

1.0 

.8735834 

.6902622 

0.782502 

( %o30) 

.8803515 

.8735834 

1.0 

.7764065 

.8323358 

.8239624 

.6902622 

.7764065 

1.0 

.7293848 

\.  75 19506 

0.782502 

.8323358 

.7293848 

1.0  ) 

Las  funciones  global_variances  y  list_correlations  calculan,  respectivamente, 
ciertas  medidas  globales  de  variación  y  correlación;  para  más  información,  pídase  ayuda  a 
Maxima  con  el  comando  describe  o  con  ?. 

Hay  varios  tipos  de  diagramas  estadísticos  programados  en  el  paquete  descriptive: 
scatterplot  para  gráficos  de  dispersión,  histogram,  barsplot,  piechart  y  boxplot  para 
el  análisis  homocedástico.  A  continuación  se  presentan  un  par  de  ejemplos  cuyos  salidas 
se  presentan  en  la  Figura  9. 

(°/0i31)  /*  Comparando  variabilidades  */ 
boxplot ( 
s2 , 

terminal  =  eps, 

title  =  "Velocidades  del  viento")$ 

(°/0i32)  /*  Diagramas  de  dispersión  para  las 

tres  últimas  estaciones  meteorológicas  */ 
scatterplot ( 

submatrix(s2, 1 ,2) , 
point_type  =  circle, 
terminal  =  eps)$ 
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Velocidades  del  viento 


Figura  9:  Gráficos  para  muestras  multivariantes:  a)  diagramas  de  cajas;  b)  diagramas  de 
dispersión  para  tres  variables. 
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6.3.  Estadística  inferencial 

El  paquete  stats  contiene  algunos  procedimientos  clásicos  de  estadística  inferencial. 
Cada  uno  de  estos  procedimientos  devuelve  un  objeto  de  tipo  inf erence_result,  el  cual 
guarda  todos  los  resultados  obtenidos,  aunque  sólo  muestre  una  selección  de  ellos  por 
defecto,  los  que  se  consideran  de  más  relevancia;  el  resto  permanecen  ocultos  y  sólo  serán 
visibles  a  petición  del  usuario.  Los  siguientes  ejemplos  aclararán  su  uso. 

Partamos  de  una  muestra  de  tamaño  20  y  tratemos  de  contrastar  la  hipótesis  nula  de 
que  la  media  de  la  población  de  la  cual  procede  es  igual  a  32,  frente  a  la  alternativa  de 
que  la  media  es  diferente  a  este  valor.  Para  aplicar  el  test  de  Student  nos  interesa  saber 
si  podemos  considerar  la  población  normal;  la  función  test_normality  nos  puede  ser  de 
ayuda, 

C7.il)  load(stats)$ 

(°/0i2)  x:  [28.9,35.0,30.6,31.3,31.9,31.7,29.3,37.8,29.9,36.8, 
28.0,25.3,39.6,30.4,23.3,24.7,38.6,31.3,29.8,31.9]$ 

(°/0i3)  test_normality(x) ; 


( %o3) 


SHAPIRO  -  WILK  TEST 
statistic  =  0.9508091 
p.value  =  0.3795395 


Según  este  resultado,  no  hay  evidencias  suficientes  para  rechazar  la  hipótesis  de  nor¬ 
malidad.  Hagamos  ahora  el  contraste  sobre  la  media, 


(°/0i4)  z:  test_mean(x,mean=32) ; 


( %o4) 


MEAN  TEST 
mean.estimate  =  31.305 
confJevel  =  0.95 

conf  .interval  =  [29.21482, 33.39518] 
method  =  Exact  t-test.  Unknown  variance. 
hypotheses  =  H0:  mean  =  32  ,  Hl:  mean  ^  32 
statistic  =  0.6959443 
distribution  =  [student.t,  19] 

P-value  =  0.4948895 


En  la  llamada  a  la  función  test_mean  se  ha  indicado  mediante  la  opción  mean  el  valor 
de  la  media  en  la  hipótesis  nula.  El  objeto  devuelto  por  la  función  consta,  como  se  ve,  de 
varios  elementos:  el  título  del  objeto,  la  media  muestral,  el  nivel  de  confianza,  el  intervalo 
de  confianza,  información  sobre  el  método  utilizado,  hipótesis  a  contrastar,  estadístico  de 
contraste,  su  distribución  y  el  p- valor  correspondiente.  Si  ahora  quisiésemos  extraer  la 
media  muestral  para  utilizar  en  otros  cálculos, 


(°/0i5)  take_inference(mean_estimate,z) ; 
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(%o5)  31.305 

En  caso  de  querer  extraer  más  de  un  resultado, 

(°/0i6)  take_inf  erence  ( [p_value ,  statistic]  ,  z) ; 

( %o6)  [0.4948895, 0.6959443] 

Para  saber  con  certeza  qué  resultados  guarda  el  objeto  devuelto  liaremos  uso  de  la 
función  items_inference,  con  lo  que  comprobamos  que  la  función  test_mean  no  oculta 
ninguno  de  sus  resultados, 

(°/0i7)  items_inf  erence  (z) ; 

( %o7) 

[mean. estímate,  confJevel,  conf  .interval ,  method ,  hypotheses,  statistic,  distribution,p. valué] 

En  caso  de  que  la  muestra  no  procediese  de  una  población  normal,  y  si  el  tamaño 
muestral  fuese  suficientemente  grande,  se  podría  hacer  uso  de  la  opción  asymptotic=true, 
si  la  desviación  típica  de  la  población  fuese  conocida,  por  ejemplo  a  =  4.5,  podría  añadirse 
la  opción  dev=4.5.  En  fin,  las  opciones  que  admite  esta  función  permiten  al  usuario  tener 
cierto  control  sobre  el  procedimiento;  para  más  información  consúltese  el  manual. 

Si  la  población  de  referencia  no  es  normal  y  la  muestra  es  pequeña,  siempre  se  puede 
acudir  a  un  test  no  paramétrico  para  realizar  un  contraste  sobre  la  mediana, 

(°/0i8)  signed_rank_test(x,median=32)  ; 


SIGNED  RANK  TEST 
med.estimate  =  30.95 
method  =  Asymptotic  test.  Ties 
hypotheses  =  H0:  med  =  32  ,  Hl:  med  ^  32 
statistic  =  75 

distribution  =  [normal,  104.5,  26.78152] 
p.value  =  0.2706766 


Existe  en  Maxima  la  posibilidad  de  ajustar  por  mínimos  cuadrados  a  datos  empíricos 
los  parámetros  de  curvas  arbitrarias  mediante  las  funciones  del  paquete  lsquares.  Vamos 
a  simular  una  muestra  de  tamaño  100  de  valores  x  en  el  intervalo  [0, 10].  Con  estas  abscisas 
calculamos  después  las  ordenadas  a  partir  de  la  suma  de  una  señal  determinista  (/)  más  un 
ruido  gaussiano.  Para  hacer  la  simulación  necesitamos  cargar  en  primer  lugar  al  paquete 
distrib. 


(°/0i9)  load  (distrib)  $ 

(°/0ilO)  abs:  makelist (random_continuous_unif orm(0 , 10)  ,k,  1, 100)$ 

(°/„ill)  f  (x)  :=3*(x-5)~2$ 

(°/0il2)  data:  apply (matrix, makelist ( [x,  f  (x) +random_normal (0,1)3  ,x, abs)) $ 
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Aunque  la  señal  determinista  obedece  a  una  función  polinómica  de  segundo  grado,  si 
no  lo  sabemos  a  priori  intentamos  ajustar  un  polinomio  cúbico: 

(°/0il3)  load(lsquares) $ 

(°/0il4)  param:  lsquares_estimates(data,  [x,y]  , 

y=a*x~3+b*x"'2+c*x+d,  [a,b,c,d]),  numer; 


(%ol4) 


[[a  =  —0.002705800223881305,  b  =  3.018798873646606  , 
c  =  -29.94151342602112,  d  =  74.78603431944423]] 


Vemos  claramente  que  el  término  de  tercer  grado  es  supérfluo,  por  lo  que  reajustamos 
al  de  segundo  grado, 


(°/0il5)  param:  lsquares_estimates  (data,  [x,y]  ,y=b*x~2+c*x+d,  [b,  c,d] ), numer ; 


( %ol5)  [[ b  =  2.979110687882263,  c  =  -29.78353057922009,  d  =  74.64523259993118]] 

Ahora  estamos  en  disposición  de  estudiar  los  residuos  para  valorar  el  ajuste.  En  primer 
lugar  calculamos  el  error  cuadrático  medio  del  residuo,  definido  como 

MSE  =  ^  ¿  (rhs(e¿)  -  lhs(e¿))2  , 

siendo  n  el  tamaño  de  la  muestra  y  rhs(e¿)  y  lhs(e¿)  los  miembros  derecho  e  izquierdo, 
respectivamente,  de  la  expresión  a  ajustar.  Para  el  caso  que  nos  ocupa,  el  valor  calculado 
es 


(°/0il6)  lsquares_residual_mse(data,  [x,y]  ,  y=b*x~2+c*x+d,  f ir st (param) )  ; 


(%ol6)  1.144872557335554 

También  podemos  calcular  el  vector  de  los  residuos,  definidos  como  lhs(e¿)  —  rhs(e¿),  Vi, 
para  a  continuación  representarlos  gráficamente  junto  con  los  datos  y  la  curva  ajustada, 
tal  como  se  aprecia  en  la  Figura  10. 

(°/0il7)  res:  lsquares_residuals(data,  [x,y]  ,  y=b*x~2+c*x+d,  f  ir  st  (param) )  $ 
(°/0il8)  load(draw)$  /*  necesitaremos  las  rutinas  gráficas  */ 

(°/0il9)  scenel:  gr2d(points(data) , 

explicit (ev(b*x~2+c*x+d,f irst (param) ) ,x,0 , 10) )$ 

(°/0i20)  scene2:  gr2d(points_joined  =  true, 
points(res) )$ 

(°/0i21)  draw(terminal  =  eps,  scenel,  scene2)$ 
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Figura  10:  Resultados  gráficos  del  ajuste  por  mínimos  cuadrados. 


Nosotros  ya  sabemos  que  el  ruido  del  proceso  simulado  es  gaussiano  de  esperanza  nula. 
En  una  situación  real,  este  dato  lo  desconoceremos  y  nos  interesará  contrastar  la  hipótesis 
de  normalidad  y  la  de  la  nulidad  de  la  media  de  los  residuos.  El  paquete  stats  tiene 
algunos  procedimientos  inferenciales,  entre  los  cuales  están  las  funciones  test_normality 
y  test_mean;  primero  cargamos  el  paquete  y  contrastamos  la  hipótesis  de  normalidad, 

(°/0i22)  load(stats)$ 

(°/0i23)  test_normality  (res)  ; 

SHAPIRO  -  WILK  TEST 

( %o31)  statistic  =  0.986785359052448 

p. valué  =  0.423354600782769 

El  p- valor  es  lo  suficientemente  alto  como  para  no  rechazar  la  hipótesis  nula  de  nor¬ 
malidad;  ahora  contrastamos  la  nulidad  de  la  esperanza  del  ruido, 


(°/0i24)  test_mean(res) ; 
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(%o32) 


MEAN  TEST 

mean-estimate  =  —1.002451472320587  x  10-7 
confJevel  =  0.95 

conf  interval  =  [-.2133782738324136,-2133780733421191] 
method  =  Exact  t-test.  Unknown  variance. 
hypotheses  =  HO:  mean  =  0  ,  Hl:  mean  ^  0 
statistic  =  9.321855844715968  x  lO"7 


distribution  =  [student-t,  99] 
pjvalue  =  .9999992583830712 


Sin  duda  admitiríamos  la  hipótesis  nula  como  válida. 

Consúltese  el  manual  de  referencia  de  Maxima  para  ver  qué  otros  tests  están  progra¬ 
mados  en  el  paquete  stats,  así  como  las  opciones  que  se  pueden  utilizar  en  cada  caso. 


6.4.  Interpolación 

El  paquete  interpol  permite  abordar  el  problema  de  la  interpolación  desde  tres  enfo¬ 
ques:  lineal,  polinomio  de  Lagrange  y  splines  cúbicos. 

A  lo  largo  de  esta  sección  vamos  a  suponer  que  disponemos  de  los  valores  empíricos  de 
la  siguiente  tabla: 

x  II  7  I  8  I  1  I  3  I  6 
y  '  2  2  5  2  7 

Nos  planteamos  en  primer  lugar  el  cálculo  de  la  función  de  interpolación  lineal,  para 
lo  cual  haremos  uso  de  la  función  linearinterpol, 

C/.il)  load(interpol) $ 

C/.Í2)  datos:  [[7,2]  ,  [8,2]  ,  [1,5]  ,  [3,2]  ,  [6,7]]$ 


(°/0i3)  linearinterpol  (datos)  ; 

(0¡  .  (t  —  if)  charfuns  (x,  —  oo,  3)  +  2  charfuns  (x,  7,  oo)  + 

(37  —  5  x)  charfun2  (x,  6,  7)  +  —  3)  charfun2  (x,  3, 6) 

(°/oi4)  f(x)  :  =  ”•/.$ 

Empezamos  cargando  el  paquete  que  define  las  funciones  de  interpolación  y  a  conti¬ 
nuación  introducimos  los  pares  de  datos  en  forma  de  lista.  La  función  linearinterpol 
devuelve  una  expresión  definida  a  trozos,  en  la  que  charfun2(x,a,b)  devuelve  1  si  el 
primer  argumento  pertence  al  intervalo  [a,  b)  y  0  en  caso  contrario.  Por  último,  defini¬ 
mos  cierta  función  f  previa  evaluación  (dos  comillas  simples)  de  la  expresión  devuelta  por 
linearinterpol.  Esta  función  la  podemos  utilizar  ahora  tanto  para  interpolar  como  para 
extrapolar: 
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(0/.i5)  map  (f ,  [7.3, 25/7 ,  °/0pi]  )  ; 


(%o5) 

(0/.i6)  f  loat  C/J  ; 


'  62  5? r 

2’21’T 


( %o6)  [2.0, 2.952380952380953, 2.235987755982989] 

Unos  comentarios  antes  de  continuar.  Los  datos  los  hemos  introducido  como  una  lista 
de  pares  de  números,  pero  también  la  función  admite  una  matriz  de  dos  columnas  o  una 
lista  de  números,  asignándole  en  este  último  caso  las  abscisas  secuencialmente  a  partir  de 
la  unidad;  además,  la  lista  de  pares  de  la  variable  datos  no  ha  sido  necesario  ordenarla 
respecto  de  la  primera  coordenada,  asunto  del  que  ya  se  encarga  Maxima  por  cuenta 
propia. 

El  polinomio  de  interpolación  de  Lagrange  se  calcula  con  la  función  lagrange;  en  el 
siguiente  ejemplo  le  daremos  a  los  datos  un  formato  matricial  y  le  indicaremos  a  Maxima 
que  nos  devuelva  el  polinomio  con  variable  independiente  w, 

(0/„i7)  datos2 :  matrix(  [7,2]  ,  [8,2]  ,  [1,5]  ,  [3,2]  ,  [6,7] ) ; 


( %o7) 


/7  2\ 
8  2 
1  5 
3  2 

\6  V 


(°/0i8)  lagrange  (datos2 ,  varname=  ’  w)  ; 


73  w* 
420 


701  w3  8957  w2 
210  +  420 


5288  w  186 

105  +  t~ 


C7.19)  g(w)  :  =  ”*/.$ 

(°/0ilO)  map(g,  [7.3,25/7,°/0pi] ) ,  numer; 


( %ol0)  [1.043464999999799, 5.567941928958199, 2.89319655125692] 

Disponemos  en  este  punto  de  dos  funciones  de  interpolación;  representémoslas  gráfica¬ 
mente  junto  con  los  datos  empíricos, 


(lili)  load(draw)$ 
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(°/0il2)  draw2d( 

key  =  "Interpolador  lineal", 

explicit (f (x) ,x,0, 10) , 
line_type  =  dots, 

key  =  "Interpolador  de  Lagrange", 

explicit (g(x) ,x,0, 10) , 

key  =  "Datos  empiricos", 

points (datos) , 

terminal  =  eps)$ 

cuyo  resultado  se  ve  en  el  apartado  a)  de  la  Figura  11. 

El  método  de  los  splines  cúbicos  consiste  en  calcular  polinomios  interpoladores  de 
tercer  grado  entre  dos  puntos  de  referencia  consecutivos,  de  manera  que  sus  derivadas 
cumplan  ciertas  condiciones  que  aseguren  una  curva  sin  cambios  bruscos  de  dirección.  La 
función  que  ahora  necesitamos  es  cspline, 

(°/0il3)  cspline  (datos)  ; 

(1Sf!  -  -Hir  -  W  +  lü)  charfnn2  (*, -oo,  3)  + 

+  +  charfun2(x, 7,po)  + 

+  charlun^x, 6.7)4- 


( %o!3) 


(°/0il4)  si  (x)  : 


”/.$ 


(°/0il5)  map(sl,  [7.3,25/7,°/0pi] ) ,  numer; 


( %ol5)  [1.438224452554664, 3.320503453379974, 2.227405312429507] 

La  función  cspline  admite,  además  de  la  opción  ’  varname  que  ya  se  vió  anterior¬ 
mente,  otras  dos  a  las  que  se  hace  referencia  con  los  símbolos  ’dl  y  ’dn,  que  indican  las 
primeras  derivadas  en  las  abscisas  de  los  extremos;  estos  valores  establecen  las  condicio¬ 
nes  de  contorno  y  con  ellas  Maxima  calculará  los  valores  de  las  segundas  derivadas  en 
estos  mismos  puntos  extremos;  en  caso  de  no  suministrarse,  como  en  el  anterior  ejemplo, 
las  segundas  derivadas  se  igualan  a  cero.  En  el  siguiente  ejemplo  hacemos  uso  de  estas 
opciones, 

(°/0il6)  cspline  (datos,  ,varname=’z,dl=l  ,dn=0) ; 


(  613, 


¡Ü  -  l4í»f  +  tÜt  -  Si)  (z, -oo,3)  + 

+  Mgi  _  Mgp  +  llgis)  chaTfuní  (,,  7,  oo)  + 

+  charfunz  (z,  6,7)  + 

4S8¡Í  +  -  T1T  +  T+)  zharfima  (o,  3, 6) 
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b) 


Figura  11:  Interpolación:  a)  lineal  y  de  Lagrange;  b)  Splines  cúbicos. 


C/.Í17)  s2(z)  :  =  ’  ,0/0$ 

(°/0il8)  map(s2,  [7.3,25/7//0pi] ) ,  numer; 


( %ol8)  [1.595228723404261, 2.88141531519733, 2.076658794432369] 

Con  esto  hemos  obtenido  dos  interpoladores  distintos  por  el  método  de  los  splines  cúbi¬ 
cos;  con  el  siguiente  código  pedimos  su  representación  gráfica,  cuyo  resultado  se  observa 
en  el  apartado  b)  de  la  Figura  11. 

(°/0il9)  draw2d ( 

key  =  "si", 

explicit (si (x) ,x,0, 10) , 

line_type  =  dots, 

key  =  "s2", 

explicit (s2(x) ,x,0, 10) , 

key  =  "Datos  empíricos", 

points (datos) , 

terminal  =  eps)$ 
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7.  Gráficos 

Maxima  no  está  habilitado  para  realizar  él  mismo  gráficos,  por  lo  que  necesitará  de 
un  programa  externo  que  realice  esta  tarea.  Nosotros  nos  limitaremos  a  ordenar  qué  tipo 
de  gráfico  queremos  y  Maxima  se  encargará  de  comunicárselo  a  la  aplicación  gráfica  que 
esté  activa  en  ese  momento,  que  por  defecto  será  Gnuplot.  La  otra  herramienta  gráfica  es 
Openmath,  un  programa  Tcl-tk  que  se  distribuye  conjuntamente  con  Maxima. 

7.1.  El  módulo  ”plot” 

Las  funciones  plot2d  y  plot3d  son  las  que  se  utilizan  por  defecto.  Veamos  un  ejemplo 
de  cada  una  de  ellas. 

C/.il)  xy:  [[10,  .6]  ,  [20, .9],  [30,1.1],  [40,1.3],  [50,1.4]]$ 

(°/0i2)  plot2d(  [ [discrete ,xy]  ,  2*°/0pi*sqrt (u/980) ]  ,  [u,0,50]  , 

[style ,  [points ,5,2,6] ,  [lines ,1,1]] , 

[legend, "Datos  experimentáis" , "Predicion  da  teoria"] , 

[xlabel, "Lonxitude  do  péndulo  (cm)"],  [ylabel , "periodo  (s) "] , 
[gnuplot_preamble , 

"set  terminal  postscript  eps;set  out  ,plot2d.epsJ "] )$ 

Se  ha  definido  en  el  ejemplo  una  lista  de  puntos  empíricos  a  representar,  junto  con  su 
función  teórica,  pidiendo  su  representación  gráfica  en  el  dominio  [0,  50].  El  resto  de  código 
hace  referencia  a  diversas  opciones:  style,  legend,  xlabel,  ylabel  y  gnuplot_preamble, 
esta  última  permite  escribir  código  de  Gnuplot  para  que  éste  lo  ejecute;  aquí  se  le  pide 
que  devuelva  el  gráfico  en  formato  PostScript.  El  ejemplo  siguiente  genera  una  superficie 
explícita.  Las  dos  salidas  gráficas  se  representan  en  la  Figura  12 

(0/„i3)  plot3d  (2~(-u~2  +  v~2) ,  [u,  -3,  3],  [v,  -2,  2], 

[gnuplot_preamble , 

"set  terminal  postscript  eps;set  out  ’plot3d. eps ’ "] )$ 

El  control  de  las  opciones  gráficas  se  consigue  manipulando  la  variable  global  plot_options, 
cuyo  estado  por  defecto  es 

(°/0i4)  plot_options ; 

(°/.o4)  [[x,  -  1 . 755559702014e+305 ,  1 . 755559702014e+305]  , 

[y,  -  1 . 755559702014e+305 ,  1 . 755559702014e+305] , 

[t,  -  3,  3],  [GRID,  30,  30],  [VIEW_DIRECTI0N,  1,  1,  1], 

[C0L0UR_Z,  FALSE] ,  [TRANSF0RM_XY ,  FALSE] , 

[RUN_VIEWER,  TRUE] ,  [PL0T_F0RMAT,  GNUPLOT] , 

[GNUPL0T_TERM ,  DEFAULT] ,  [GNUPL0T_0UT_FILE,  FALSE] , 

[NTICKS ,  10],  [ADAPT_DEPTH ,  10],  [GNUPL0T.PM3D ,  FALSE], 

[GNUPL0T_PREAMBLE ,  ],  [GNUPL0T_CURVE_TITLES ,  [DEFAULT]], 

[GNUPL0T_CURVE_STYLES ,  [with  lines  3,  with  lines  1, 
with  lines  2,  with  lines  5,  with  lines  4,  with  lines  6, 
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Figura  12:  Gráficos  generados  por  las  funciones  plot2d  y  plot3d. 


with  lines  7]],  [GNUPLOT_DEFAULT_TERM_COMMAND ,  ], 

[GNUPLOT_DUMB_TERM_COMMAND ,  set  term  dumb  79  22], 

[GNUPLOT_PS_TERM_COMMAND ,  set  size  1.5,  1.5; set  term  postsc# 
ript  eps  enhanced  color  solid  24] ] 

Para  mayor  información  sobre  el  significado  de  cada  uno  de  los  elementos  de  esta  lista, 
así  como  de  las  funciones  plot2d  y  plot3d,  consúltese  la  documentación. 

Por  defecto,  Maxima  invocará  al  programa  Gnuplot  para  la  generación  de  gráficos, 
pero  quizás  prefiramos  el  programa  Openmath,  que  forma  parte  de  la  distribución  de 
Maxima;  en  tal  caso  tendríamos  que  modificar  previamente  las  opciones  guardadas  en 
plot_options  y  a  continuación  solicitar  el  gráfico  deseado,  como  en  este  caso  en  el  que  se 
representa  la  función  gamma  y  su  inversa. 

(°/0i5)  set_plot_option(  [plot_f  ormat ,  openmath]  )$ 

(°/0i6)  plot2d(  [gamma(x)  ,l/gamma(x)]  ,  [x,-4.5,5]  ,  [y, -10, 10]  )$ 

El  resto  de  esta  sección  lo  dedicaremos  a  hacer  una  somera  descripción  del  paquete 
draw,  un  proyecto  consistente  en  el  desarrollo  de  un  interfaz  que  permita  aprovechar  al 
máximo  las  habilidades  gráficas  de  Gnuplot. 

7.2.  El  módulo  ”draw” 

Aquí,  las  funciones  a  utilizar  son  draw2d,  draw3d  y  draw,  para  escenas  en  2d,  3d  y 
para  gráficos  múltiples  y  animaciones,  respectivamente.  Puesto  que  se  trata  de  un  módulo 
adicional,  será  necesario  cargarlo  en  memoria  antes  de  hacer  uso  de  él.  Empecemos  por  un 
ejemplo  comentado. 

C7.il)  load  (draw)  $ 

(°/0i2)  draw2d( 

key  =  "Cubic  poly" , 

explicit  C/0pi*x~3+sqrt  (2)  *x~2+10 ,  x ,  0 , 2)  , 
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color  =  blue, 

key  =  "Parametric  curve", 

line_width  =  3, 

nticks  =  50, 

parametric (2*cos(rrr)  +3,  rrr,  rrr,  0,  6*°/0pi) , 

line_type  =  dots, 

points_joined  =  true, 

point_type  =  diamant , 

point_size  =  3, 

color  =  red, 

line_width  =  1 , 

key  =  "Empiric  data" , 

points (makelist (random(40 .0) ,k,l,5)) , 

title  =  "DRAWING  CURVES", 

terminal  =  eps  )$ 

Los  argumentos  de  draw2d  se  dividen  en  tres  grupos:  objetos  gráficos  (en  el  ejem¬ 
plo,  explicit,  parametric  y  points),  opciones  locales  (que  afectan  directamente  a  la 
representación  de  los  objetos  gráficos,  como  key,  color,  line_width,  nticks,  line_type, 
points_joined  y  line_width)  y  opciones  globales  (que  hacen  referencia  a  aspectos  ge¬ 
nerales  del  gráfico,  como  title  y  terminal).  Todos  estos  argumentos  se  interpretan  se- 
cuencialmente,  de  forma  que  al  asignar  un  cierto  valor  a  una  opción  local,  ésta  afectará  a 
todos  los  objetos  gráficos  que  le  sigan.  En  este  ejemplo,  line_width  comienza  teniendo 
valor  1,  que  es  el  asignado  por  defecto,  luego  se  le  da  el  valor  3  para  la  representación 
de  la  curva  paramétrica  y  finalmente  se  le  devuelve  el  valor  original  antes  de  representar 
los  segmentos  que  unen  los  puntos  aleatoriamente  generados.  En  cambio,  las  dos  opciones 
globales,  title  y  terminal,  aunque  se  colocaron  al  final,  podrían  haberse  ubicado  en 
cualquier  otro  lugar. 

Siguiendo  con  los  gráficos  en  dos  dimensiones,  el  siguiente  ejemplo  muestra  una  escena 
en  la  que  intervienen  los  objetos  ellipse,  image,  label,  vector  y,  ya  lo  conocemos, 
explicit.  Antes  de  ejecutar  esta  instrucción  es  necesario  leer  el  fichero  gráfico  gatos  .  xpm9 

(°/0i3)  cats:  read_xpm( "gatos. xpm" )$ 

(°/0i4)  draw2d( 

terminal  =  eps , 
yrange  =  [-4 , 10] , 

ellipse (5 , 3 , 8 , 6 , 0 , 360) , 
image (cats, 0,0, 10,7) , 
line_width  =  2, 
head_length  =  0.3, 
color  =  blue, 

label (["This  is  Francisco", -1,-0. 5] ) , 
vector ([-1,0] ,  [2,4]), 


'El  formato  gráfico  XPM  es  el  único  que  puede  leer  Maxima. 
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color  =  green, 

vector ( [11,7] ,  [-2,-1] ) , 
label( ["This  is  Manolita" , 11 ,8] ) , 
explicit(sin(x)-2,x,-4,15)  )$ 

Junto  con  los  objetos  gráficos  introducidos  en  los  ejemplos,  cuyos  resultados  se  pueden 
ver  en  los  apartados  a)  y  b)  de  la  Figura  13,  también  existen  polygon,  rectangle,  polar, 
implicit  y  geomap,  este  último  para  mapas  cartográficos. 

Mostramos  a  continuación  algunas  escenas  tridimensionales.  En  primer  lugar,  el  valor 
absoluto  de  la  función  F  de  Euler,  junto  con  sus  líneas  de  contorno. 

(°/0i5)  gamma2  (x ,  y) :  = 

block(  [re, im,g: gamma (x+°/0i*y)]  , 
re : realpart (g) , 
im : imagpart (g) , 
sqrt (re~2+im~2) ) $ 

(°/0i6)  draw3d ( 

zrange  =  [0,6] , 

xu_grid  =  50, 

yv_grid  =  50, 

surface_hide  =  true, 

contour  =  sur face, 

contour_levels  =  [0,0. 5, 6] ,  /*  de  0  a  6  en  saltos  de  0.5  */ 

color  =  cyan, 

terminal  =  eps, 

explicit (gamma2 (x , y) , x , -4 , 4 , y , -2 , 2) ) $ 

Una  superficie  paramétrica,  que  junto  con  la  anterior  escena,  se  pueden  ver  ambas  en 
la  Figura  13,  apartados  c)  y  d). 

(°/0i7)  draw3d( 

terminal  =  eps, 

title  =  "Figure  8  -  Klein  bottle", 

rot_horizontal  =  360, 

xrange  =  [-3 .4,3.4] , 

yrange  =  [-3 .4,3.4] , 

zrange  =  [-1.4, 1.4], 

xtics  =  none, 

ytics  =  none, 

ztics  =  none, 

axis_3d  =  false, 

surface_hide  =  true, 

parametric_surf ace ( (2+cos (u/2) *sin(v) -sin (u/2) *sin(2*v) ) *cos (u) , 
(2+cos (u/2) *sin(v) -sin (u/2) *sin(2*v) ) *sin(u) , 
sin (u/2) *sin(v)  +  cos(u/2)*sin(2*v) , 
u,  -°/0pi,  36O*°/opi/18O-0/opi,  v,  0,  2*%pi)  )$ 
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En  el  gráfico  siguiente  se  combinan  varios  objetos:  una  superficie  explícita,  dos  curvas 
paramétricas  y  dos  etiquetas,  cuyo  aspecto  es  el  mostrado  en  el  apartado  e)  de  la  Figura  13. 

(°/0i8)  draw3d  ( 

color  =  green, 

explicit(exp(sin(x)+cos(x~2) ) ,x,-3,3,y,-3,3) , 
color  =  blue, 

parametric(cos(5*u)~2,sin(7*u) ,u-2,u,0,2) , 

color  =  brown, 

line_width  =  2, 

parametric(t~2,sin(t) ,2+t,t,0,2) , 

surface_hide  =  true, 

title  =  "Surface  &  curves", 

color  =  red, 

label( ["UP" ,-2,0,3] ) , 

label( ["DOWN" ,2, 0,-3] ) , 

rot_horizontal  =  10, 

rot_vertical  =  84, 

terminal  =  eps  )$ 

En  el  siguiente  ejemplo  hacemos  una  proyección  esférica  del  hemisferio  sur,  que  se  ve 
en  el  apartado  f)  de  la  Figura  13.  El  paquete  worldmap  carga  en  memoria  las  coordenadas 
latitud-longitud  de  las  líneas  fronterizas  y  costeras  de  todos  los  países  del  mundo,  así  como 
algunas  funciones  necesarias  para  su  manipulación  y  procesamiento, 

(°/0i9)  load  (worldmap)  $ 

(°/0ilO)  draw3d( 

surface_hide  =  true, 
rot_horizontal  =  60, 
rot_vertical  =  131, 
color  =  cyan, 

parametric_surf ace ( 

eos (phi) *cos (theta) , 
cos(phi)*sin(theta) , 
sin(phi) , 
theta ,  -°/0pi ,  %pi , 
phi , -%pi/2 , %pi/2) , 
color  =  red, 

geomap( [South_America, Africa, Australia] , 

[spherical_projection, 0,0,0, 1] ) , 
color  =  blue, 

geomap( [South_America, Africa, Australia] , 

[cylindrical_projection,0,0,0,l,2] ) , 
terminal  =  eps)$ 
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9) 

Figura  13:  Gráficos  generados  con  el  paquete  draw:  a)  y  b)  con  draw2d;  c),  d),  e)  y  f)  con 
draw3d;  g).  un  gráfico  múltiple. 
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Además  de  los  objetos  gráficos  tridimensionales  ya  vistos,  también  se  hayan  definidos 
points,  vector  y,  a  partir  de  la  versión  5.14  de  Maxima,  implicit. 

También  es  posible  generar  múltiples  gráficos  en  un  mismo  fichero  o  hacer  animaciones 
en  formato  GIF  o  en  la  ventana  gráfica.  Para  ver  más  ejemplos  de  gráficos  generados  con 
el  paquete  draw,  se  recomienda  acceder  a  la  dirección 

http://www.telefonica.net/web2/biomates/maxima/gpdraw 

o  consultar  el  sistema  de  ayuda  de  Maxima. 

Ya  como  colofón,  un  ejemplo  de  gráfico  múltiple  en  el  que  se  muestra  también  cómo 
dar  sombreado  a  las  superficies  tridimensionales.  El  resultado  en  el  apartado  g)  de  la 
Figura  13. 

CY.ill)  escenal : 

gr3d(surf ace_hide  =  true, 
enhanced3d  =  true, 
palette  =  gray, 

explicit (sin(sqrt (x~2+y~2) ) ,x, -5, 5, y, -5,5))$ 

(°/0il2)  escena2: 

gr3d(surf ace_hide  =  true, 
enhanced3d  =  true , 
palette  =  gray, 

user_preamble  =  "set  pm3d  map" , 
explicit (sin(sqrt (x~2+y~2) ) ,x,-5 ,5 ,y ,-5 ,5) ) $ 

(°/0il3)  draw( 

columns  =  2, 
terminal  =  eps_color, 
eps_width  =  20, 
eps_height  =  8, 
escenal,  escena2)$ 

Se  comienza  definiendo  las  dos  escenas  a  representar,  la  función  explícita  y  el  gráfico 
de  densidades,  en  sendos  objetos  gr3d,  ya  que  ambos  son  de  naturaleza  tridimensional. 
Estos  objetos  se  pasan  luego  a  la  función  draw  como  argumentos,  junto  con  algunas  op¬ 
ciones  globales.  La  función  draw  es  realmente  la  que  ha  generado  también  los  gráficos 
anteriores,  siendo  draw2d  y  draw3d  sinónimos  de  draw(gr2d( .  .  . )  )  y  draw(gr3d(  ...)), 
respectivamente. 
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8.  Programando  en  Maxima 

8.1.  Programación  a  nivel  de  Maxima 

Esencialmente,  programar  consiste  en  escribir  secuencialmente  un  grupo  de  sentencias 
sintácticamente  correctas  que  el  intérprete  pueda  leer  y  luego  ejecutar;  la  manera  más 
sencilla  de  empaquetar  varias  sentencias  es  mediante  paréntesis,  siendo  el  resultado  del 
programa  la  salida  de  la  última  sentencia: 

(°/„il)  (a: 3,  b:6,  a+b) ; 

(%ol)  9 

(°/.i2)  a; 

( %o2)  3 

(°/.i3)  b; 

( %o3)  6 

Como  se  ve  en  las  salidas  %o2  y  °/,o3,  este  método  conlleva  un  peligro,  que  consiste  en 
que  podemos  alterar  desapercibidamente  valores  de  variables  que  quizás  se  estén  utilizando 
en  otras  partes  de  la  sesión  actual.  La  solución  pasa  por  declarar  variables  localmente,  cuya 
existencia  no  se  extiende  más  allá  de  la  duración  del  programa,  mediante  el  uso  de  bloques; 
el  siguiente  ejemplo  muestra  el  mismo  cálculo  anterior  declarando  cyd  locales,  además 
se  ve  cómo  es  posible  asignarles  valores  a  las  variables  en  el  momento  de  crearlas: 

(°/0i4)  block(  [c,d:6]  , 
c  :3, 
c+d  ) ; 


(%o4)  9 

(°/.i5)  c; 

(%o5)  c 

(°/.i6)  d; 

( %o6)  d 

A  la  hora  de  hacer  un  programa,  lo  habitual  es  empaquetarlo  como  cuerpo  de  una 
función,  de  forma  que  sea  sencillo  utilizar  el  código  escrito  tantas  veces  como  sea  necesario 
sin  más  que  escribir  el  nombre  de  la  función  con  los  argumentos  necesarios;  la  estructura 
de  la  definición  de  una  función  necesita  el  uso  del  operador  :  = 
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f (<argl> , <arg2> ,...): =<expr> 

donde  <expr>  es  una  única  sentencia  o  un  bloque  con  variables  locales;  véanse  los 
siguientes  ejemplos: 

(°/.i7)  loga(x,a)  :=  float(log(x)  /  log(a))  $ 

(0/„i8)  loga(7,  4); 


( %o8)  1 .403677461028802 

(°/0i9)  f  act  (n)  :  =block(  [prod:  1]  , 

for  k:l  thru  n  do  prod:prod*k, 
prod  )$ 

(°/„il0)  f act  (45); 

(%ol0)  119622220865480194561963161495657715064383733760000000000 

En  el  primer  caso  (loga)  se  definen  los  logaritmos  en  base  arbitraria  (Maxima  sólo 
tiene  definidos  los  naturales);  además,  previendo  que  sólo  los  vamos  a  necesitar  en  su 
forma  numérica,  solicitamos  que  nos  los  evalúe  siempre  en  formato  decimal.  En  el  segundo 
caso  (f  act)  hacemos  uso  de  un  bucle  para  calcular  factoriales.  Esta  última  función  podría 
haberse  escrito  recursivamente  mediante  un  sencillo  condicional, 

(Xill)  fact2(n) :=  if  n=l  then  1 

else  n*fact2(n-l)  $ 


(°/0il2)  fact2(45) ; 

(%ol2)  119622220865480194561963161495657715064383733760000000000 

O  más  fácil  todavía, 

(0/.il3)  45!; 

(%ol3)  119622220865480194561963161495657715064383733760000000000 

Acabamos  de  ver  dos  estructuras  de  control  de  flujo  comunes  a  todos  los  lenguajes  de 
programación,  las  sentencias  if-then-else  y  los  bucles  for.  En  cuanto  a  las  primeras, 
puede  ser  de  utilidad  la  siguiente  tabla  de  operadores  relaciónales 
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= 

...igual  que... 

# 

...diferente  de... 

> 

...mayor  que... 

< 

...menor  que... 

>= 

...mayor  o  igual  que... 

<= 

...menor  o  igual  que... 

He  aquí  un  ejemplo  de  uso. 

(°/0il4)  makelist(is(log(x)<x-l) ,  x,  [0 . 9 , 1 , 1 . 1]  ) ; 

( %ol4)  [true,  false,  true] 

Los  operadores  lógicos  que  frecuentemente  se  utilizarán  con  las  relaciones  anteriores 
son  and,  or  y  not,  con  sus  significados  obvios 

p  q  p  and  q  p  or  q  not  p 

false  false  false  false  true 

false  true  false  true  true 

true  false  false  true  false 

true  true  true  true  false 

(°/0il5)  if  sin(4)<  9/10  and  2~2  =  4  then  1  else  -1  ; 

( %ol5)  1 

(°/0il6)  if  sin(4)<  -9/10  and  2~2  =  4  then  1  else  -1  ; 

(%ol6)  -1 

Se  ilustra  en  el  ejemplo  anterior  (°/0il4)  cómo  se  evalúa  con  la  función  is  el  valor  de 
verdad  de  una  expresión  relacional  ante  la  ausencia  de  un  condicional  o  de  un  operador 
lógico;  la  siguiente  secuencia  aclara  algo  más  las  cosas: 

(0/.il7)  sin(4)<  9/10  ; 

9 

( %ol7)  sin  4  <  — 

(°/0il8)  is(sin(4)<  9/10); 

( %ol8)  true 

Pero  ante  la  presencia  de  un  operador  lógico  o  una  sentencia  condicional,  i  s  ya  no  es 
necesario: 
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(0/„il9)  sin(4)<  9/10  and  2~2  =  4; 

( %ol9)  true 

(°/0i20)  if  sin(4)<  9/10  then  1; 

( %o20)  1 

En  cuanto  a  los  bucles,  for  es  muy  versátil;  tiene  las  siguientes  variantes: 

for  <var>:<vall>  step  <val2>  thru  <val3>  do  <expr> 
for  <var>:<vall>  step  <val2>  while  <cond>  do  <expr> 
for  <var>:<vall>  step  <val2>  unless  <cond>  do  <expr> 

Cuando  el  incremento  de  la  variable  es  la  unidad,  se  puede  obviar  la  parte  de  la 
sentencia  relativa  a  step,  dando  lugar  a  los  esquemas 

for  <var>:<vall>  thru  <val3>  do  <expr> 
for  <var>:<vall>  while  <cond>  do  <expr> 
for  <var>:<vall>  unless  <cond>  do  <expr> 

Cuando  no  sea  necesaria  la  presencia  de  una  variable  de  recuento  de  iteraciones,  tam¬ 
bién  se  podrá  prescindir  de  los  for,  como  en 

while  <cond>  do  <expr> 
unless  <cond>  do  <expr> 

Algunos  ejemplos  que  se  explican  por  sí  solos: 

(°/0i21)  for  z :  -5  while  z+2<0  do  print(z)  ; 

-  5 

-  4 

-  3 

( %o21)  done 

(°/0i22)  for  z:-5  unless  z+2>0  do  print(z)  ; 

-  5 

-  4 

-  3 

-  2 


(%o22) 


done 
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(°/0i23)  for  contri  thru  3  step  0.5  do  (var:  cont~3,  print(var))  ; 

1 

3.375 

8.0 

15.625 

27.0 

( %o23)  done 

(°/0i24)  [z ,  cont ,  var]  ; 

( %o24)  [z,  cont ,  27.0] 

(°/0i25)  while  random(20)  <  15  do  print("qqq")  ; 

qqq 

qqq 

qqq 

( %o25)  done 

Véase  en  el  resultado  °/0o24  cómo  las  variables  z  y  cont  no  quedan  con  valor  asignado, 
mientras  que  var  sí;  la  razón  es  que  tanto  z  como  cont  son  locales  en  sus  respectivos 
bucles,  expirando  cuando  éste  termina;  esto  significa,  por  ejemplo,  que  no  sería  necesario 
declararlas  locales  dentro  de  un  bloque. 

Otra  variante  de  for,  sin  duda  inspirada  en  la  propia  sintaxis  de  Lisp,  es  cuando  el 
contador  recorre  los  valores  de  una  lista;  su  forma  es 

for  <var>  in  <lista>  do  <expr> 


y  un  ejemplo: 

(°/0i26)  for  z  in  [sin(x) ,  exp(x) ,  x~(3/4)]  do  print(diff  (z,x))  $ 
cos(x) 
x 

°/„e 


3 


1/4 


4  x 
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Cuando  una  función  debe  admitir  un  número  indeterminado  de  argumentos,  se  colo¬ 
carán  primero  los  que  tengan  carácter  obligatorio,  encerrando  el  último  entre  corchetes 
para  indicar  que  a  partir  de  ahí  el  número  de  argumentos  puede  ser  arbitrario.  La  siguiente 
función  suma  y  resta  alternativamente  las  potencias  n-ésimas  de  sus  argumentos,  siendo 
el  exponente  su  primer  argumento, 

(°/0i27)  sumdif  (n,  [x] )  :  =  x~n  .  makelist  ( (-1)  ~  (k+1)  ,  k,  1,  length(x))  $ 

(°/0i28)  sumdif (7,a,b,c,d,e,f ,g) ; 

( %o28)  g7  -  f +  e7  -d7 +  c7  -b7  +  a7 

(°/0i29)  sumdif  (°/0pi  ,u+v) ; 

(%o29)  (v  +  uf 

En  la  entrada  °/„i27,  dentro  del  cuerpo  de  la  función,  la  variable  x  es  la  lista  de  los 
argumentos  que  podríamos  llamar  restantes ;  como  tal  lista,  la  función  length  devuelve  el 
número  de  elementos  que  contiene.  Recuérdese  también  que  una  lista  elevada  a  un  expo¬ 
nente  devuelve  otra  lista  con  los  elementos  de  la  anterior  elevados  a  ese  mismo  exponente. 
Por  último,  el  operador  .  calcula  el  producto  escalar. 

En  otras  ocasiones,  puede  ser  de  interés  pasar  como  argumento  una  función,  en  lugar 
de  un  valor  numérico  o  una  expresión.  Esto  no  tiene  problema  para  Maxima,  pero  debe 
hacerse  con  cierto  cuidado.  A  continuación  se  define  una  función  de  Maxima  que  toma 
como  argumento  una  función  y  la  aplica  a  los  cinco  primeros  números  enteros  positivos, 

(°/„i30)  fun_a(G)  :=  map(G,  [1,2, 3, 4,5])  $ 

Pretendemos  ahora  aplicar  a  estos  cinco  elementos  la  función  seno,  lo  cual  no  da  ningún 
problema, 

(°/0i31)  fun_a(sin) ; 

( %o31)  [sin  1 ,  sin  2,  sin  3,  sin  4,  sin  5] 

Pero  el  problema  lo  vamos  a  tener  cuando  queramos  aplicar  una  función  algo  más 
compleja, 

(°/0i32)  fun_a(sin(x)+cos(x))  ; 

Improper  ñame  or  valué  in  functional  position: 
sin(x)  +  cos(x) 

#0:  fun_a(g=sin(x)+cos (x) ) 

—  an  error.  To  debug  this  try  debugmode(true) ; 
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lo  cual  se  debe  al  hecho  de  que  la  función  map  no  reconoce  el  argumento  como  función.  En 
tal  caso  lo  apropiado  es  definir  la  función  objetivo  separadamente,  lo  cual  podemos  hacer 
como  función  ordinaria  o  como  función  lambda: 

(°/0i33)  sc(z):=  sin(z)  +  cos(z)  $ 

(°/0i34)  fun_a(sc) ; 

( %o34)  [sin  1  +  eos  1 ,  sin  2  +  eos  2,  sin  3  +  eos  3,  sin  4  +  eos  4,  sin  5  +  eos  5] 

(°/0i35)  fun_a(  lambda([z],  sin(z)+cos (z) )  ); 

( %o35)  [sin  1  +  eos  1 ,  sin  2  +  eos  2,  sin  3  +  eos  3,  sin  4  +  eos  4,  sin  5  +  eos  5] 

Cuando  se  hace  un  programa,  lo  habitual  es  escribir  el  código  con  un  editor  de  texto  y 
almacenarlo  en  un  archivo  con  extensión  mac.  Una  vez  dentro  de  Maxima,  la  instrucción 

load ( "ruta/fichero . mac " ) $ 

leerá  las  funciones  escritas  en  él  y  las  cargará  en  la  memoria,  listas  para  ser  utilizadas. 

Si  alguna  de  las  funciones  definidas  en  el  fichero  fichero  .mac  contiene  algún  error,  de¬ 
volviéndonos  un  resultado  incorrecto,  Maxima  dispone  del  modo  de  ejecución  de  depurado 
(debugmode),  que  ejecuta  el  código  paso  a  paso,  pudiendo  el  programador  chequear  inter¬ 
activamente  el  estado  de  las  variables  o  asignarles  valores  arbitrariamente.  Supongamos  el 
siguiente  contenido  de  cierto  fichero  prueba. mac 


foo(y)  := 

block  ( [u: y~2] , 
u:  u+3, 
u:  u~2, 
u); 


Ahora,  la  siguiente  sesión  nos  permite  analizar  los  valores  de  las  variables  en  cierto 
punto  de  la  ejecución  de  la  función: 

(°/0i26)  load( "prueba. mac" ) ;  /*  cargamos  fichero  */ 

( %o26)  prueba.mac 

(°/0i27)  :break  foo  3  /*  declaro  punto  de  ruptura  al  final  de  línea  3*/ 

Turning  on  debugging  debugmode (true) 

Bkpt  0  for  foo  (in  prueba.mac  line  4) 
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(°/0i27)  foo(2); 

Bkpt  O : (prueba . mac  4) 
prueba . mac : 4 : : 

(dbm:l)  u; 

7 

(dbm: 1)  y; 

2 

(dbm: 1)  u:  1000; 

1000 

(dbm:l)  : continué 


/*  llamada  a  función  */ 

/*  se  detiene  al  comienzo  de  la  línea  4  */ 
/*  ¿valor  de  u?  */ 

/*  ¿valor  de  y?  */ 

/*  cambio  valor  de  u  */ 

/*  continúa  ejecución  */ 


( %o27)  1000000 

8.2.  Programación  a  nivel  de  Lisp 

A  veces  se  presentan  circunstancias  en  las  que  es  mejor  hacer  programas  para  Maxima 
directamente  en  Lisp;  tal  es  el  caso  de  las  funciones  que  no  requieran  de  mucho  procesa¬ 
miento  simbólico,  como  funciones  de  cálculo  numérico  o  la  creación  de  nuevas  estructuras 
de  datos. 

Sin  dar  por  hecho  que  el  lector  conoce  el  lenguaje  de  programación  Lisp,  nos  limitare¬ 
mos  a  dar  unos  breves  esbozos.  Empecemos  por  ver  cómo  almacena  Maxima  internamente 
algunas  expresiones: 

(°/.il)  3/4; 

(%ol)  l 

(°/0i2)  :lisp  $%o  1 
((RAT  SIMP)  3  4) 

(°/0i2)  :lisp  ($num  $°/„ol) 

3 

La  expresión  |  se  almacena  internamente  como  la  lista  ( (RAT  SIMP)  3  4) .  Una  vez 
se  ha  entrado  en  modo  Lisp  mediante  :  lisp,  le  hemos  pedido  a  Maxima  que  nos  devuelva 
la  expresión  %ol;  la  razón  por  la  que  se  antepone  el  símbolo  de  dólar  es  que  desde  Lisp,  los 
objetos  de  Maxima,  tanto  variables  como  nombres  de  funciones,  llevan  este  prefijo.  En  la 
segunda  instrucción  que  introducimos  a  nivel  Lisp,  le  indicamos  que  aplique  a  la  fracción 
la  función  de  Maxima  num,  que  devuelve  el  numerador  de  una  fracción. 

Del  mismo  modo  que  desde  el  nivel  de  Lisp  ejecutamos  una  instrucción  de  Maxima, 
al  nivel  de  Maxima  también  se  pueden  ejecutar  funciones  de  Lisp.  Como  ejemplo,  veamos 
el  comportamiento  de  las  dos  funciones  de  redondeo,  la  definida  para  Maxima  y  la  propia 
del  lenguaje  Lisp: 
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(°/0i3)  round  (°/„pi) ; 

( %o3)  3 

(°/0i4)  ?round(°/0pi)  ; 

Maxima  encountered  a  Lisp  error: 

ROUND:  $%PI  is  not  a  real  number 
Automatically  continuing. 

To  reenable  the  Lisp  debugger  set  *debugger-hook*  to  nil. 

En  el  primer  caso  llamamos  a  la  función  de  redondeo  de  Maxima  y  nos  devuelve  un 
resultado  que  reconocemos  como  correcto.  En  el  segundo  caso,  al  utilizar  el  prefijo  ?, 
estamos  haciendo  una  llamada  a  la  función  round  de  Lisp  y  obtenemos  un  error;  la  razón 
es  que  esta  función  de  Lisp  sólo  reconoce  números  como  argumentos.  En  otras  palabras, 
a  nivel  de  Maxima  se  ha  redefinido  round  para  que  reconozca  expresiones  simbólicas. 

A  efectos  de  saciar  nuestra  curiosidad,  veamos  cómo  se  almacenan  internamente  otras 
expresiones: 

(°/0i5)  x+y*2 ; 


(%o5)  2  y  +  x 

C/.Í6)  [1,2]; 


( %o6)  [1,2] 

(°/0i7)  :  lisp  $%o5 

( (MPLUS  SIMP)  $X  ( (MTIMES  SIMP)  2  $Y) ) 

(°/0i7)  :  lisp  $%o6 
( (MLIST  SIMP)  1  2) 

Otra  forma  de  mostrar  la  representación  interna  de  expresiones  de  Maxima  es  invo¬ 
cando  la  función  print  de  Lisp.  El  siguiente  ejemplo  nos  muestra  que  Maxima  almacena 
internamente  una  expresión  negativa  como  un  producto  por  -1. 

(°/0i8)  ?print(-a)$ 


((MTIMES  SIMP)  -1  $A) 
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Para  terminar,  un  ejemplo  de  cómo  definir  a  nivel  de  Lisp  una  función  que  cambia  de 
signo  una  expresión  cualquiera  que  le  pasemos  como  argumento: 

(°/0i8)  :lisp  (defun  $opuesto  (x)  (list  ’ (mtimes)  -1  x)) 

$0PUEST0 

(°/0i9)  opuesto  (a+b)  ; 

(%o9)  —b  —  a 

Se  ha  hecho  un  gran  esfuerzo  para  que  Maxima  sea  lo  más  universal  posible  en  cuanto 
a  entornos  Common  Lisp,  de  tal  manera  que  en  la  actualidad  Maxima  funciona  perfecta¬ 
mente  con  los  entornos  libres  clisp,  gcl,  cmucl  o  sbcl. 

También  es  posible  llamar  desde  el  nivel  de  Lisp  funciones  que  previamente  hayan  sido 
definidas  en  Maxima,  para  lo  cual  será  necesario  utilizar  la  función  Lisp  mfuncall. 

(°/0ilO)  cubo(x):=  x~3$ 

(lili)  : lisp  (mfuncall  ’ $cubo  3) 

27 


Con  :lisp  sólo  podemos  ejecutar  una  única  sentencia  Lisp;  pero  si  lo  que  queremos 
es  abrir  una  sesión  Lisp  completa,  haremos  uso  de  la  función  to_lisp,  terminada  la  cual 
volveremos  a  Maxima  ejecutando  (to-maxima). 

(°/0il2)  to_lisp()$ 

Type  (to-maxima)  to  restart,  ($quit)  to  quit  Maxima. 

MAXIMA>  (+  #c(3  4)  #c(5  6/8))  ;  suma  complejos 
#C(8  19/4) 

MAXIMA>  (if  (<  2  0) 

"2  es  negativo" 

"2  no  es  negativo")  ;  ejemplo  de  condicional 
"2  no  es  negativo" 

MAXIMA>  (to-maxima) 

Returning  to  Maxima 
(°/.il3) 


índice  alfabético 

;  (punto  y  coma),  7 
,  (coma),  43 
$,  9 
’,  9 
”,  10 
%,  11 
%in,  10 
%on,  10 
%e,  11 
%i,  11 
%pi,  11 
%gamma,  11 
%phi,  11 

+,  18,  25,  26,  59,  87 
%rn,  53 
-,  26,  59 
*,  25,  26,  59,  87 
18 

*,  18 
.,  59 

/,  26,  59 
/,  18 
"  ,  26,  59 
*  ,  18 
**,  18 
.,  26 
,  58 
.,  87 
=,  115 
#,  115 
>,  115 
<,  115 
>=,  115 
=<,  115 
:=,  71 
:=,  56 
:lisp,  22 
23 

/*-*/>  19 


[...],  23 

",  89 

:  (operador  de  asignación),  8 
?  (operador  interrogación),  11 
??  (operador  doble  interrogación),  12 

abs,  22,  69 

acos,  69 

acosh,  69 

addcol,  57 

addrow,  57 

adjacency_matrix,  38 

adjoion,  35 

airy,  71 

algsys,  52 

allroots,  54 

and,  116 

append,  25 

apply,  25 

applyl,  65,  67 

applybl,  66 

array,  28 

asin,  69 

asinh,  69 

assume,  78,  90 

atan,  69 

atan2,  69 

atanh,  69 

atvalue,  82 

B,  19 

barsplot,  97 
bc2,  81 
bessel,  71 
beta,  69 
bfloat,  18 
binomial,  69 
block, 114 
boxplot,  97 

cardinality,  36 
cartesian.product,  36 
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cdfmormal,  90 
central  mornent.  95 
charat,  31 
charfun2,  103 
charlist,  31 
charpoly,  60 
chromatic_number,  38 
col,  58 
color,  109 

complement_graph,  39 
complete_graph,  40 
concat,  31 
conjúgate,  22 
cons,  24 

continuous  f'req,  95 

cor,  97 

eos,  69 

cosh,  69 

cot,  69 

cov,  97 

create_graph,  37 
esc,  69 
cspline,  105 
curl,  88 
cv,  95 

eyele  grapli,  40 

debugmode,  120 
declare,  68 
defrule,  65,  67 
del,  75 
delete,  23 
denom,  47 
dependencies,  74 
depends,  74 
descriptive,  93 
desolve,  82 
determinant,  60 
detout,  60 
diag,  62 
diagmatrix,  57 
diameter,  39 
diff,  73,  80 
diffeq,  84 


discrete  jreq,  95 
disjoin,  35 
disp Jordán,  62 
distrib,  90,  92 
div,  88 
divisors,  20 

dodecahedron_graph,  40 
draw,  101,  108 
draw2d,  108 
draw3d,  108 
draw_graph,  38 
dynamics,  86 

E,  19 

eigenvalues,  60 
eigenvectors,  61 
elementp,  35 
elimínate,  53 
ellipse,  109 
elliptic,  71 
emptyp,  35 
endeons,  24 
entermatrix,  55 
erf,  69,  71 
ev,  13,  43,  88 
evaLstring,  34 
evenp,  19 
exp,  69 

expand,  43,  44,  46 
explicit,  109 
express,  88 

factor,  19,  44-46 
fast_linsolve,  54 
file_search,  94 
fillarray,  28 
find  root,  54 
first,  23 
fixnum,  28 
flatten,  27 
float,  15,  18 
flonum,  28 
for,  117 
fpprec,  18 


ÍNDICE  ALFABÉTICO 


fpprintprec,  97 
fullratsimp,  46 

gamma,  69 
gcd,  45 
genfact,  69 
genmatrix,  56 
geomap,  110 
globaLvariances,  97 
grad,  87 
gradef,  76 
graphs,  37 

halfangles,  49 
hankel,  64 
hessian,  64 
hilbert_matrix,  64 
histogram,  97 

i,  20 
icl ,  80 
ic2,  81 
ident,  57 
if,  117 
ilt,  79 
image,  109 
imagpart,  21 
implicit,  110 
in,  118 
inchar,  10 
inf,  72 

inference_result,  99 
infix,  68 
intégrate,  77 
interpol,  103 
intersection,  36 
invert,  60 
is,  68,  116 
is  cormected,  39 
is_planar,  39 
Ítems  Jnference,  100 

jordan,  62 

key,  109 


km,  9 

kurtosis,  95 
kurtosis_gumbel,  91 

label,  109 
lagrange,  104 
lambda,  26 
laplace,  79 
laplacian,  88 
last,  23 
lcm,  45 
length,  23 
lhs,  52 
limit,  72 
line  type,  109 
line_width,  109 
linearalgebra,  65 
linearinterpol,  103 
linsolve,  54 
linsolve_by_lu,  54,  63 
lisp,  121 

list_correlations,  97 
listarith,  27 
listarray,  28 
listify,  35 
listp,  23 
load,  13 
log,  69 

logconcoeffp,  48 
logcontract,  48 
logexpand,  47 
lsquares,  100 
lsquares.estimates,  101 
lsquares_residual_mse,  101 
lsquares_residuals,  101 

rnake  array,  30 
make_random_state,  91 
makelist,  25,  91 
map,  25,  36 
matchdeclare,  65,  67 
matrix,  15,  55 
matrixmap,  62 
matrixp,  58 


ÍNDICE  ALFABÉTICO 


127 


max,  69 

maxapplydepth,  66 
maxapplyheight,  66 
mean,  95 
mean_weibull,  90 
median,  95 
member,  23 
min,  69 
minf,  72 
minor,  61 
minus,  73 
mnewton,  55 
ModeMatrix,  62 

negative,  77 
noncentral  mornent ,  95 
not,  116 
nticks,  109 
nullspace,  61 
num,  47 
numer,  15,  96 

oddp,  19 
ode2,  80 
or,  116 
outchar,  10 

parametric,  109 
parse  string,  34 
part,  23,  51 
partfrac,  46 
pdf_poisson,  90 
pdiff,  76 
piechart,  97 
plot2d,  107 
plot3d,  107 
plot_options,  107 
plotdf,  83 
plus,  73 
points,  109 
points_joined,  109 
polar,  110 
polarform,  21 
polygon,  110 
positive,  77 


primep,  19 
print  graph,  37 
product,  25 

quad.qag,  78 
quad  qagi,  78 
quad.qags,  79 
quad  qawc,  79 
quad.qawf,  79 
quad.qawo,  79 
quad.qaws,  79 
quadpack,  78 
quantile,  95 
quantile_student_t,  90 
quit,  13 
quotient,  20 

radcan,  46 

random,  91 

random.binomial,  92 

random  continuous  uniform,  100 

random.digraph,  41 

random_exp,  92 

random_graph,  41 

randomJiypergeometric,  92 

random  poisson,  92 

random.weibull,  92 

range,  95 

rank,  61 

ratsimp,  46 

readdist,  93 

readunatrix,  93 

realonly,  52 

realpart,  21 

realroots,  54 

rectangle,  110 

rectform,  21 

remainder,  20 

remove,  75 

rest,  23 

reverse,  23 

rhs,  52 

rk,  86 

row,  58 


ÍNDICE  ALFABÉTICO 


128 


rungel,  85 
runge2,  85 

save,  12 
scatterplot,  97 
sdowncase,  32 
sec,  69 
second,  23 
sequal,  32 
set,  35 

set_plot_option,  108 
set_random_state,  91 
setdifference,  36 
setify,  35 
shortest  path,  41 
showtime,  19 
signum,  69 
simp,  45,  68 
sin,  69 
sinh,  69 
skewness,  95 
skewness_rayleigh,  90 
slength,  31 
sol  ve,  51 
sort,  23 
sqrt,  69 
ssearch,  33 
ssubst,  33 
stats,  99 
std,  95 
step,  117 
stringout,  13 
stringp,  31 
submatrix,  56 
subsetp,  36 
subst,  13,  36,  44 
substring,  34 
sum,  25 
supcase,  32 
System,  16 

take Jnference,  99 
tan,  69 
tanh,  69 


taylor,  76 
tellsimp,  68 
tellsimpafter,  68 
terminal,  109 
test_mean,  99,  102 
testmormality,  99,  102 
test_signed_rank,  100 
tex,  16 
third,  23 
thru,  117 
time,  19 
title,  109 
todisp,  123 
toeplitz,  64 
tokens,  31 
transpose,  60 
trigexpand,  48 
trigexpandplus,  48 
trigexpandtimes,  48 
triginverses,  50 
trigrat,  50 
trigsign,  50 
trigsimp,  50 

und,  73 
unión,  36 

uniteigenvectors,  61 
unless,  117 

vandermonde_matrix,  64 
var,  95 

var  binomial,  90 
vect,  87 
vector,  109 
vértices,  41 

while,  117 
worldmap,  111 

zero,  77 
zeromatrix,  57 


